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PREFACE 


This  thesis  represents  the  first  investigation  into  the  compatibility  of  Field 
Emission  (FE)  cathodes  and  Electric  Propulsion  (EP)  systems.  FE  cathodes  have 
demonstrated  1000s  of  hours  of  operation  in  Ultra  High  Vacuum  (UHV)  environments. 
Hall  thrusters  have  demonstrated  1000s  of  hours  of  operation  in  10'6-10'5  Torr  with 
hollow  cathodes.  FE  cathodes  are  being  considered  to  replace  hollow  cathodes  in  Hall 
and  ion  thrusters  because  they  do  not  require  propellant,  operate  cold,  and  operate  at 
much  lower  power  levels.  These  features  are  especially  important  as  EP  systems  are 
shrinking  in  size  and  power  levels  to  meso-  and  microscales.  The  effects  of  hostile 
thruster  environments  on  the  performance  of  FE  cathodes  was  investigated  in  this  thesis  . 
work.  Issues  of  contamination,  lifetime,  and  space-charge  limited  electron  emission  were 
addressed.  Experiments  were  performed  to  determine  if  silicon  and  molybdenum  Field 
Emission  Array  (FEA)  cathodes  could  operate  in  xenon  environments  without 
performance  degradation,  what  caused  the  performance  degradation,  and  how  the 
degradation  rate  could  be  controlled.  Models  were  developed  to  predict  performance 
degradation  rates  of  Si  and  Mo  FEA  cathodes  in  xenon  environments.  They  were  used  to 
determine  cathode  lifetimes  and  suggest  cathode  configurations  that  will  satisfy  lifetime 
and  current  requirements.  A  sheath  model  was  also  developed  to  predict  space-charge 
current  limits  in  operating  configurations  where  a  plasma  provides  a  virtual  anode  to 
collect  the  emitted  current.  This  model  was  also  used  to  investigate  the  space-charge 
limits  on  FE  cathodes  in  Low  Earth  Orbit  (LEO)  on  a  tether  system.  The  sheath  and 
performance  models  were  used  to  optimize  cathode  configurations  for  thruster  and  tether 
applications. 


vii 


TABLE  OF  CONTENTS 


ACKNOWLEDGMENTS . IV 

PREFACE . VH 

LIST  OF  FIGURES . XI 

LIST  OF  TABLES . XVI 

LIST  OF  APPENDICES . XVIII 

CHAPTERS . 1 

1  INTRODUCTION . 1 

1.1  Electric  Propulsion . 1 

1.1.1  Role  of  Electric  Propulsion . 1 

1.1.2  Hall  Thrusters . 3 

1.1.3  Ion  Thrusters . 5 

1.1.4  Oxygen  Thrusters . 6 

1.1.5  Meso-  and  Microscale  Thrusters . 7 

1.1.6  Conventional  Electric  Propulsion  Cathode 

Technology . 8 

1.1.7  Advanced  Cathode  Requirements . 9 

1.2  Space  Tethers . 9 

1.2.1  Electrodynamic  Tethers . 10 

1.1.2  Cathode  Requirements . 10 

1.3  Field  Emission(FE)  Cathodes . 1 1 

1.1.1  Field  Electron  Emission . 12 

1.1.2  State-of-the-Art  FE  Cathode  Technology . 13 

1.4  Challenges . 15 

1.4.1  Cathode  Lifetime . 15 

1.1.2  Space-Charge  Limited  Emission . 18 

1.5  Role  of  this  Thesis . 21 

2  CATHODE  ENVIRONMENT  CHARACTERIZATION . 25 

2.1  Experimental  Apparatus . 25 

2.1.1  Test  Facility . 25 

2.1.2  Hall  Thruster-PEPL-70 . 27 

2.1.3  Diagnostics . 28 


vm 


2.2  Experimental  Methods  and  Results . 29 

2.2. 1  Thruster  Performance . . . 29 

2.2.2  Xenon  Pressure  and  Ion  Current  Density  in 

Cathode  Region . 32 

3  FIELD  EMISSION  CATHODE  EXPERIMENTS . 34 

3.1  Carbide  Single  Tip  Cathode  Experiments . 35 

3.1.1  Cathodes . 35 

3.1.2  Experimental  Apparatus . 36 

3.1.3  Experimental  Methods  and  Results . 37 

3.2  Mo  and  Si  Field  Emission  Array  Cathode  Experiments . 42 

3.2.1  Cathodes . 43 

3.2.2  Experimental  Apparatus . 46 

3.2.3  Si  FEA  Cathode  Experimental  Methods  and 

Results . 48 

3.2.4  Mo  FEA  Cathode  Experimental  Methods  and 

Results . 33 

3.3  Carbon  Film  FE  Cathode  Experiments . 59 

3.3.1  Carbon  Film  FE  Cathode . 60 

3.3.2  Experimental  Apparatus . .  61 

3.3.3  Experimental  Methods  and  Results . 62 

4  MICROTIP  CATHODE  PERFORMANCE  MODELING . 67 

4.1  FEA  Cathode  Performance  Degradation  Model . 67 

4.1.1  Field  Electron  Emission  Model . 68 

4.1.2  Tip  Sputtering  Model . 72 

4.1.3  Sputter  Yield  Model . 76 

4.2  Modeling-Results . 80 

4.2.1  Mo  SRI  International  Cathode . 81 

4.2.2  Si  MCNC  FE  Cathodes . 87 

4.3  Application  of  Modeling  and  Experimental  Results . 90 

5  SPACE-CHARGE  LIMITED  ELECTRON  EMISSION  FROM 

FIELD  EMISSION  CATHODES  INTO  A  PLASMA . 96 

5.1  Sheath  Model . 99 

5.1.1  Planar  Sheath  Geometry . 99 

5.1.2  Cylindrical  Geometry . 104 

5.1.3  Spherical  Geometry . 107 

5.2  Modeling  Results . 108 

6  CONCLUSIONS . 122 

6.1  Summary  of  Major  Results . 122 

6.2  Recommended  Cathode  Designs . 126 


IX 


6.3  Recommendations  for  Future  Research . 128 

6.4  Other  Applications . 129 

APPENDICES . 131 

APPENDIX  A . 132 

APPENDIX  . . 137 

APPENDIX  . . 141 

APPENDIX  . . . — . 149 

BIBLIOGRAPHY . 155 


X 


LIST  OF  FIGURES 


Figure  1.1  Cross-section  of  a  Hall  thruster  with  a  hollow  cathode . 5 

Figure  1.2  Cross-section  of  an  ion  thruster  (courtesy  of  Matt  Domonkos) . 6 

Figure  1.3  A  Hall  thruster  shown  with  hollow  and  field  emission  cathodes . 11 

Figure  1.4  Potential  energy  diagram  for  electrons  at  a  metal  surface,  without  (left)  and 

with  (right)  an  applied  electric  field,  F . 12 

Figure  1.5  Field  emission  tip  and  gate  electrode  configuration  (courtesy  of  Dr.  Kevin 

Jensen) . ^ 

Figure  1.6  Potential  diagram  of  the  cathode,  thruster,  and  plasma  configuration.  Anode,  , 
plasma  ion  beam,  sheath  boundary,  gate  electrode,  and  cathode  potentials  are  shown 

relative  to  ground . 19 

Figure  1.7  Flow  diagram  for  the  field  emission  array  cathode  performance  degradation 

model . *23 

Figure  2.1  Thrust  stand  calibration  data . 27 

Figure  2.2  The  PEPL-70  Hall  thruster  with  a  hollow  cathode,  neutral  particle  flux  (NPF) 

probe,  and  Faraday  Probe . 28 

Figure  2.3  Faraday  probe  schematic  (courtesy  of  James  Haas) . 29 

Figure  2.4  I-V  curves  from  the  PEPL-70  at  different  xenon  volumetric  flow  rates . 31 

Figure  2.5  Thruster  performance  data  with  varying  propellant  flow  rate  and  a  constant 
discharge  voltage,  300V.  The  numbers  on  the  markers  on  the  graph  correspond  to 
the  total  propellant  flow  rate  (seem).  The  table  above  the  graph  shows  the  discharge 
current  also . 31 


xi 


Figure  3.1  A  scanning  electron  microscope  (SEM)  micrograph  of  a  single  carbide 

cathode  in  a  Vogel-type  mount . 36 

Figure  3.2  An  illustration  of  how  the  field  emission  microscope  (FEM)  was  used  to 

determine  the  emitting  area  of  the  tip  from  the  angle,  cr,  which  can  be  calculated 

from  the  size  of  the  emission  pattern,  d,  on  the  screen  at  a  distance  x  from  the 

37 

emitter . Jl 

Figure  3.3  Electrical  schematic  of  the  FEM . 37 

Figure  3.4  Emission  patterns  observed  with  the  FEM  from  fairly  clean  carbide 

39 

cathodes . 

Figure  3.5  Emission  patterns  observed  with  the  FEM  from  damaged  or  contaminated 

39 

carbide . * 

Figure  3.6  Evolution  of  the  tip  emission  pattern  observed  with  the  FEM  a,b)  during  and 

c)  after  the  exposure  to  argon . 41 

Figure  3.7  Fowler-Nordheim  (F-N)  (left)  and  I-V  curves  (right)  obtained  from  the  ZrC 

cathode  before  and  after  exposure  to  argon . 42 

Figure  3.8  A  MCNC  Si  field  emission  array  single  tip  configuration  (courtesy  of 

MCNC) . 44 

Figure  3.9  A  SRI  International  Mo  field  emission  array  cathode  single  tip  and  array.  ...45 

Figure  3.10  The  ultra-high  vacuum  (UHV)  facility  at  the  Jet  Propulsion  Laboratory . 47 

Figure  3.11  The  field  emission  cathode  test  flange  in  UHV  facility . 47 

Figure  3.12  Electrical  schematic  for  field  emission  cathode  testing . 48 

Figure  3.13  I-V  curves  taken  before  and  after  the  Si  cathode  was  exposed  to  Xe . 49 

Figure  3.14  Si  cathode  current  response  to  an  increase  in  Xe  pressure  with  Vg  at  8$  V 

and  Va  at  400  V . 50 

Figure  3.15  Si  cathode  current  response  to  2xl0'5  Torr  of  Xe  with  Vg  at  70  V  and  the  Va 

at  70  V . 51 


Xll 


Figure  3.16  MCNC  Si  cathode  current  response  to  a  1  hour  Xe  exposure  with  Vg  at  50  V 

and  Va  at  70  V . 52 

Figure  3.17  I-V  data  taken  before  and  after  the  Mo  cathode  was  exposed  to  xenon . 53 

Figure  3.18  Mo  cathode  current  response  to  2x1  O'6  Torr  of  xenon  with  Vg  at  65.6  V  and 
Va  at  8C  V . . . 54 


Figure  3.19  The  Mo  cathode  current  response  to  increasing  gate  electrode  voltage  with 

Va  at  100  V,  2xl0'5  Torr  of  Xe . 55 

Figure  3.20  A  comparison  of  the  Mo  cathode  current  response  to  increasing  gate 

electrode  voltage  in  2xl0’9  Torr  and  2xl0"5  Torr  of  Xe . 56 

Figure  3.21  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2x1 0'5  Torr  of 

xenon  for  5  hours  with  Vg  at  50  V  and  Vaat  60  V . : . —  57 

Figure  3.22  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2xl0'5  Torr 

with  varying  Va  and  Vg  at  50  V . 57 

Figure  3.23  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2xl0'5  Torr  for 

2  hours  with  Vg  at  50  V  and  Va  at  100  V . 58 

Figure  3.24  A  FEPET  carbon  film  cathode  on  a  glass  stem . 61 

Figure  3.25  ..I-V  traces  showing  the  performance  of  a  FEPET  carbon  film  cathode  in  a 

UHV  environment . 62 

Figure  3.26  The  cathode  current  response  to  increasing  the  Xe  pressure  up  to  2xl0'6 
Torr  for  4  hours.  The  anode  potential  was  900  V  above  the  potential  of  the  carbon 

film . 64 

Figure  3.27  I-V  curves  taken  in  UHV  before  and  after  the  carbon  film  cathode  was 

exposed  to  xenon . 64 

Figure  3.28  Carbon  film  cathode  turn-on  and  turn-off  I-V  curves  with  an  anode  voltage 
of  900  V  in  2xl0'6  Torr  of  xenon  compared  to  turn-on  I-V  data  in  UHV . 65 


xm 


Figure  4.1  The  actual  cathode  configuration,  and  the  configuration  used  in  the  microtip 
sputtering  model  to  predict  the  potential  field  near  the  microtip  and  electric  field  at 


the  microtip . 74 

Figure  4.2  Significant  radial  positions  used  in  the  model  of  microtip  sputtering . 75 

Figure  4.3  Xe  ionization  cross-sect;ons,  Qj,  from  Brown.74 . 76 


Figure  4.4  Sputter  yields  for  a  too  target  being  bombarded  by  Xe  ions.  Shown  on  the 
graph  are  theoretical  curves  for  two  values  of  E*,  and  experimental  values  measured 


by  Rosenberg  and  Wehner . 79 

Figure  4.5  Sputter  yields  for  a  Si  target  being  bombarded  by  Xe  ions.  Shown  on  the 
graph  are  theoretical  curves  for  three  values  of  E*,  and  experimental  values 

measured  by  Rosenberg  and  Wehner . . . . . 80 

Figure  4.6  Experimental  and  theoretical  Mo  cathode  current  response  to  a  Xe  pressure 

increase  up  to  2xl0'6  Torr  with  Vg  at  65.6  V  and  Va  at  80  V  with  E^  of  39  eV . 82 

Figure  4.7  Mo  cathode  I-V  data  measured  experimentally  before  the  Xe  exposure  and 
predicted  by  Jensen’s  performance  model . 83 


Figure  4.8  Performance  decay  rate  predicted  by  the  performance  degradation  model  for  a 
Mo  cathode  operating  in  a  X*. environment  with  a  pressure  of  2xl0'5  Torr  for  a 

cathode  with  parameters  (I)  shown  in  Table  4.3 . 85 

Figure  4.9  Field  emission  array  cathode  performance  decay  predicted  by  the 

performance  degradation  model  for  a  Mo  cathode  operating  in  2xl0'5  Torr  of  Xe  for 


100  hours  at  various  gate  electrode  voltages  with  initial  cathode  parameters  (I) 

shown  in  Table  4.3 . 86 

Figure  4.10  Measured  and  predicted  Si  cathode  current  response  to  an  increase  in  Xe 

pressure  up  to  2xl0’5  Torr  with  Vg  and  Va  at  70  V . 88 

Figure  4.11  Measured  and  calculated  Si  cathode  performance  degradation  as  the  Xe 

pressure  increases  up  to  7xl0‘6  Torr  with  Vg  at  86  V  and  Va  at  400  V . 89 


xiv 


Figure  5.1  Cathode  sheath  configuration  which  was  modeled.  The  boundary  between  the 

sheath  and  pre-sheath  is  shown  at  . . 97 

Figure  5.2  A  planar  cathode  sheath  scenario  (left)  with  r^r^  and  a  spherical  cathode 

sheath  scenario  (right)  with  . . 98 

Figure  5.3  Results  of  the  cathode  sheath  model  compared  to  the  Prewett  and  Allen  sheath 

model  results . ^ 

Figure  5.4  The  effect  of  Vg  on  for  the  planar  sheath  geometry.  The  numeric 

markers  represent  the  value  of  r|0  for  each  solution . . . 1 1 1 

Figure  5.5  The  effect  of  <t>c  on  for  the  planar  sheath  geometry.  The  numeric  markers 

represent  the  value  of  ri0  for  each  solution . 1 1 1 

Figure  5.6  a,  b  Potential  (a)  and  electric  field  (b)  profiles  for  different  emission  currents 

in  the  cathode  sheath . 

Figure  5.7  Charge  density  profiles  in  a  planar  cathode  sheath  with  r|c=4  and  rig=6 . 1 13 

Figure  B.l  Cathode  Lens  and  Ion  Repeller  (CLAIR)  configuration . 138 

Figure  B.2  Equipotentials  predicted  by  MAGIC  for  the  CLAIR  configuration  described 

in  Figure  B.l  and  Table  B.l . 139 

Figure  B.3  Electron  trajectories  predicted  by  MAGIC  for  rhe  CLAIR  configuration 

described  in  Figure  B.l  and  Table  B.l . 139 


xv 


LIST  OF  TABLES 


Table  2.1  Ion  current  densities  measured  by  the  Faraday  probe  flush  with  the  exit  plane 
of  the  thruster.  Discharge  voltage  and  current,  electromagnet  inner,  1^,  and  outer, 

IM0,  currents,  and  Xe  pressures  are  also  given . 3^ 

Table  2.2  Neutral  Xe  particle  pressures  and  ion  current  densities  measured  by  the  probes 


3  mm  behind  the  physical  exit  plane  of  the  thruster . 32 

Table  4.1  Ion  and  target  material  parameters  used  in  the  sputter  yield  model . 78 

Table  4.2  Cathode  parameters  before  (I)  and  after  (II)  Xe  exposure  obtained  from  I-V 

data  and  predicted  (HI)  by  the  performance  degradation  model . 83 

Table  4.3  Cathode  parameters  before  (I)  and  after  (H)  xenon  exposure  which  are  used  in 

the  performance  decay  model  for  V,  at  100  V  in  2x1  O'5  Torr  of  xenon. . . 85 

Table  4.4  MCNC  Si  cathode  parameters  before  the  exposure  to  xenon  during 

operation . 90 

Table  4.5  SRI  International  #52M  Mo  cathode  parameters  (in  bold)  and  predicted" 


performance  using  different  parameters,  with  and  without  a  CEX  ion  bombardment 

during  100  hours . 9 

Table  4.6  Performance  limitations  of  Si  and  HfC  cathodes  with  and  without  a  charge- 


exchange  (CEX)  ion  population . 94 

Table  5.1  Characteristics  of  the  thruster  and  tether  plasma  environments  investigated 

using  the  cathode  sheath  model . 99 

Table  5.2  Space-charge  limited  currents  estimated  by  the  planar  sheath  model  for 

thruster  and  tether  environments . 109 


xvi 


Table  5.3  Space-charge  limited  currents  predicted  by  the  planar,  cylindrical,  and 

spherical  sheath  models  in  normalized  units . 114 

Table  5.4  Field  emission  current  limitations  predicted  by  the  sheath  model  in 

dimensional  parameters  for  planar,  cylindrical,  and  spherical  geometries  for  the  Hall 
thruster  environment  with  an  external  cathode.  Vg  is  at  30  V  and  Vc  is  at  20  V. ..  1 15 
Table  5.5  Field  emission  current  limitations  predicted  by  the  sheath  model  in  normalized 

parameters  for  a  spherical  sheath  in  a  tether  environment  where  Te  is  0.1  eV . 116 

Table  5.6  Field  emission  current  limitations  predicted  by  the  spherical  sheath  model  for 

the  tether  environment . 116 

Table  5.7  Current  limitations  and  virtual  cathode  potentials  predicted  by  the  Wang 
simulation  model  for  a  field  emission  cathode  in  a  Hall  thruster  environment  with 
rc=  4  Xo,  Vg  at  30  V,  <(>c  at  -20  V,  Te  at  5  eV,  and  neo  at  8xl08/cm3 . 119 

Table  5.8  Current  limitations  and  virtual  cathode  potentials  predicted  by  the  Wang 

simulation  model  for  a  field  emission  cathode  operating  in  a  tether  environment  with 
rc=  Vg  at  30  V,  <J>C  at  -20  V,  Te  at  0.1  eV,  and  neo  at  5x107cm3 . 120 

Table  6.1  Recommended  number  of  cathodes,  operating  voltage,  cathode  material,  and 
emitting  area  required  in  different  applications  if  packing  density  is  108/cm2,  As=  2, 

rg=1000,  and  r,=40 . 127 

Table  6.2  Recommended  number  of  cathodes,  operating  voltage,  cathode  material,  and 
emitting  area  required  in  different  applications  if  packing  density  is  108/cm2,  As=  2, 

rg=1000,  r,=40,  and  CLAIR  is  used . 127 

Table  B.l  Electrode  potentials,  thicknesses,  and  spacing  for  the  CLAIR  configurgation 
simulated  by  MAGIC . 138 


xvn 


LIST  OF  APPENDICES 


Appendix 

A.  Nonemclature . . . 

B.  Cathode  Lens  and  Ion  Repeller  (CLAIR) . 

C.  Fea  Cathode  Performance  Degradation  Model  Macro 

D.  Cathode  Sheath  Model  Macro . 


132 

137 

141 

149 


xvm 


INTRODUCTION 


The  compatibility  of  Field  Emission  (FE)  cathodes  in  Electric  Propulsion  (EP) 
systems  is  explored  in  this  dissertation.  Some  of  the  results  of  this  thesis  work  were  also 
used  to  address  compatibility  issuu.  between  electrodynamic  space  tethers  and  FE 
cathodes.  The  plasma  environments  of  these  systems  place  challenging  demands  on 
cathode  performance.  Field  emission  cathode  performance  in  simulated  environments  of 
interest  is  theoretically  and  experimentally  evaluated.  This  chapter  provides  an 
introduction  to  the  configurations,  operation,  and  applications  of  Hall  and  ion  thrusters, 
electrodynamic  tethers  and  FE  cathodes.  The  motivation  for  using  FE  cathodes  with 
thrusters  and  tethers  and  the  associated  challenges  are  presented.  Previously  obtained 
experimental  results  which  motivated  the  experimental  performance  evaluations 
employed  in  this  thesis  work  is  discussed. 

1.1  Electric  Propulsion 

EP  systems  have  been  under  investigation  since  the  1950s  as  alternatives  to 
chemical  propulsion  systems  on  spacecraft.  They  offer  advantages  including  higher 
specific  impulse,  and,  therefore,  lower  wet  system  mass  for  the  same  mission,  resulting  in 
significantly  reduced  launch  costs.  EP  systems  fall  into  three  categories:  electrothermal, 
electrostatic,  and  electromagnetic.  These  systems  have  been  optimized  for  different 
thrust,  specific  impulse,  and  power  regimes. 

1.1.1  Role  of  Electric  Propulsion 

The  primary  advantage  of  EP  systems  over  chemical  propulsion  systems  is  their 
higher  effective  exhaust  velocity,  ue.  The  amount  of  propellant  required  by  the  rocket 
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engine  to  impart  the  necessary  change  in  velocity,  AV,  to  a  payload  is  determined  by  the 
rocket  equation, 

Ml 

Mf 

Mi  is  the  initial  mass  of  the  spacecraft  including  fuel,  and  Mf  is  the  final  mass  of  the 
spacecraft  after  the  fuel  has  been  spent.  The  amount  of  fuel  required  to  attain  the  desired 
AV  exponentially  increases  with  AV.  Chemical  rockets  are  limited  to  ue  below  ~5000 

m/s.  Particle  exhaust  velocities  can  be  an  order  of  magnitude  higher  in  EP  systems  than 
in  chemical  systems  because  electrical  energy  is  added  to  the  system  to  heat  or  ionize  and 

accelerate  the  propellant  instead  of  relying  on  chemical  energy  of  the  propellant.  The 
rocket  equation  shows  how  much  less  propellant  is  needed  for  the  same  mission  AV  when 

ue  is  significantly  increased  with  the  use  of  electric  propulsion. 

EP  systems  have  been  used  for  North-South  Station-Keeping  (NSSK)  on  satellites 
and  primary  propulsion  for  interplanetary  probes.  They  could  also  be  used  for  orbit 
raising,  comet  sample  return,  and  piloted  missions  to  Mars.  Interplanetary  probes  require 
significant  periods  of  continuous  thrusting  during  a  mission  while  satellites  require  short 
duration  thrusting  periodically  during  the  satellite  lifetime.  Lifetimes  in  excess  of  1000s 
of  hours  are  required  of  EP  systems  for  satellites  or  interplanetary  probes,  and  piloted 
spacecraft.  Power  is  supplied  to  EP  systems  with  solar  (<100  kW)  or  nuclear  sources, 
depending  on  the  amount  of  power  required. 

There  are  several  types  of  EP  systems,  each  with  unique  capabilities  and  roles. 
Currently,  the  most  widely  used  and  simplest  EP  system  is  the  resistojet,  an 
electrothermal  device.  It  is  used  on  communication  satellites  for  NSSK,  attitude  control, 
on-orbit  maneuvering,  and  is  capable  of  running  on  waste  water.  The  MR501  hydrazine 
resistojet  operates  at  350-500  W  to  generate  0.18-0.36  N  of  thrust  at  280-300  s.1  The 
electrothermal  arcjet  is  also  commonly  used  on  communication  satellites,  operating  at  1.8 
kW  using  hydrazine  propellant.  Spacecraft  equipped  with  megawatts  of  power  could 


employ  magnetoplasmadynamic  (MPD)  thrusters,  electromagnetic  devices,  for  sample 
return  missions  and  piloted  missions  to  Mars.  This  system  can  operate  at  60%  efficiency 
on  lithium  propellant.  Field  Emission  Electric  Propulsion  Systems  (FEEP)  are  being 
considered  for  drag  compensation  for  precise  spacecraft  positioning.2  These  systems 
operate  at  10  kV,  ~1  mA,  with  maximum  thrust  approximately  1  mN.  FEEPs  hav-  been 
built  to  operated  at  only  30  jiN.  Hall  thrusters  have  been  used  on  communication 

satellites  that  employ  solar  power.  These  1.35  kW  systems  have  been  optimized  to 
operate  at  -53%  efficiency  on  xenon  to  generate  -83  mN  of  thrust  at  1600  s  specific 
impulse.3  Ion  thrusters  have  been  used  on  communications  satellites  and  interplanetary 
probes.  A  30  cm  ion  thruster  is  currently  being  used  for  primary  propulsion  on  DS-1,  a 
mission  to  fly-by  an  asteroid  and  comet.  A  30  cm  ion  thruster  can  generate  91  mN  of 
thrust  at  a  specific  impulse  of  3380s  and  66%  efficiency  with  xenon  propellant.4  Hall  and 
ion  thrusters  are  both  electrostatic  propulsion  systems.  Constellations  of  hundreds  of 
satellites  for  global  telecommunication  systems,  are  being  developed  which  will  employ 
these  systems. 

1.1.2  Hall  Thrusters 

* 

Development  of  the  Hall  thruster  was  initiated  independently  in  both  the  USSR 
and  the  U.S.  in  the  early  1960s.  Research  efforts  in  the  U.S.  were  more  heavily  focused 
on  the  ion  thruster  by  1970,  while  development  of  the  Hall  thruster  flourished  in  the 
former  USSR.5  Following  the  first  flight  of  a  Hall  thruster,  the  SPT-60,  in  Feb.  1972  on 
the  Meteor  satellite,  approximately  100  Hall  thrusters  have  been  used  on  Russian 
satellites.  Engineers  in  the  U.S.  have  been  working  on  integrating  Russian  Hall  thrusters 
and  U.S.  satellites  since  1992.  Thrusters  developed  for  llkW  communication  satellites 
operate  at  1.35  kW  to  generate  83  mN  of  thrust  and  1600  s  specific  impulse.  Recent 
discussions  about  deploying  constellations  of  —300  satellites  have  included  scenarios  with 
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5-6  kW  of  on-board  power  that  could  support  4.5  kW  propulsion  systems  with  2-6 
systems  on-board.  Since  1993,  there  has  been  a  flurry  of  research  activity  in  the  U.S. 
focusing  on  characterizing  their  performance,  modeling  their  performance,  and 
understanding  the  issues  associated  with  integrating  Hall  thruster  and  satellite  systems. 
High  energy  ions  emitted  by  the  propulsion  systems  can  damage  a  spacecraft  and  plasma 
noise  can  interfere  with  communication  systems.  Hall  thruster  design  development  has 
been  on-going  for  more  than  20  years  to  improve  performance  and  minimize  interference 
with  communication  systems  and  damage  to  spacecraft  surfaces.  These  thrusters  operate 
most  efficiently  on  xenon,  however,  they  have  also  used  Ar,  Kr,  Air,  02 ,  N2,  Bi,  and  Cs 
propellants. 

A  Hall  thruster  developed  at  PEPL  is  shown  in  Figure  1.1.  The  PEPL-70,  based 
on  Russian  designs,  was  developed  as  a  laboratory  model  and  used  only  to  generate  a 
plasma  environment  to  characterize.  Hall  thrusters  employ  an  anode  as  the  positive 
electrode  and  gas  distributor  in  the  annular  ceramic  discharge  chamber.  Multiple 
chambers  are  typically  used  in  the  anode  to  improve  the  azimuthal  uniformity  of  the 
propellant  injection.  The  discharge  chambers  are  usually  composites  of  BN,  Si,  and 
A102.  Hall  thrusters  employ  inner  and  outer  electromagnets  with  poles  at  the  thruster 
exit  plane  to  support  a  radial  magnetic  field  (-300  Gauss  in  the  center  of  the  discharge 
chamber)  with  large  gradients.  The  cathode  is  commonly  positioned  external  to  the 
thruster.  When  the  discharge  voltage  is  applied,  the  anode  and  cathode  support  an  axial 
electric  field.  The  cathode  emits  electrons  to  ionize  the  propellant  and  to  neutralize  the 
ion  beam  emitted  from  the  thruster.  These  thrusters  require  0.2-10  A  of  current, 
depending  on  the  size  of  the  thruster.  Electrons  emitted  by  the  cathode  become  trapped 
in  the  crossing  electric  and  magnetic  fields  in  a  closed  azimuthal  drift.  Xenon  atoms  are 
injected  into  the  drifting  electron  cloud  from  the  anode,  and  are  then  ionized  by  electron 
bombardment.  The  axial  electric  field  accelerates  ions  axially  out  of  the  discharge 
chamber  to  produce  thrust.  In  the  Hall  thruster,  the  ionization  and  acceleration  occur  in 
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the  same  region;  therefore  plasma  quasineutrality  is  maintained  everywhere  in  the 
thruster. 


Figure  1.1  Cross-section  of  a  Hall  thruster  with  a  hollow  cathode. 


1.1.3  Ion  Thrusters 

An  ion  thruster  is  another  form  of  an  electrostatic  EP  system  with  performance 
similar  to  Hall  thrusters.  They  typically  operate  at  a  higher  specific  impulse  and  higher 
efficiency  than  do  Hall  thrusters.  This  thruster  also  ionizes  the  propellant  by  electron 
bombardment.  It  employs  permanent  bar  magnets  which  run  along  the  discharge 
chamber,  as  shown  in  Figure  1.2.  The  electrons  become  trapped  in  the  magnetic  field 
between  the  cathode  and  anode.  This  thruster  has  separate  regions  for  ionization  and 
acceleration.  The  thruster  configuration  is  shown  in  Figure  1.2.  This  thruster  requires  a 
cathode  inside  of  the  discharge  chamber  for  propellant  ionization  and  another  cathode 
outside  of  the  thruster  for  ion  beam  neutralization.  The  propellant  is  ionized  in  the 
discharge  chamber  and  then  extracted  by  the  screen  grid.  Ion  beam  acceleration  occurs 
between  the  acceleration  and  screen  grids.  Ion  thrusters  operate  at  lower  current  densities 
than  Hall  thrusters  because  space-charge  effects  between  the  grids  limit  the  ion  current 
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density;  however,  an  ion  thruster  can  accelerate  ions  through  much  higher  voltages  than 
the  single  stage  Hall  thruster  discussed  in  this  thesis.  A  typical  30  cm  ion  engine  is 
capable  of  generating  92  mN  of  thrust  at  62%  efficiency  with  a  specific  impulse  of  3100 
s.  Specific  impulses  above  10,000  s  have  been  achieved  by  ion  thrusters. 


Figure  1.2  Cross-section  of  an  ion  thruster  (courtesy  of  Matt  Domonkos). 

1.1.4  Oxygen  Thrusters 

Hall  thrusters  that  use  oxygen  as  the  propellant  are  being  developed  for  missions 
in  which  in-situ  resource  utilization  is  a  necessity.  Lunar  surface  material  is  50  % 
oxygen  by  weight,6  and  95  %  of  the  Martian  atmosphere  is  COj.  If  oxygen  is  extracted 
from  these  environments,  it  can  be  used  as  propellant  for  lunar  cargo  vehicles  or  Mars 
sample  return  missions  when  necessary.  Acquiring  propellant  at  the  mission  target  for 
the  return  trip  could  significantly  reduce  the  launch  costs  of  a  mission.  In  the  future,  if 
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and  when  these  missions  occur  regularly,  using  this  propulsion  system  could  prove  to  be 
rewarding. 

Before  oxygen  Hall  thrusters  can  be  utilized,  their  efficiency  must  be  improved7 
and  the  hollow  cathodes  typically  used  with  Hall  thrusters  must  be  replaced  with  cold 
cathodes.  Hollow  cathodes  employ  thermionic  electron  emitters,  which  are  very  sensitive 
to  the  operating  environment.  In  an  oxygen  environment,  'these  electron  emitting  cathode 
inserts  quickly  oxidize  and  their  work  function  increases.  To  overcome  the  increase  in 
work  function  and  maintain  a  constant  current,  higher  power  levels  are  required  for  the 
heaters  to  attain  higher  operating  temperatures  for  the  insert.  At  high  enough 
temperatures  surface  oxidation  rates  significantly  decrease;  however,  the  lifetime  of  the 
emitters  will  also  significantly  decrease.  If  Hall  thrusters  are  to  successfully  operate  on 
oxygen  propellant,  cold  cathodes  must  replace  hollow  cathodes  in  these  systems. 

1.1.5  Meso-  and  Microscale  Thrusters 

MesoO  and  microscale  thrusters  are  being  developed  for  spacecraft  in  the  1-200 
kg  wet  mass  range  and  for  fine  attitude  control  of  larger  spacecraft.  The  interest  in  such 
small  spacecraft  is  motivated  by  the  reduction  of  launch  msts  and  mission  risk.  It  has 
been  proposed  that  particular  missions,  like  exploring  the  rings  of  Saturn,  be 
accomplished  by  fleets  of  microcraft  possibly  in  communication  with  a  larger  mother 
spacecraft.8  Such  spacecraft  configurations  provide  a  high  degree  of  mission  flexibility.9 
Larger  spacecraft  could  use  microscale  electric  propulsion  systems  for  continuous 
disturbance  torque  compensation  and  drag  make-up.  Such  missions  include  laser 
interferometer  space  antenna  (LISA)10  for  the  detection  of  gravitational  waves,  Gravity 
Recovery  and  Climate  Experiment  (GRACE)  follow-on  mission,  and  Terrestrial  Planet 
Finder.  Incredibly  precise  position  maintenance  is  required  for  both  of  these  missions. 
Microscale  propulsion  systems  are  a  necessity  for  the  success  of  this  type  of  mission. 
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While  microscale  propulsion  system  technology  is  still  very  immature,  a  mesoscale  50  W 
Hall  thruster  has  been  developed.11  The  primary  technological  hurtle  for  this  propulsion 
system  is  now  finding  a  compatible  cathode.  The  cathodes  commonly  used  with  larger 
scale  Hall  thrusters  are  not  presently  compatible  with  this  thruster.  A  1/4”  cathode 
operates  at  25  W,  and  a  1/8”  hollow  cathode  operates  at  7-10  W.  These  cathodes 
consume  power  that  can  be  as  high  as  50  %  of  the  total  power  of  a  mesoscale  system. 
Some  mesoscale  thrusters  at  50  W  will  require  ~1.5  seem  of  propellant,  and  an  1/8” 
cathode  will  require  0. 5-1.0  seem.  Therefore,  hollow  cathodes  can  also  consume  50%  of 
the  propulsion  system  propellant.  This  problem  will  be  even  more  severe  with 
microscale  propulsion  systems  requiring  an  electron  source. 

1.1.6  Conventional  Electric  Propulsion  Cathode  Technology 

Cathodes  conventionally  used  with  Hall  and  ion  thrusters  are  too  inefficient  to  be 
used  with  meso-  and  microscale  propulsion  systems  and  are  easily  contaminated  in 
oxygen  environments.  Hollow  and  filament  cathodes  are  most  commonly  used  with  Hall 
and  ion  thrusters.  Hollow  cathodes  employ  lanthanum  hexaboride,  LaB6,  or  barium 
oxide,  Ba02,  impregnated  tungsten  inserts  which  require  heaters  to  increase  their 
temperature  up  to  -1000  °C  for  sufficient  electron  emission  before  thruster  start-up. 

They  also  require  propellant  to  create  plasma  which  ‘contacts’  the  electron  beam  to  the 
ion  beam.  Ions  created  in  this  plasma  supply  thermal  energy  to  the  insert  by  collisional 
energy  transfer,  and  neutralize  the  negative  space-charge  of  the  electron  beam.  This  local 
plasma  increases  the  electron  transmission  from  the  cathode  to  the  thruster  discharge 
chamber  and  ion  beam  at  voltages  -20  V.  State-of-the-art  1/8”  hollow  cathodes  under 
development  at  NASA  Glenn  Research  Center  (GRC)  and  the  Plasmadynamics  and 
Electric  Propulsion  Laboratory  (PEPL)  at  the  University  of  Michigan  consume  -7  W  and 
0.1  mg/s  of  xenon  for  -0.1  A.12  Filament  cathodes  can  also  be  used  with  Hall  and  ion 
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thrusters,  and  field  emission  electric  propulsion  systems  (FEEP);  however,  they  also 
consume  significant  power.  An  oxide  thermionic  cathode  is  being  developed  for  FEEP 
systems  at  ESTEC  with  an  anticipated  performance  of  5  mA/W.13 

1.1.7  Advanced  Cathode  Requirements 

The  development  of  small,  meso-,  and  microscale  propulsion  systems,  and 
oxygen  thrusters  for  spacecraft  will  be  enabled  by  the  development  of  cold  cathodes, 
which  do  not  require  propellant.  Such  cathodes  should  be  easily  scalable  in  size  to  less 
than  1  cm2  and  power  to  less  than  1  W/A  to  be  used  by  small,  meso-,  and  microscale 
thrusters.  These  cathodes  must  survive  for  1000s  of  hours  eventually;  however,  initial 
performance  objectives  are  only  100  raA/cm2,  for  100  hr.  in  simulated  thruster 
environments. 


1.2  Space  Tethers 

Space  tethers  are  conducting  or  non-conducting  cables  which  are  used  to  couple 
spacecraft  together  for  energy  transfer  or  for  propulsive  maneuvers.  They  typically 
consist  of  a  center  conducting  wire  surrounded  by  several  layers  of  thin  strands  of  high 
strength  fibers.  It  has  been  proposed  to  use  tethers  for  orbit  transfer  and  for  orbit 
insertion  from  the  Earth’s  or  Moon’s  surface  using  rotating  tethers.14,15  The  action  of 
tethers  depends  on  orbital  mechanics.  A  tether  dropped  down  from  a  spacecraft  is  pulled 
down  to  the  Earth  while  a  tether  extended  up  from  the  spacecraft  in  orbit  is  pulled  up  by 
centrifugal  forces.  Electrodynamic  tethers  use  interactions  between  Earth’s  magnetic  field 
and  current  in  the  tether  for  power  or  propulsion.16  SEDS-2  mission  demonstrated  that  a 
20  km  long  tether  could  successfully  de-orbited  at  payload  from  LEO.  The  TSS-1R 
mission  demonstrated  that  tethers  can  be  used  to  boost  satellites  into  higher  orbits.  The 
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terminator  tether  is  one  type  of  electrodynamic  tether  that  can  be  used  to  de-orbit 
satellites. 

1.2.1  Electrodynamic  Tethers 

It  iias  been  proposed  to  use  the  electrodynamic  (ED)  space  tethers  to  de-orbit 
satellites  at  the  end-of-life.17  One  concept  for  using  ED  tethers  to  de-orbit  spacecraft  is 
the  Terminator  Tether™.17  These  ED  tethers  use  passive  electrodynamic  tether  drag  to  de¬ 
orbit  a  satellite  from  LEO.  As  the  vertical  tether  is  pulled  through  the  mostly  horizontal 
Earth  magnetic  field  by  the  spacecraft,  a  voltage  is  induced  in  the  tether  wire.  Wire  ends 
make  contact  with  the  space  plasma  with  one  end  emitting  electrons.  The  current  is 
conducted  through  the  space  plasma  to  the  other  end  of  the  tether  where  electron 
collection  occurs.  The  Earth’s  magnetic  field  induces  a  JxB  force  on  the  wire  which  . 
opposes  the  host  spacecraft  motion.  Electrodynamic  drag  is  many  orders  of  magnitude 
greater  than  aerodynamic  drag  for  orbits  above  500  km.17  A  700  kg  satellite  at  1475  km 
could  be  de-orbited  in  6  months  by  a  10-km-long  tether.17  The  mass  of  the  tether  system 
is  less  than  a  third  of  the  propellant  mass  required  by  a  resistojet  if  the  electron  emitter  on 
the  tethr.  system  does  not  require  a  heater  or  consumable  gas  supply.17 

1.2.2  Cathode  Requirements 

Electron  emitters  used  with  an  ED  tether  de-orbiting  system  must  operate  cold 
and  without  propellant,  and  be  capable  of  as  much  as  2  A  to  make  it  competitive  with 
other  satellite  de-orbiting  options.  It  is  also  important  that  cathode  materials  are  resistant 
to  oxidation  in  Low  Earth  Orbit  (LEO).  The  cathode  current  requirement  could  be 
satisfied  by  using  several  cathodes.  The  tether  system  could  use  a  single  hollow  cathode 
as  the  electron  source,  however,  the  mass  penalty  of  this  cathode  system  may  negate  the 
advantage  of  the  ED  over  the  resistojet  to  de-orbit  satellites.  The  use  of  cold  cathodes 
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that  do  not  require  a  consumable  gas  could  return  the  competitive  edge  to  the  ED  tether 
concept. 


1.3  Field  Emission(FE)  Cathodes 

A  FE  cathode  is  a  promising  alternative  to  thermionic  hollow  cathodes  because 
they  are  small,  cold,  efficient,  chemically  inert,  and  do  not  require  an  ionizable  gas. 
Figure  1.3  shows  a  Hall  thruster  with  FE  cathodes  in  two  possible  positions:  in  the  center 
of  the  thruster  (internal)  and  on  the  periphery  (external).  FE  cathodes  employ  strong 
electric  fields  instead  of  high  temperatures  to  emit  electrons.  FE  cathodes  require  much 
simpler  and  smaller  systems  to  operate  them  than  conventional  cathodes.  They  can  be 
easily  scaled  in  size  and  power  level  to  be  compatible  with  microscale  propulsion 
systems,  and  have  demonstrated  significantly  higher  efficiency  than  thermionic 
cathodes.18  Two  types  of  FE  cathodes  were  investigated  in  this  thesis:  thin  film  carbon 
FE  cathodes,  and  microcone  field  emission  array  (FEA)  cathodes. 


Figure  1.3  A  Hall  thruster  shown  with  hollow  and  field  emission  cathodes. 
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1.3.1  Field  Electron  Emission 


Field  electron  emission  occurs  when  applied  electric  fields,  F,  at  the  cathode 
surface  are  high  enough  (0.3-0.6  V/A)  to  deform  the  potential  barrier,  <|>w,  between  the 

cathode  surface  and  vacuum  to  allow  electrons  to  tunnel  out  of  the  cathode  surface  and 
through  the  deformed  potential  barrier.19  Potential  barrpr  deformation  by  an  electric  field 
is  shown  in  Figure  1.4.  The  applied  electric  field  is  F,  the  work  function  is  <j)w,  the  Fermi 

energy  is  Ep,  and  e  is  the  electron  charge.  The  Fermi  energy  is  the  maximum  energy 

level  occupied  by  electrons  if  all  of  the  lower  energy  levels  are  filled  in  the  quantum 
structure.  The  height  of  the  potential  barrier  between  the  surface  and  vacuum  is  <)),  and  the 

width  of  the  barrier  is  <j )JF.  The  electrons  emitted  by  field  emission  originate  below  the 

Fermi  energy  level,  therefore  these  low  energy  electrons  are  considered  cold.  Field 
emission  differs  significantly  from  thermionic  or  photoemission  where  electrons  in  the 
solid  must  be  given  enough  energy  to  go  over  the  potential  barrier.  These  electrons  have 
energy  considerably  greater  than  the  Fermi  energy,  and  are  considered  hot. 


Figure  1.4  Potential  energy  diagram  for  electrons  at  a  metal  surface,  without  (left) 
and  with  (right)  an  applied  electric  field,  F. 


Spindt-type  FEA  cathodes  with  microtip  structures  require  electric  field 
enhancement  by  the  sharp  tip  structure.  Electric  fields  at  the  emitting  tips  in  excess  of 
4xl07  V/cm  are  required  for  field  emission.  The  field  emission  tip  configuration,  which 
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is  typically  employed,  is  shown  in  Figure  1.5;  dimensions  and  materials  are  processor 
dependent.  The  Fowler-Nordheim  (F-N)  equation  for  field  emission  current  density  at  an 
emission  site  is  Eqn.  1.1.  The  F-N  coefficients  are  a*,  and  bfE  are  given  in  Eqn.  1.2  and 
Eqn.  1.3  where  t(y)=1.05657,  vo=0. 93685,  m=511000  eV/c2 ,  c=2998.79  A/fs, 

Q=a  h c/4=3.6eV-A,  a=l/137.04,  and  ^=0.655eV-fs.  A  more  detailed  derivation  of 


these  relationships  is  given  elsewhere.20,21 


Figure  1.5  Field  emission  tip  and  gate  electrode  configuration  (courtesy  of  Dr. 
Kevin  Jensen). 
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1.3.2  State-of-the-Art  FE  Cathode  Technology 

The  applications  discussed  in  this  thesis  require  state-of-the-art  FE  cathode 
technology  for  materials,  cathode  structure,  and  current  limiting  architectures.  This 
technology  is  evolving  as  microfabrication  techniques  mature  to  deposit  more  robust 
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materials,  and  smaller  and  more  complicated  structures.  Microfabrication  techniques  are 
applied  to  deposit  and  etch  insulating  and  conducting  films  and  arrays  of  microcones  on 
silicon  wafers.  Electron  beam  and  interferometric  lithography  techniques  are  used  for 
patterning  the  cathodes.  Tip  and  gate  aperture  radii  are  on  the  order  of  1  and  100  nm, 
respectively.22  Packing  densities  greater  than  108  tips/cm2  have  been  achieved. 

Materials  most  commonly  used  are  silicon  and  molybdenum,  however,  diamond-like 
carbon  films  and  coatings  on  the  microtips  have  been  used  to  improve  robustness  of  the 

tips,  work  functions,  chemical  inertness,  and  emission  stability.  ’  '  ’  Molybdenum  FEA 
cathodes  fabricated  at  the  Massachusetts  Institute  of  Technology  Lincoln  Laboratory 
(MIT/LL)  have  attained  1  pA  at  25  V  with  an  array  of  900  tips.  The  same  cathodes  with 
cesiated  tips  demonstrated  1  pA  at  10  V.24  Spindt  at  SRI  International  has  reported  dc 
field  emission  up  to  100  pA/tip  from  a  100-tip  array.29  Emission  current  densities  greater 
than  2000  A/cm2  have  been  achieved.30,31  The  current  collected  by  the  gate  electrodes  is 
typically  1/1000  of  the  emitted  current.  For  optimum  cathode  performance  typical  tip 
operating  pressure  is  ~10'9  Torr  after  careful  conditioning  in  hydrogen  and  neon 
environments.32,33'34  Microtip  FEA  cathodes  at  SRI  International  have  demonstrated 
performance  of  100  mA  with  less  than  1  mW  consumed  by  the  gate  electrode.35 

Alternatives  to  the  micro-fabricated  PEA  cathodes  are  diamond,  diamond-like 
carbon,  or  carbon  thin  film  FE  cathodes.  The  unique  advantage  to  using  doped  carbon  is 
its  low  electron  affinity.36  These  cathodes  can  emit  electrons  at  macroscopic  vacuum 
electric  field  strengths  that  are  two  orders  of  magnitude  lower  than  those  fields  required 
by  silicon  and  molybdenum  Spindt-type  cathodes  and  performance  is  not  as  sensitive  to 
sputtering  by  ion  bombardment.36,37,38  Diamond  cathodes  are  very  resistant  to  oxidation 
while  some  carbon  films  have  demonstrated  performance  degradation  in  oxygen 
environments. 

Carbon  FE  cathodes  are  more  robust  than  the  Spindt-type  FEA  cathodes,  however 
the  Spindt-type  cathodes  have  typically  demonstrated  higher  current  densities  and 
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efficiency.  Various  types  of  carbon  cathodes  have  been  investigated.  It  is  common  for 
carbon  film  cathodes  to  operate  with  only  20%  of  the  emitted  current  escaping  through 
the  gate  electrode  and  being  collected  by  the  anode.  A  40  mA  current  was  measured  at 
Field  Emission  Picture  Element  Technology  (FEPET)  Inc.  from  a  25  mm  carbon  film 
cathode  with  an  extraction  field  of  6.7xl04  V/cm.37  An  emission  current  density  of  100 
mA/cm4  was  achieved  with  an  extraction  field  of  less  than  1.2x10s  V/cm.  Blyablin  et  al. 
demonstrated  500  mA/cm2  from  diamond  film  cathodes.  Cathodes  with  diamond  gnt  in 
microfabricated  gated  structures  have  achieved  100  mA/cm2  at  20  V.36 

1.4  Challenges 

State-of-the-art  FE  cathodes  have  demonstrated  excellent  performance  in  UHV 
environments  with  large  anode  voltages  and  promising  results  in  some  higher  pressure 
environments.  Thruster  and  tether  environments  could  significantly  limit  the 
performance  of  these  cathodes.  FE  cathodes  have  demonstrated  lifetimes  of  1000’ s  of 
hours  in  UHV  environments.  In  higher-pressure  environments  cathode  lifetime  can  be 
shortened  by  arcing  between  the  gate  electrodes  and  tips,  and  by  the  cathode  getting 
sputter  .d  by  ion  bombardment.  In  low-pressure  environments  without  a  close-spaced 
anode,  electron  emission  can  be  space-charge  limited. 

1.4.1  Cathode  Lifetime 

In  UHV  environments,  the  primary  cathode  failure  mode  has  been  arcing  between 
the  gate  and  tip,  shorting  out  the  cathode  and  ending  its  life.  In  higher  pressure 
environments,  this  threat  is  even  greater.  Recent  advances  in  microfabrication  should 
eliminate  this  cathode  failure  mode.  Catastrophic  arcing  events  can  be  prevented  by 
using  current  limiting  structures  such  as  resistive  layers,40  field  effect  transistors  (FET),41 
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and  vertical  current  limiting  (VECTL)42  architectures.  Robust  materials  can  be  used  to 
improve  cathode  sputter  resistance  and  resistance  to  oxidation. 43,44,45 

Cathode  performance  will  also  be  affected  by  ion  bombardment.  There  are  two 
populations  of  ions  which  will  bombard  the  cathode  tips  in  electric  propulsion  systems 
such  as  Hall  thrusters.  One  population  is  created  between  the  tips  and  gate  electrode,  and 
the  other  population  originates  ftear  the  thruster  where  charge  exchange  (CEX)  collisions 
occur  between  beam  ions  and  neutrals.  The  trajectories  of  beam  ions  will  not  intersect 
the  cathode,  however,  they  will  create  CEX  ions  that  can  be  accelerated  into  the  cathode 
region.  Measurements  of  this  CEX  current  are  discussed  in  Chapter  2.  The  energy  of 
these  ions  depends  on  the  potential  between  the  cathode  and  local  plasma  as  shown  in 
Figure  1.6. 

The  second  population  of  ions  is  generated  near  the  electron  emitting  surface 
between  the  tip  and  gate  electrode  where  emitted  electrons  ionize  neutrals  by  ion 
bombardment.  If  ions  are  accelerated  into  the  microtips  through  large  enough  voltage 
differences,  these  ions  will  damage  the  cathode  surface.  The  flux  of  these  ions  is  directly 
dependent  on  the  local  pressure.  FEA  cathodes  have  demonstrated  stable  emission  in 
elevated  pressure  environment*!  vnen  start-up  occurs  at  ~10'9  Torr  and  the  pressure  is 
slowly  increased.18  Spindt  et  al.18  showed  that  these  cathodes  could  operate  in  ~10‘5  Torr 
of  02,  Ar,  He,  H2,  and  N2  for  hundreds  of  hours.  During  these  exposure  tests,  the 
emission  current  dropped  by  about  an  order  of  magnitude,  if  at  all.  When  the  chamber 
returned  to  the  original  UHV  environment,  the  cathodes  demonstrated  full  recoveries. 
Similar  experiments  at  the  Microelectronics  Center  of  North  Carolina  (MCNC)  yielded 
similar  results  for  silicon  cathodes.46,47  The  results  of  these  experiments  showed  that 
increases  in  pressure  increased  cathode  work  function  temporarily.  Exposure 
experiments  performed  with  xenon  at  SRI  International  showed  that  cathodes  were 
severely  damaged  in  a  xenon  environment.  Tests  at  Linfield  Research  Institute  (LRI)  of 
single  ZrC  and  HfC  tips  with  kilovolt  gap  voltages  showed  that  emission  continued  for 
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several  minutes  at  KT4  Torr  of  Ar  and  02.  Although  a  glow  discharge  operating  mode 
was  possibly  attained,  the  tips  were  not  destroyed  during  operation  for  a  few  minutes 
with  ions  bombarding  the  tips  with  several  keV  energies.  The  results  of  these 
experiments  in  elevated  pressure  environments  are  very  promising,  especially  because 
some  of  the  cathodes  seemed  to  be  undamaged  during  operation. 

Carbon  cathodes  are  more  robust  than  the  Spiiidt-type  cathodes,  however  the 
Spindt-type  cathodes  have  typically  demonstrated  higher  current  densities  and  efficiency. 
These  cathodes  can  operate  at  400  V  at  pressures  greater  than  750  mTorr  of  nitrogen 
without  being  damaged.  While  performance  degradation  has  been  observed  for  the 
diamond  and  carbon  film  cathodes  while  operating  in  oxygen-rich  environments38,^  is 
anticipated  that  the  cathode  performance  will  be  much  more  stable  in  the  inert  gas 
environment  of  a  xenon  plasma  thruster. 

The  effect  of  higher  pressure  environments  on  the  lifetime  of  FEA  cathodes  has 
only  been  investigated  to  determine  the  removal  rate  of  material  from  FE  microtips.  This 
model  was  developed  by  Brodie.48  His  model  is  used  to  calculate  the  number  of  ions 
bombarding  the  tips  and  rate  of  material  removed  from  the  emitting  area  of  the  tips,  and 
depends  on  the  sputter  yield  of  the  ion-target  combination,  and  ionization  energy, 
pressure,  and  temperature  of  the  ambient  gas.  He  assumed  an  average  ion  energy  and 
corresponding  ionization  cross-section,  Q,  and  sputter  yield,  Y,  and  employed  operating 

voltages  around  100  V.  He  showed  that  a  Mo  tip  with  a  500  A  radius  of  curvature, 
operating  at  100  V,  at  10'10  Torr  of  hydrogen,  emitting  1  pA,  will  lose  10  layers  of 

material  from  ion  bombardment  in  2xl07  hours.  With  a  linear  relationship  with  pressure, 
the  cathode  will  lose  10  layers  of  material  in  200  hours  at  10‘5  Torr.  The  material  lost 
rate  will  increase  significantly  with  increasing  ion  mass.  This  model  was  developed  for 
single  tip  configurations;  however,  a  similar  model  could  also  be  used  for  microtip 
structures  operating  at  much  lower  voltages.  At  lower  operating  voltages,  average  ion 
energies  should  not  be  assumed  because  of  the  exponential  nature  of  the  relationship 
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between  Q  and  electron  energy  and  S  and  ion  energy.  Many  models  have  been  developed 
to  predict  the  performance  of  FEA  cathodes  as  a  function  of  tip  radius.  Performance 
models  have  never  been  coupled  with  Brodie’s  model  to  determine  the  effect  of  cathode 
sputtering  in  higher  pressure  environments  on  cathode  performance. 

1.4.2  Space-Charge  Limited  Emission 

The  second  concern  with  integrating  FE  cathode  with  EP  and  tether  technologies 
is  space-charge  limited  electron  emission  when  emitting  into  plasma.  FE  cathodes 
typically  employ  a  close-spaced  triode  configuration  with  a  physical  anode  only 
millimeters  away  from  the  gate  electrode  to  collect  the  emitted  current.  Hall  and  ion 
thrusters  employ  anodes  and  cathodes  separated  by  a  plasma  with  pressures  up  to  10'4 
Torr  in  the  ionization  regions.  In  this  environment,  a  virtual  anode  will  collect  the 
emitted  electron  current  at  the  cathode  sheath  boundary.  The  thickness  of  this  sheath  and 
potential  drop  through  it  depend  on  the  environment,  cathode  potentials,  and  properties  of 
the  emitted  current.  Hall  thrusters  typically  operate  in  a  configuration  where  the  cathode 
is  floating.  The  cathode  sheath  voltage  is  approximately  20  V  in  a  1.35  kW  Hall  thruster 
system  with  a  300  V  anode-cathode  voltage  difference.  Hollow  cathodes  us**  propellant 
in  addition  to  the  primary  supply  to  the  anode  to  generate  a  plasma  at  the  emitter  surface 
to  reduce  space-charge  effects  and  the  cathode  sheath  voltage.  A  reduction  in  sheath 
voltage  improves  the  efficiency  of  the  thruster  since  less  of  the  discharge  voltage  is  tied 
up  in  the  cathode  fall  and  a  higher  voltage  is  available  to  accelerate  the  ions  to  produce 
thrust.  It  is  more  crucial  that  the  cathode  potential  is  low  with  respect  to  plasma  potential 
when  FE  cathodes  are  used  because  of  the  lifetime  limitations  induced  by  ions  sputtering 
the  microtip  structures  and  gate  electrode.  It  is  important  that  the  gate  electrode  potential 
is  lower  than  the  plasma  potential  to  retard  the  collection  of  plasma  electrons  while  the 
cathode  beam  electrons  escape  to  the  plasma.  It  is  anticipated  that  the  field  emitter  gate 
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electrode  will  be  set  at  approximately  20  V  below  plasma  potential,  <j)c.  A  potential 

diagram  is  shown  in  Figure  1.6.  The  thruster  discharge  voltage  between  the  anode  and 
cathode,  Vd,  and  the  gate  voltage  with  respect  to  the  cathode  FE  surface,  <|>g,  are  shown. 

In  a  plasma  environment  with  virtual  anodes,  space-charge  effects  will  limit  the  electron 
current  that  escapes  through  the  sheath  to  the  plasma.  This  concern  must  be  properly 
addressed  and  current  limits  must  be  determined  using  sheath  modeling  before  FE 
cathodes  can  be  designed  and  used  in  plasma  environments. 


Figure  1.6  Potential  diagram  of  the  cathode,  thruster,  and  plasma  configuration. 
Anode,  plasma  ion  beam,  sheath  boundary,  gate  electrode,  and  cathode  potentials 
are  shown  relative  to  ground. 


Space-charge  limited  current  problems  have  been  studied  extensively.  The 
Childs-Langmuir  relationship  describes  space-charge  current  limitations  in  a  vacuum 
diode.  In  the  derivation,  the  current  was  limited  by  a  zero  electric  field  at  the  cathode 
surface  because  initial  electron  velocities  were  negligible.49,50  Fay  et  al.  determined 
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emission  current  limitations  in  a  parallel  plate  geometry  when  electrons  are  emitted  with 
considerable  velocities  in  a  vacuum.51  In  that  case,  positive  electric  fields  at  the  cathode 
surface  are  stable  solutions  to  the  Poisson  equation,  permitting  larger  currents  to  escape 
from  the  cathode.  Virtual  cathodes  can  eventually  form  between  the  cathode  and  anode 
to  limit  the  emitted  current.  Langmuir  considered  space-charge  limitations  for  parallel 
plates  separated  by  a  plasma,  where  ions  contribute  to  negative  space-charge 
neutralization  to  further  increase  space-charge  limited  current.  He  claimed  that  electron 
and  ion  currents  were  limited  in  proportion  to  ion  and  electron  masses  when  ion  and 
electron  initial  velocities  are  negligible,  but  showed  that  these  currents  can  increase  with 
non-zero  velocities  at  the  sheath  boundary.52  Bohm  modeled  the  sheath  near  the  cathode 
and  determined  that  for  a  cathode  sheath  to  be  stable,  the  ions  must  enter  it  with  non-zero 
initial  velocity.53  To  ensure  that  potential  minima  occur  only  at  sheath  boundaries  for  a 
non-emitting  cathode,  it  was  determined  that  vi>(Te/mi)Al/2  at  the  sheath  boundary. 
Crawford  and  Cannara54  considered  the  case  for  space-charge  limited  emission  of  a  hot 
cathode  into  plasma.  They  determined  a  valid  velocity  range  for  ions  entering  the  sheath, 
the  velocity  being  limited  by  charge-exchange  and  momentum  exchange  collisions  in 
pre-sKaths.  Prewett  and  Allen55  studied  the  same  configuration,  and  included  a  criterion 
derived  by  Andrews  and  Allen56  to  construct  an  analytical  expression  for  the  initial 
energy  of  ions  entering  a  sheath.  Using  that  expression  in  the  sheath  model,  they 
achieved  excellent  correlation  between  theory  and  experiment.  Goodfellow  augmented 
this  model  by  considering  the  initial  thermal  velocities  of  the  electrons  emitted  at  the 
cathode.  He  used  this  sheath  model  combined  with  thermal  models  to  predict  cathode 
temperatures  in  a  magnetoplasmadynamic  thruster  system.  He  also  achieved  good 
correlation  between  experiment  and  theory.57 
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1.5  Role  of  this  Thesis 


It  is  anticipated  that  the  successful  integration  of  FE  cathodes  with  EP  systems 
will  require  the  modification  of  currently  existing  Hall  thrusters,  more  advanced  FE 
cathode  technology,  and  an  understanding  of  how  these  cathodes  will  operate  in  plasma 
environments.  Thrusters  will  require  different  start-up  sequences  than  commonly  used, 
to  protect  the  cathode.  Lower  work  function  materials,  smaller  cathode  dimensions,  and 
current  limiting  structures  will  significantly  improve  the  performance  and  compatibility 
of  these  cathodes  with  electric  propulsion  systems.  Testing  of  FE  cathodes  in  Xe 
environments  similar  to  Hall  thruster  environments  is  necessary  to  study  the  effect  of 
higher  pressure  environments  and  operating  voltages  on  cathode  lifetime.  This  process 
should  also  be  modeled  to  expedite  the  development  and  integration  process,  reduce 
costs,  and  to  determine  how  to  optimize  cathode  design  for  different  applications.  Also, 
space-charge  limited  emission  should  be  modeled  for  the  FE  cathode  and  virtual  anode 
configuration  to  develop  and  understanding  of  the  performance  and  cathode  current 
limitations  due  to  space-charge  effects  in  plasma  environments.  Some  of  the  modeling 
results  should  be  experimentally  validated.  This  thesis  tackles  a  few  of  these 
technological  hurtles  including  performance  assessment  of  available  FE  cathodes  in 
simulated  thruster  environments,  the  development  of  a  model  of  the  performance 
degradation  due  to  cathode  sputtering  by  ion  bombardment,  and  a  model  to  predict  space- 
charge  limited  emission  from  a  FE  cathode  in  a  plasma  environment. 

The  primary  objectives  of  this  work  were  to  determine  if  field  emission  cathodes 
can  operate  in  thruster  environments  and  what  limitations  should  be  expected.  Cathode 
lifetime  and  space-charge  limited  electron  emission  are  focal  points  in  this  investigation. 
It  is  assumed  that  FE  cathodes  with  current  limiting  configurations  will  soon  be  available 
to  eliminate  arcing  between  tips  and  gate  electrodes  as  a  cathode  failure  mode. 
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The  first  stage  of  this  research  focused  on  characterizing  the  cathode  environment 
of  a  Hall  thruster.  A  70  mm  Hall  thruster,  PEPL-70,  was  built  and  used  in  the  cathode 
environment  characterization.  The  performance  of  the  thruster  was  evaluated  to  ensure 
nominal  operating  conditions  were  achieved.  Ion  current  densities  were  measured  in  the 
cathode  region  of  the  thruster  using  a  Faraday  probe.  4  neutral  particle  flux  (NPF)  probe 
was  used  to  measure  the  local  pressure.  These  measurements  were  used  to  design 
experiments  and  model  the  effect  of  the  thruster  environment  on  FE  cathode 
performance.  The  PEPL-70  thruster  configuration  represents  only  one  of  the  possible 
applications  of  FE  cathodes  to  electric  propulsion..  Initial  FE  cathode  and  thruster  testing 
will  be  conducted  with  a  thruster  of  this  size  and  power  level,  therefore  this  thruster 
environment  was  simulated  for  initial  FE  cathode  testing. 

Experiments  were  conducted  to  evaluate  the  performance  of  FE  cathodes  in 
environments  representative  of  Hall  thruster  environments.  The  performance  of  carbide 
cathodes  was  evaluated  in  argon  and  oxygen  environments.  The  performance  response  of 
molybdenum,  silicon,  and  carbon  film  cathodes  to  increases  in  xenon  pressure  was 
investigated  at  different  operating  voltages  and  pressures.  The  objectives  of  these 
experiments  were  to  determine  if  xenon  exposure  duri-g  operation  would  change  the 
cathode  work  function  and  if  the  cathodes  could  be  operated  at  low  enough  voltages  to 
avoid  damaging  them  during  ion  bombardment. 

A  model  was  developed  to  predict  performance  degradation  of  microtip  cathodes 
in  xenon  environments.  This  model  uses  a  cathode  performance  model  developed  by 
Jensen  to  predict  effective  tip  radius  and  spread  in  radii  of  an  array  of  cathode  tips  from 
an  I-V  trace  to  predict  the  cathode  performance.  The  performance  degradation  model 
also  uses  a  model  developed  by  Brodie  to  predict  the  rate  of  removal  of  material  from  the 
microtip  emitting  area.  These  models  are  combined  with  models  which  predict  changes 
in  effective  tip  radius  and  spread  in  tip  radii  in  time  to  estimate  the  change  in  current 
emitted  by  a  field  emission  array  cathode  during  xenon  ion  bombardment.  The  model 
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flow  diagram  is  shown  in  Figure  1.7.  Experimental  data  and  this  performance 
degradation  model  were  used  together  to  estimate  the  energy  thresholds  for  sputtering  of 
Mo  and  Si  by  xenon  ions.  With  this  model  the  performance  degradation  can  be  predicted 
for  any  Mo  or  Si  FEA  cathode  operating  in  a  Xe  environment.  Good  agreement  was 
achieved  between  modeling  and  experimental  results.  This  model  can  also  be  used  as  a 
tool  to  optimize  the  cathode  configuration  for  various  applications. 


Figure  1.7  Flow  diagram  for  the  field  emission  array  cathode  performance 
degradation  model. 


A  model  was  developed  to  predict  space-charge  current  limitations  of  FE 
cathodes  emitting  into  thruster  and  tether  environments.  This  model  is  based  on  the 
planar  sheath  models  developed  by  Goodfellow  and  Prewett  and  Allen.  The  criterion 
derived  by  Andrews  and  Allen  was  used  to  estimate  the  ion  velocity  at  the  boundary 
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between  the  sheath  and  pre-sheath.  The  model  developed  in  this  research  considers 
electrons  emitted  with  energies  greater  than  20  eV.  This  model  also  considers  cylindrical 
and  spherical  sheath  geometries  to  model  situations  in  which  electron  beam  expansion 
occurs  in  the  sheath  as  1/r  and  1/r2,  respectively.  Beam  expansion  must  be  considered  in 
the  cases  when  the  sheath  dimensions  are  larger  than  the  cathode  dimensions.  These 
models  were  used  to  investigate  the  effect  of  cathode  environment  and  dimensions  on  the 
space-charge  current  limit. 

Finally,  the  modeling  and  experimental  results  are  applied  to  design  cathodes  for 
thruster  and  tether  applications.  Preliminary  cathode  design  objectives  are  -100  mA 
from  an  area  smaller  than  1  cm2  for  100  hrs.  Cathode  lifetime  and  space-charge  current 
limitations  are  considered  in  the  designs.  Cathode  geometrical  and  operating 
configurations  are  suggested  for  the  different  applications. 
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CHAPTER  2 


CATHODE  ENVIRONMENT  CHARACTERIZATION 

The  cathode  environment  of  a  Hall  thruster  must  be  understood  to  determine  the 
compatibility  of  Hall  thrusters  and  FE  cathodes.  A  Hall  thruster  was  built  and  operated 
with  xenon  propellant  to  provide  a  plasma  environment  which  is  similar  to  the  plasma 
environment  of  a  commercial  system.  The  performance  of  the  thruster  was  evaluated 
with  I-V  and  thrust  measurements.  Faraday  and  Neutral  Particle  Flux  (NPF)  probes  were 
used  to  measure  ion  current  density  and  neutral  particle  pressure  in  the  region  of  a  Hall 
thruster  where  FE  cathodes  could  be  positioned.  Electron  number  densities  and 
temperatures  on  a  similar  Hall  thruster  system  were  measured  elsewhere.58  Parameters 
obtained  from  cathode  environment  characterization  were  used  to  design  FE  cathode 
experiments  in  xenon  environments  and  model  the  effect  of  the  thruster  environment  on 
cathode  performance  and  lifetime.  The  thruster  developed,  diagnostics  used  to 
characterize  the  cathode  environment,  and  the  measurements  obtained  in  the  cathode 
environment  characterization  are  discussed  in  this  chapter. 

2.1  Experimental  Apparatus 


2.1.1  Test  Facility 

These  experiments  were  conducted  at  the  Jet  Propulsion  Laboratory  (JPL).  The 
test  facility  consisted  of  a  2.3  x  4.6  m  vacuum  chamber  with  a  xenon  pumping  speed  of 
30,000 1/s  and  base  pressure  of  4xl0'6  Torr  maintained  by  both  diffusion  and  cryogenic 
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pumps.  The  pressure  gauge  was  calibrated  against  a  xenon  calibrated  ionization  gauge. 

A  JPL  hollow  cathode  was  used  to  supply  electrons  for  propellant  ionization  and  ion 
beam  neutralization.  Thruster  performance  was  evaluated  using  a  NASA  GRC  designed 
inverted  pendulum  thrust  stand.59  Calibration  data  are  shown  in  Figure  2.1.  The  thrust 
stand  is  an  inverted  pendulum  type  with  a  Linear  Voltage  Digital  Transducer  (LVDT) 
that  detects  movement  of  the  thrust  stand.  An  inclinometer  and  motor  are  used  to 
maintain  a  constant  thrust  stand  inclination  because  changes  in  inclination  can  be 
detected  by  the  LVDT  as  artificial  thrust.  A  motor  drives  the  deployment  of  a  chain  with 
calibration  weights  distributed  at  several  positions.  The  weights  are  deployed 
periodically  and  the  LVDT  displays  the  voltage  corresponding  to  the  thrust  stand 
displacement  during  calibrations.  During  the  calibrations,  the  weights  deployed  shift  the 
LVDT  from  its  zero  weight  (thrust)  position  associated  with  a  voltage  output,  V|,VDTo-  A 
change  in  voltage  from  the  zero  position  is  then  associated  with  each  calibration  weight 
to  estimate  thrust  during  thruster  operation.  This  thrust  stand  was  calibrated  for  much 
larger  thrust  levels  than  the  PEPL-70  generated;  therefore,  the  thrust  estimates  had  to  be 
extrapolated  from  the  calibration  curve.  Thrust  stand  response  has  been  shown  to  be 
linear  to  very  low  thrust  lev^L. 
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Figure  2.1  Thrust  stand  calibration  data. 


2.1.2  Hall  Thruster-PEPL-70 

Limited  Hall  thruster  availability  necessitated  the  fabrication  of  a  laboratory 
model  at  the  PEPL.  Fabrication  of  this  thruster  was  also  necessary  to  investigate 
performance  of  Hall  thrusters  with  oxygen  propellant.  Oxygen  can  be  very  destructive  to 
the  thruster.  Different  thruster  designs  and  material'’  are  optimal  for  operation  on  oxygen 
and  xenon.  The  boron  nitride  discharge  chamber  has  a  70  mm  average  diameter.  The 
titanium  anode  has  one  propellant  inlet  into  a  gas  distributor  with  48  outlets  in  the  first 
chamber  and  24  outlets  in  the  second  chamber.  Two  chambers  with  different  radial 
locations  of  the  propellant  outlets  improve  the  azimuthal  uniformity  of  xenon  in  the 
discharge  chamber.  The  magnetic  poles  in  this  thruster  are  slightly  thicker  than  the  pole 
pieces  commonly  used  in  xenon  thrusters  to  create  a  larger  ionization  region  because 
ionization  of  oxygen  requires  a  longer  path  length.  The  ionization  region  is  confined  by 
the  magnetic  field  generated  by  the  electromagnets.  The  thruster,  shown  in  Figure  1.3 
and  Figure  2.2,  employs  outer  and  inner  electromagnets  to  generate  a  radial  magnetic 
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field.  Hall  thrusters  commonly  used  4  outer  electromagnets;  however,  this  thruster 
employs  8  outer  electromagnets  to  improve  azimuthal  uniformity  in  the  magnetic  field. 
An  axial  electric  field  is  maintained  between  the  anode  and  cathode.  The  axial  electric 
field  and  the  radial  magnetic  field  create  a  closed  Hall  current  of  electrons.  Propellant 
atoms  injected  into  this  electron  cloud  through  the  anode-gas  distributor  are  ionized  by 
electron  bombardment.  The  ions  are  then  accelerated  axially  out  of  the  discharge 
chamber  by  the  electric  field  to  generate  thrust. 


Figure  2.2  The  PEPL-70  Hall  thruster  with  a  hollow  cathode,  neutral  particle  flux 
(NPF)  probe,  and  Faraday  Probe. 


2.1.3  Diagnostics 

Faraday  and  neutral  particle  flux  (NPF)  probes  were  used  to  characterize  the 
cathode  environment.  A  Faraday  button  probe  which  was  developed  at  JPL,  shown  in 
Figure  2.3,  was  used  to  measure  ion  current  densities  in  a  region  external  to  the  thruster 
where  a  FE  cathode  could  be  used.  A  schematic  for  the  Faraday  probe  is  shown  in  Figure 
2.3.  A  NPF60  probe  was  used  to  measure  neutral  particle  pressures  in  the  cathode  region. 
This  probe  measures  neutral  particle  pressure  by  using  a  series  of  retarding  grids  to 
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deflect  ions  entering  the  ionization  gauge.  A  thermocouple  provides  the  probe  wall 
temperature.  It  is  assumed  that  the  temperature  of  the  gas  is  in  equilibrium  with  the 
probe  walls.  The  NPF  probe  is  shown  in  Figure  2.3.  The  actual  positions  of  the  cathode 
and  probes  during  the  measurements  presented  in  this  chapter  are  different  than  the 
positions  shown  in  Figure  2.3.  The  cathode  was  at  a  3  o’clock  position,  the  NPF  at  10 
o’clock,  and  the  Faraday  probe  at  8  o’clock. 

Shield 


Figure  2.3  Faraday  probe  schematic  (courtesy  of  James  Haas). 


2.2  Experimental  Methods  and  Results 

2.2.1  Thruster  Performance 

During  several  hours  of  testing,  the  thruster  demonstrated  I-V  trends  comparable 
to  similar  Hall  thrusters  used  for  satellite  NSSK.  Figure  2.4  shows  I-V  data  at  different 
propellant  flow  rates.  A  properly  operating  Hall  thruster  will  demonstrate  a  steep  I-V 
curve  at  low  voltages,  followed  by  a  steep  drop,  and  then  a  plateau  in  current  with  further 
increases  in  discharge  voltage.  The  thruster  is  operating  optimally  in  the  plateau  region 
of  the  I-V  curve;  increases  in  discharge  voltage  increase  the  ion  exit  velocities  while  the 
ionization  fraction  of  the  propellant  does  not  change. 
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During  several  hours  of  testing,  the  thruster  demonstrated  short-term  performance 
comparable  to  similar  Hall  thrusters  used  for  satellite  north-south  station  keeping.  A  1.4 
kW  commercial  Hall  thruster  system  generates  4.5  A  at  300  V  and  53  seem  of  xenon  with 

83  mN  of  thrust  at  53  %  efficiency.3  The  PEPL-70  generated  4.6  A  and  65  mN  of  thrust 
at  36%  efficiency,  300  V,  and  43  seem  of  xenon.  Thrust  efficiency,  T]eff,  is  determined 

using  Eqn.  2.1,  where  m  is  the  mass  flow  rate  of  propellant  (mg/s),  TN  is  thrust  (N),  Id  is 
the  discharge  current  (A),  and  Vd  is  the  discharge  voltage  (V).  Performance  data  are 
shown  in  Figure  2.5.  The  error  bars  on  these  measurements  are  fairly  large  because  the 
thrust  stand  was  configured  to  measure  thrust  on  much  larger  propulsion  systems  with 
much  higher  thrust  levels,  and  the  zero  drift  of  the  LVDT  was  significant  during  the 
experiments.  During  the  experiments,  the  thrust  is  estimated  by  the  thrust  stand  position 
with  respect  to  the  initial  zero  thrust  position.  A  shift  in  the  thrust  stand  zero  position 
during  the  experiments  can  result  in  significant  over-  or  underestimates  in  thrust.  The 
most  accurate  method  of  determining  thrust  is  to  run  the  thruster  at  the  operating  point  of 
interest  for  at  least  0.5  hour  and  then  turn  it  off  to  determine  the  zero  point  of  the  thrust 
stand  and  the  thrust  associated  with  the  LVDT  voltage  reading  before  the  thruster  was 
turn'd  off.  Time  limitations  prevented  the  acquisition  of  such  accurate  thrust 
measurements  in  these  experiments.  Operating  points  were  only  explored  for  a  few 
minutes.  The  thrust  stand  was  calibrated  before  and  after  the  experiments.  The 
calibration  taken  after  the  experiments  is  shown  in  Figure  2.1.  The  thrust  stand  zero 
thrust  point  shifted  significantly  during  the  experiment.  Therefore,  the  actual  thrust  can 
be  bound  using  the  initial  and  final  zero  points  and  the  final  calibration.  The  error  bars 
on  the  measurements  shown  in  Figure  2.5  reflect  the  spread  in  thrust  values  estimated, 
which  is  a  result  of  the  thrust  stand  zero  position  shift. 


30 


100  150  200  250  300  350  400 

Vd(V) 


Figure  2.4  I-V  curves  from  the  PEPL-70  at  different  xenon  volumetric  flow  rates. 


Vd(V)  V  (seem)  Id(A)  T(mN)  rieff(%) 

300  23.4  1.8  26.1  27 

300  27.4  2.4  33.8  29 

300  34.0  3.2  45.4  32 

300  38.6  3.8  55.3  35 

300  43.0  4.6  65.4  36 


1  2  3  4  5 

Id  (A) 


Figure  2.5  Thruster  performance  data  with  varying  propellant  flow  rate  and  a 
constant  discharge  voltage,  300V.  The  numbers  on  the  markers  on  the  graph 
correspond  to  the  total  propellant  flow  rate  (seem).  The  table  above  the  graph 
shows  the  discharge  current  also. 
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Eqn.  2.1 


T2 

n  =  N  — 
**  2  mIdVd 


2.2.2  Xenon  Pressure  and  Ion  Current  Density  in  Cathode  Region 

Xenon  pressure  and  ion  current  density  were  measured  at  several  thruster 
operating  points.  The  NPF  probe  showed  that  the  pressure  in  the  cathode  region  of  a  Hall 
thruster  at  each  operating  point  is  primarily  attributed  to  the  neutral  particles;  retarding 
the  ions  had  virtually  no  effect  on  the  measured  pressures.  Table  2.1  shows  data  acquired 
with  the  Faraday  probe  flush  with  the  exit  plane  of  the  thruster.  Corresponding  NPF  data 
are  not  available.  Table  2.2  shows  data  taken  with  both  probes  positioned  3  mm  behind 
the  physical  exit  plane  of  the  thruster. 


Table  2.1  Ion  current  densities  measured  by  the  Faraday  probe  flush  with  the  exit 
plane  of  the  thruster.  Discharge  voltage  and  current,  electromagnet  inner,  IMI,  and 
outer,  IMO,  currents,  and  Xe  pressures  are  also  given.  _ 


V  (seem) 

I„  (A) 

V„(V) 

IUA) 

Imo(A) 

Pmik  (Torr) 

J  (mA/cm2) 

23.8 

1.9 

259 

19 

8 

mssm 

0.010 

23.8 

1.9 

300 

19 

8 

■KSEB 

0.010 

28.6 

2.7 

298 

18 

5 

1.3xl0'5 

0.013 

33.5 

3.5 

298 

18 

5 

1.5x1 0"5 

0.018 

38.0 

3.5 

300 

28 

15 

1.7xl0'5 

0.014 

42.7 

4.5 

300 

26 

7 

1.9xl0"5 

0.019 

43.1 

4.9 

300 

27 

0 

1.9xl0"5 

0.020 

Table  2.2  Neutral  Xe  particle  pressures  and  ion  current  densities  measured  by  the 
probes  3  mm  behind  the  physical  exit  plane  of  the  thruster. 


V  (seem) 

h  (A) 

V„(V) 

IUA) 

WA) 

P«*  (Torr) 

Pnpf  (Torr) 

J  (mA/cm2) 

18.7 

1.3 

260 

25 

12 

8.0x10-* 

2.1xl0"5 

0.001 

19.1 

1.3 

280 

25 

12 

8.3x10-* 

2.1xl0"5 

0.001 

19.4 

1.3 

302 

25 

12 

8.3x10-* 

2.1xl0"5 

0.001 

23.5 

2.0 

200 

27 

7 

9.8x10-* 

8.9x10-* 

0.003 

23.4 

1.8 

300 

27 

7 

l.OxlO"5 

9.1x10-* 

0.002 

27.4 

2.4 

300 

27 

7 

1.2xl0"5 

9.1x10-* 

0.003 

34.0 

3.1 

300 

27 

8 

1.4xl0"5 

1.3xl0"3 

0.005 

38.6 

3.8 

300 

27 

8 

1.9xl0"5 

1.7xl0"5 

0.006 

43.0 

4.6 

300 

26 

8 

1.9x1 0"5 

1.6xl0"5 

0.008 
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Data  in  Table  2.1  and  Table  2.2  show  that  xenon  pressure  measurements  by  the 
NPF  probe  near  the  thruster  could  be  higher  or  lower  than  tank  pressure.  This 
relationship  depends  on  the  ionization  fraction  of  the  propellant  exiting  the  thruster,  and 
the  ion  optics  of  the  thruster,  which  depend  on  the  discharge  voltage  and  currents  in  the 
electromagnets.  The  ion  optics  were  typically  optimized  for  this  thruster  each 
propellant  flow  rate  and  discharge  voltage  by  adjusting  the  currents  in  the 
electromagnets. 

The  measurements  obtained  in  the  experiments  discussed  in  this  chapter  will  be 
used  in  the  following  chapters  of  the  thesis  to  determine  the  effect  of  this  thruster 
environment  on  cathode  performance  and  lifetime.  The  pressure  measured  by  the  NPF 
probe  is  strongly  influenced  by  vacuum  chamber  pressure.  Therefore,  lower  pressures 
are  expected  in  the  cathode  region  in  lower  pressure  environments.  Ion  current  densities 
of  0.02  mA/cm2,  and  local  pressures  up  to  2xl0'5  Torr  were  used  in  the  models  and 
experimental  performance  evaluations  to  estimate  lower  limits  on  the  cathode  lifetime. 
These  measurements  must  be  conducted  on  any  thruster  that  could  use  FE  cathodes  to 
assess  their  compatibility. 
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CHAPTER  3 


FIELD  EMISSION  CATHODE  EXPERIMENTS 

FE  cathodes  are  typically  used  in  UHV  environments  where  lifetimes  of 
thousands  of  hours  have  been  demonstrated.  The  experimental  results  discussed  in  the 
previous  chapter  show  how  hostile  the  thruster  environment  is,  with  pressures  several 
orders  of  magnitude  higher  than  the  UHV  environments  preferred  for  these  cathodes.  In 
this  environment,  FE  cathodes  can  become  contaminated  and  can  be  sputtered  by  ions 
which  are  created  by  the  electrons  emitted  from  the  cathode  and  CEX  ions  originating 
near  the  thruster.  Cathode  environment  characterization  was  conducted  only  with  xenon 
propellant;  however,  there  is  also  interest  in  using  these  cathodes  in  oxygen 
environments,  as  discussed  in  Chapter  1.  Preliminary  experiments  were  conducted  in 
hostile  environments  at  the  Linfield  Research  Institute  (LRI)  to  determine  if  single  FE 
cathode  tips  can  operate  in  argon  and  oxygen  environments  with  pressures  as  high  as 
~10'5  Torr.  The  results  of  these  preliminary  experiments  were  promising;  therefore,  the 
investigation  into  the  compatibility  of  these  two  technologies  continues.  FEA  and  carbon 
film  FE  cathodes  were  later  obtained  for  performance  characterization  in  UHV 
environments  and  performance  degradation  experiments  in  xenon  environments. 
Experiments  were  conducted  on  single  FE  tip,  FEA,  and  carbon  film  FE  cathodes  to 
study  the  emitting  area  of  single  tips,  emission  stability  in  higher  pressure  environments, 
changes  in  work  function  due  to  contamination  by  the  environment,  and  the  effects  of 
gate  and  anode  electrode  voltages  on  cathode  sputtering.  The  results  of  these 
experiments  are  discussed  in  this  chapter. 
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3.1  Carbide  Single  Tip  Cathode  Experiments 


Single  carbide  FE  cathode  testing  conduced  at  LRI  focused  on  determining  if  FE 
cathodes  could  operate  in  oxygen  and  argon  environments  at  pressures  as  high  as  10'5 
Torn  This  investigation  also  focused  on  gaining  an  understanding  of  the  effect  of  argon 
and  oxygen  exposures  on  the  performance  of  single  carbide  FE  tips. 

3.1.1  Cathodes 

Single  HfC  and  ZrC  FE  cathodes  were  fabricated  and  tested  at  the  Linfield 
Research  Institute  (LRI).  Carbide  FE  cathode  work  has  been  the  primary  focus  at  LRI 
for  the  past  several  years.  FEA  cathodes  have  been  coated  with  carbides  at  LRI,  however 
most  of  the  cathode  fabrication  work  at  LRI  focuses  on  single  tips  that  are  fabricated 
from  commercially  available  sintered  carbide  rods.  The  rods  are  crystallized  at  LRI 
through  a  process  of  floating  zone  arc  refinement61  with  carbon  to  metal  ratios  of  the 
single  crystal  specimen  typically  less  than  one.  The  zone  refined  rods  are  ground  to  0.75 
mm  diameters  and  then  sandwiched  between  two  0.5  mm  pyrolytic  carbon  blocks  in  a 
Vogel-type  mount.  The  rod  and  blocks  are  in  spring  tension  between  two  Mo-Re  posts 
on  .  ceramic  base.  A  single  tip  carbide  cathode  in  the  Vogel-type  mount  is  shown  in 
Figure  3.1.  The  FE  cathode  is  then  formed  by  electrochemical  polishing  of  the  rod  into  a 
sharp  tip.  This  fabrication  process  produces  cathodes  of  varying  tip  radii  and  cone 
angles.  After  fabrication,  the  cathodes  are  imaged  in  a  scanning  electron  microscope 
(SEM)  to  determine  the  cathode  structure.  In  UHV  environments  the  cathodes  are 
‘flashed’  to  temperatures  greater  than  1900  K  by  resistively  heating  the  cathode  base. 

This  process  removes  oxygen  and  other  impurities  from  the  cathode  surface.  Field 
desorption  is  commonly  used  to  smooth  an  uneven  tip  surface  by  positively  biasing  the 
tip  to  ~10kV. 
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Figure  3.1  A  scanning  electron  microscope  (SEM)  micrograph  of  a  single  carbide 
cathode  in  a  Vogel-type  mount. 


3.1.2  Experimental  Apparatus 

Carbide  cathode  testing  was  performed  at  LRI  in  UHV  facilities.  Two  vacuum 
chambers  were  primarily  used  in  the  experiments.  They  employed  mechanical,  turbo, 
and  ionization  pumps  to  attain  pressures  as  low  as  10'10  Torr.  The  vacuum  chambers 
employ  Field  Emission  Microscopes  (FEM)  which  have  phosphor  coated  glass  with  a  tin 
oxide  conductive  coating  that  serves  as  the  anode.  The  cathodes  are  positioned  a  few 
centimeters  behind  the  ancuc;  therefore,  several  kilovolts  are  required  from  the  anode 
power  supply  to  provide  the  electric  fields  at  the  cathode  tips  necessary  for  field 
emission.  As  the  anode  collects  the  current  emitted  by  the  cathode,  the  film 
phosphoresces  from  electron  bombardment.  The  image  observed  on  the  phosphor  coated 
glass  is  a  highly  magnified  representation  of  the  emission  pattern  on  the  cathode  surface. 
The  FEM  configuration  is  illustrated  in  Figure  3.2  showing  a  pattern  that  has  been 

observed.  An  electrical  schematic  for  the  FEM  is  shown  in  Figure  3.3.  A  Keithly  485 
picoammeter  was  used  to  measure  emission  current.  A  1-20  MQ  resistor  was  used  to 

limit  the  cathode  current  to  prevent  cathode  destruction  as  a  result  of  arcing  between  the 

tip  and  anode. 
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Figure  3.2  An  illustration  of  how  the  field  emission  microscope  (FEM)  was  used  to 
determine  the  emitting  area  of  the  tip  from  the  angle,  O,  which  can  be  calculated 
from  the  size  of  the  emission  pattern,  d,  on  the  screen  at  a  distance  x  from  the 
emitter. 


Figure  3.3  Electrical  schematic  of  the  FEM. 


3.1.3  Experimental  Methods  and  Results 

Twenty-seven  HfC  and  ZrC  cathodes  were  fabricated  and  tested  in  FEMs  at  LRI. 
SEMs  were  used  to  image  the  cathodes  after  they  were  fabricated  to  determine  the  cone 
angles  of  the  cathodes  and  radii  of  curvature  of  the  tips.  The  tips  were  also  imaged  after 
the  experiments  to  quantify  changes  in  tip  structure.  The  resolution  of  the  SEM  was 
approximately  10  nm.  The  radii  of  the  tips  fabricated  and  tested  were  typically  much 
larger  than  10  nm.  The  FEM  was  used  to  observe  emission  patterns  from  the  cathodes. 
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Because  different  crystal  planes  in  the  carbides  have  different  work  functions,  electron 
emission  is  not  uniform  from  the  cathode  tip.  Instead,  the  majority  of  the  emission 
current  originates  from  the  low  work  function  planes.  The  pattern  observed  on  the  FEM 
screen  indicates  where  the  tip  is  centered  in  the  crystal  structure  and  which  planes  are 
most  strongly  emitting.  Organized  patterns  are  only  observed  for  clean  cathodes.  The 

size  of  the  pattern  is  indicative  of  the  area  of  the  cathode  that  is  emitting,  wnich  is  a 
function  of  half-cone  emission  angle,  a,  and  the  tip  radius,  r,.  FEMs  were  used  to 

measure  a,  with  tan(c)=d/2x,  as  shown  in  Figure  3.2.  Current-voltage  (I-V)  traces  give 

insight  into  the  emission  stability,  and  changes  in  work  function  and  tip  radius. 

The  SEM  and  FEM  were  used  to  gain  a  better  understanding  of  the  cathode 
dimensions  and  characteristics  required  for  performance  modeling.  SEM  images  of 
cathodes  used  in  this  investigation  showed  that  the  average  half-cone  angle  of  the 
emitting  tip  is  0.26  radians  (16°).  Tip  radii,  measured  in  the  SEM,  ranged  from  40-270 
nm.  The  half  cone  emission  angle,  a,  ranged  from  8-44°,  averaging  28°  for  eight 

cathodes  tested.  FEM  images  show  that  electron  emission  from  the  area,  defined  by  the 
tip  radius  and  half  cone  emitting  angle,  was  not  uniform.  In  fact,  only  a  small  fraction  of 
that  area  is  emitting.  Figure  3.4  shows  an  emission  pattern  observed  in  the  FEM.  The 
patterns  shown  in  Figure  3.4  are  representative  of  a  fairly  clean  tip,  while  the  patterns 
shown  in  Figure  3.5  are  representative  of  an  emission  pattern  from  a  cathode  before 
flashing  it  or  performing  field  desorption.  Organized  patterns  were  never  observed  from 
some  cathodes,  even  after  flash  cleaning  and  field  desorption. 
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Figure  3.4  Emission  patterns  observed  with  the  FEM  from  fairly  clean  carbide 
cathodes. 
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Figure  3.5  Emission  patterns  observed  with  the  FEM  from  damaged  or 
contaminated  carbide. 

The  effect  of  contamination  and  ion  bombardment  on  cathode  performance  was 
investigated  in  argon,  oxygen,  and  air  environments  with  promising  results.  These 
experiments  were  conducted  in  much  more  hostile  environments  than  the  cathode  will  be 
subjected  to  when  used  with  propulsion  systems.  Operating  voltages  were  several 
kilovolts  in  the  FEM  because  the  distances  between  the  tips  and  anodes  were  ~1  cm; 
however,  they  will  be  only  a  few  tens  of  volts  in  field  emitter  array  configurations  when 


used  with  Hall  thrusters.  The  cathodes  demonstrated  more  stable  performance  in  UHV 
environments  than  in  higher  pressure  environments.  Using  a  ~1M£2  resistor  behind  the 

tip  did  prevent  catastrophic  arcing  between  the  tip  and  anode  during  unstable  operation  in 
argon  and  oxygen  environments;  however,  it  did  not  regulate  the  current  enough  in  every 

experiment  to  prevent  tip  melting.  Experimental  results  showed  that  cathodes  operating 
at  ~1  p  A  at  3xl0'8  Torr  were  rarely  stable;  the  current  typically  slowly  degraded  with 

time.  However,  in  one  experiment  the  cathode  stability  improved  significantly  when  the 
current  was  increased  up  to  80  pA  at  3  kV.  Electron  emission  from  this  cathode  was 

incredibly  stable  for  more  than  one  hour.  In  other  experiments,  the  cathodes  were 
exposed  to  oxygen  to  determine  if  oxygen  radicals  and  ions  would  sputter  and 

contaminate  the  cathodes  to  quickly  destroy  them.  The  HfC  cathode  was  conditioned  at 
9xlO'10  Torr  by  running  it  and  flashing  it.  With  the  cathode  emitting  1  pA,  the  oxygen 

pressure  was  increased  up  to -KT6  Torr.  Some  current  reduction  was  observed;  however 
the  emission  pattern  on  the  FEM  screen  improved,  showing  enhanced  emission  from  sites 
not  previously  detected  in  the  FEM.  The  cathode  continued  to  operate  when  pressures 
reached  10  mTorr  by  accident.  At  this  pressure  the  cathode  was  voluntarily  turned  off, 
b!’v  was  not  destroyed. 

In  another  set  of  experiments,  a  ZrC  cathode  was  exposed  to  3xl0'5  Torr  of  argon 
for  ~0.5  hr.  The  cathode  was  turned  on  at  3xlO‘10  Torr.  It  was  flashed  several  times; 
however,  the  pattern  observed  with  the  FEM  was  not  representative  of  a  clean  emitter. 
After  the  argon  pressure  was  increased  up  to  3x1  O'5  Torr  with  -20  pA  and  -6  kV,  a 

pattern  started  to  develop,  as  shown  in  Figure  3.6a.  After  several  minutes,  it  is  believed 
that  a  glow  discharge  developed.  The  emission  pattern  grew  into  a  solid  disk,  as  shown 
in  Figure  3.6b,  and  the  current  jumped  up  to  a  few  hundred  microamperes.  The  argon  gas 
was  then  evacuated  from  the  vacuum  chamber  and  the  pressure  returned  to  10‘9Torr. 
After  the  exposure,  returning  to  UHV,  and  flashing  the  cathode,  the  emission  current  was 
incredibly  stable  and  the  emission  pattern  was  representative  of  a  clean  cathode.  The 


40 


post-exposure  FEM  image  is  shown  in  Figure  3.6c.  Post-exposure  examination  of  the 
cathode  with  the  SEM  showed  no  measurable  damage  by  ion  bombardment  from 
exposure  to  argon  ~10 5  Torr  and  6  kV.  I-V  curves  taken  before  and  after  the  cathode 
was  exposed  to  argon  while  operating  are  shown  in  Figure  3.7.  This  figure  also  shows  a 
Fowler-Nordheim  (F-N)  plot  of  the  Current-Voltage  (I-V)  data,  on  which  ln(I/V2)  is 
plotted  vs.  1/V.  The  argon  ioit  bombardment  actually  cleaned  the  cathode  and  improved 
the  emission  stability,  as  shown  by  the  increased  linearity  of  the  F-N  curves.  The  change 
in  slope  of  the  F-N  curves  shown  in  Figure  3.7  indicates  a  change  in  tip  radius.  Before 
the  exposure,  emission  could  have  been  originating  from  nanoprotrusions  on  the  cathode 
tip.  During  the  exposure  the  surface  was  smoothed  and  cleaned  to  increase  the  effective 
tip  radius.  The  emission  current  was  reduced  at  the  low-voltage  end  of  the  I-V  trace,  and 
increased  at  the  high-voltage  end  of  the  trace. 


Figure  3.6  Evolution  of  *ue  tip  emission  pattern  observed  with  the  FEM  a,b)  during 
and  c)  after  the  exposure  to  argon. 
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Figure  3.7  Fowler-Nordheim  (F-N)  (left)  and  I-V  curves  (right)  obtained  from  the 
ZrC  cathode  before  and  after  exposure  to  argon. 


HfC  and  ZrC  had  already  demonstrated  superior  performance  over  Mo  and  Si 
cathodes.  Carbide  cathodes  can  emit  higher  currents  at  the  same  voltage  because  of 
lower  work  functions.62’63,64’65  In  this  investigation  they  demonstrated  incredible 
robustness  in  oxygen  and  argon  environments  and  resistance  to  oxygen  poisoning. 

Resistance  to  oxygen  poisoning  is  critical  to  the  successful  integration  of  FE 
cathodes  and  oxygen  fueled  propulsion  systems  tethered  spacecraft.  This  cathode 
characteristic  is  also  valuable  when  they  are  used  with  xenon  thrusters  because  the 
cathodes  could  be  frequently  exposed  to  air  during  ground  testing,  and  tum-on  will  occur 
in  air  partial  pressures  between  10’7  and  10 8  Torr.  HfC  and  ZrC  materials  are  preferred 
for  thruster  and  tether  applications;  however,  carbide  FEA  cathodes  were  not  available 
for  testing  in  this  investigation. 

3.2  Mo  and  Si  Field  Emission  Array  Cathode  Experiments 

Field  emission  array  cathodes  with  arrays  of  more  than  10,000  tips  were  used  to 
determine  the  effect  of  xenon  exposure  on  Mo  and  Si  FEA  cathode  performance.  The 
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advantage  of  these  microtip  cathode  arrays  is  that  the  operating  voltages  can  be  less  than 
100  V  to  get  1  mA  of  current.  The  objective  of  these  experiments  is  to  determine  if 
xenon  adsorption  will  affect  the  work  function  of  the  cathodes  and  if  the  cathodes  can  be 
operated  at  low  enough  voltages  to  avoid  sputter  damaging  the  microtip  structures.  These 
experiments  also  resulted  in  estimates  for  energy  threshold,  E*.  for  sputtering  of  Mo  and 
Si  targets  by  Xe  ions  which  are  better  than  those  previously  available.  Experimental 
results  discussed  in  the  previous  chapter  showed  that  the  xenon  pressure  in  the  cathode 
environment  of  a  Hall  thruster  at  1.4  kW  is  around  2x1  O'5  Torr  at  1.4  kW.  Therefore, 
most  of  the  performance  evaluations  were  conducted  at  this  pressure.  The  performance 
of  these  cathodes  was  not  optimized  in  these  experiments  because  of  limited  time  and 
cathode  availability.  Much  better  performance  is  expected  if  the  cathodes  are  tested  in 
cleaner  environments  and  operated  at  higher  voltages.  These  experiments  focused  on 
low-voltage  operation  where  emission  currents  were  typically  less  than  1  pA. 

3.2.1  Cathodes 

Silicon  cathodes  from  MCNC 

The  MCNC  FEA  cathodes  used  in  these  experiments  employed  J  f  ,000  silicon  tips 
with  2  nm  aperture  diameters  on  10  nm  centers.  Four  cathode  arrays  were  fabricated  on  a 

single  chip  with  four  gate  electrodes  and  a  cathode  base  common  to  each  of  the  arrays. 
The  chip  was  mounted  to  a  T0-5  head.  Each  tip  in  the  arrays  used  in  the  experiments  is 
similar  to  the  tip  shown  in  Figure  3.8. 

The  cathodes  are  fabricated  at  MCNC  using  a  self-aligning  tip-on-post  and  gate 
electrode  process.66  The  process  begins  with  deposition  and  then  reactive  ion  etching 
(RTF.)  of  a  nitride  film  to  leave  circular  caps  on  a  silicon  wafer.  The  caps  are  used  as 
masks  for  etching  silicon  to  form  silicon  posts  and  tips.  Silicon  posts  are  formed  using  an 
anisotropic  RIE  with  a  mixture  of  BC13,  Cl2, 02,  and  He  gases.  The  wafers  are  then 
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oxidized  to  sharpen  the  posts  into  tips,  with  the  nitride  caps  in  place.  Si02  is  then 
deposited  onto  the  structure  with  a  protective  photoresist  layer  around  the  cone  structures. 
Chromium,  tantalum,  and  platinum  gate  electrode  film  is  then  deposited.  The  liftoff 
material  is  removed  and  a  self-aligned  gate  electrode  remains.  The  nitride  caps  are 
removed  in  a  buffered  oxide  etch,  and  a  hemispherical  cavity  remains  around  the  post,  as 
shown  in  Figure  3.8. 


Figure  3.8  A  MCNC  Si  field  emission  array  single  tip  configuration  (courtesy  of 
MCNC). 

Molybdenum  cathodes  from  SRI 

The  molybdenum  FEA  cathode,  provided  by  SRI  International,  consisted  of  an 
array  of  50,000  tips  with  0.9  urn  gate  aperture  diameters  on  4  nm  centers  fabricated  on  a 

2000  Ohm-cm  silicon  wafer.  The  high  resistivity  wafer  protects  the  cathodes  against  arc 
damage  by  effectively  providing  an  isolated  resistor  in  series  with  each  emitter.  Figure 
3.9  shows  a  FEA  cathode  from  SRI  that  was  cracked  in  liquid  nitrogen  and  imaged  in  a 
SEM.  The  cathode  shown  has  been  damaged  by  xenon  ion  bombardment.  An 
undamaged  cone  is  much  sharper,  with  the  tip  flush  with  the  top  surface  of  the  gate 
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electrode.  The  structure  of  the  conical  SRI  cathode  is  significantly  different  from  the 
MCNC  cathode,  which  has  a  tip-on-post  structure. 

The  molybdenum  cathode  used  in  these  experiments  was  fabricated  at  SRI  using  a 

process  in  which  the  gate  electrode  is  deposited  before  the  tip  structure.  First  the  silicon 
w?fer  was  covered  with  a  1  pm  thick  insulating  layer  of  thermally  grown  Si02.  A 
molybdenum  film  with  a  0.25  pm  thickness  Was  deposited  on  the  insulating  layer.  An 
array  of  50,000  apertures  with  0.9  pm  diameters  on  4  pm  centers  were  then  patterned  on 

the  Mo  film  with  electron  beam  lithography  and  then  etched.  Cavities  in  the  Si02  were 
then  etched  using  the  gate  electrode  as  the  etch  mask.  A  lift-off  layer  was  then  deposited 
on  the  gate  electrode.  Mo  was  then  deposited  on  the  film  forming  a  cone  structure  in  the 
cavity  as  the  gate  aperture  closed.  The  excess  Mo  was  then  lifted  off  and  the  cone 
structure  remained  in  the  cavity  as  shown  in  Figure  3.9.  The  cathodes  are  then  cleaned  in 
a  H+Ne  plasma. 
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Figure  3.9  A  SRI  International  Mo  field  emission  array  cathode  single  tip  and 
array. 
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3.2.2  Experimental  Apparatus 


Silicon  and  molybdenum  FEA  cathode  testing  was  performed  at  the  JPL.  The  test 
facility  consisted  of  a  0.0664  m3  vacuum  chamber  that  was  pumped  by  turbomolecular, 
ion,  and  sublimation  pumps  to  pressures  as  low  as  7xl0‘n  Torr.  The  vacuum  facility  is 
shown  in  Figure  3.10.  The  cathodes  were  mounted  to  electrical  feedthrus  on  conflat 
flanges.  The  cathode  test  flange  is  shown  in  Figure  3.11.  Xenon  pressure  was  increased 
using  a  Granville-Phillips  leak  valve,  “The  Leak”.  The  cathodes  and  the  ionization  gauge 
were  equidistant  from  the  xenon  inlet  orifice  and  vacuum  pumps  to  minimize  pressure 
gradients  between  the  cathodes  and  pressure  gauge.  Three  Keithley  480  picoammeters 
were  used  to  measure  the  current  emitted  through  the  base  and  collected  by  the  gate 
electrode  and  anode.  The  picoammeters  on  the  anode  and  gate  electrodes  floated 
electrically  above  ground  to  minimize  the  measurement  of  stray  currents.  Pressure  was 
measured  using  a  Varian  UHV-24  nude  ionization  gauge.  Pressure,  current,  and  voltage 
signals  were  recorded  using  a  National  Instruments  Lab  VIEW  data  acquisition  system. 
The  electrical  configuration  employed  during  the  experiments  is  shown  in  Figure  3.12.  A 
10  Nto  resistor  was  used  on  the  gate  electrode  for  microtip  cathode  testing.  This  resistor 

was  used  to  significantly  reduce  Vg  when  several  microamperes  of  current  were  collected 
by  the  gate  electrode,  to  prevent  arcing  between  the  tip  and  gate  electrode. 
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Figure  3.10  The  ultra-high  vacuum  (UHV)  facility  at  the  Jet  Propulsion 
Laboratory. 


Figure  3.11  The  field  emission  cathode  test  flange  in  UHV  facility. 
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Figure  3.12  Electrical  schematic  for  field  emission  cathode  testing. 


3.2.3  Si  FEA  Cathode  Experimental  Methods  and  Results 

The  first  experiments  investigated  the  effect  of  xenon  exposure  on  the  work 
function  of  a  MCNC  Si  FEA  cathode  with  16,000  tips.  The  cathode  was  exposed  to 
xenon  at  2x1  O'4  Torr  for  1  hour  to  provide  a  dose  of  xenon  to  the  surface  of  the  cathode 
greater  than  it  experienced  during  the  exposure  tests  at  2x10 5  Torr  with  the  cathode 
operating.  Data  obtained  during  this  and  repeated  experiments  show  that  xenon  coverage 
did  not  change  the  work  function  of  the  silicon  cathode  surface.  These  data  are  shown  in 
Figure  3.13. 
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Figure  3.13  I-V  curves  taken  before  and  after  the  Si  cathode  was  exposed  to  Xe. 

The  objective  of  the  next  experiments  was  to  determine  if  the  cathode  could 
operate  at  low  enough  voltages  to  avoid  being  sputter  damaged.  An  experiment  similar 
to  those  conducted  at  MCNC  was  conducted  at  JPL  where  the  cathode  was  exposed  to 
7x1  O'6  Torr  of  xenon  with  an  initial  current  of  20  pA.  The  gate  voltage  in  this  experiment 

was  86  V  and  the  anode  voltage  was  400  V.  The  response  of  the  cathode  current  to  the 
increase  in  pressure  is  shown  in  Figure  3.14.  During  the  2  hours  of  operation  in  UHV 
after  the  exposure  test,  the  cathode  current  continued  to  decay.  The  experiment  was 
repeated  on  another  cathode  with  the  same  configuration.  During  this  exposure  the  initial 
current  was  22  pA  with  75  V  on  the  gate  electrode  and  500  V  on  the  anode.  The  same 

results  were  obtained;  during  the  10  hr.  at  UHV  that  followed  the  second  exposure,  the 
current  never  recovered  to  the  original  value.  The  results  of  these  experiments  suggest 
that  the  tips  were  damaged  from  ion  bombardment.  Experiments  conducted  with  anode 
voltages  below  90  V  also  showed  performance  degradation  due  to  sputtering;  therefore, 
high-energy  ions  generated  between  the  gate  and  anode  electrodes  are  not  responsible  for 
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the  tip  sputtering  observed.  A  model  for  estimating  the  was  developed  by  Yamamura 
et  al.67  This  model  was  used  to  predict  that  the  energy  threshold  for  xenon  ions  sputtering 
a  silicon  target  is  91  eV.  Because  damage  was  done  to  the  cathode  at  anode  and  gate 
voltages  below  that  value,  the  true  threshold  value  must  be  lower  than  91  eV. 
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Figure  3.14  Si  cathode  current  response  to  an  increase  in  Xe  pressure  with  Vg  at  86 
V  and  Va  at  400  V. 


Experiments  were  then  conducted  at  lower  gate  voltages  to  determine  the  xenon 
ion  energy  threshold  for  sputtering  silicon.  Irreversible  performance  degradation  was 
observed  during  exposures  to  xenon  at  2xlO'5Torr  with  gate  voltages  of  75  (Va  at  100  V), 
70  (Va  at  70  V)  and  60  V  (Va  at  80  V).  Typically,  the  current  dropped  in  half  during  the 
first  hour  of  the  exposure.  Data  acquired  from  the  xenon  exposure  experiment  at  2xl0'5 
Torr  with  Va  and  Vg  at  70  V  is  shown  in  Figure  3.15.  An  initial  increase  in  current  was 
commonly  observed  when  xenon  was  introduced.  This  response  could  be  attributed  to 
cleaning  of  the  cathode  by  ion  bombardment  before  the  ions  begin  to  damage  the 
cathode.  Cathode  cleaning  was  observed  in  the  FEM  when  carbide  cathodes  were 
exposed  to  argon  during  operation  at  LRI.  With  50  V  applied  to  the  gate  electrode  no 
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performance  degradation  was  observed  during  a  one  hour  exposure  to  2x10  Ton-  of 
xenon  with  70  V  applied  to  the  anode;  the  current  actually  increased  slightly.  The  results 
of  this  experiment  are  presented  in  Figure  3.16.  During  other  experiments  with  MCNC 
cathodes,  the  current  through  the  gate  electrodes  was  usually  less  than  a  factor  of  100 
lower  than  the  anode  current.  The  efficiency  of  this  cathode  was  poor;  about  50%  of  the 
current  was  measured  through  the  gate  electrode.  After  the  exposure,  the  emitted  current 
increased  while  the  efficiency  decreased. 
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Figure  3.15  Si  cathode  current  response  to  2xl0  5  Torr  of  Xe  with  Vg  at  70  V  and  the 
Va  at  70  V. 
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Figure  3.16  MCNC  Si  cathode  current  response  to  a  1  hour  Xe  exposure  with  Vg  at 
50  V  and  V,  at  70  V. 


These  MCNC  cathodes  operated  for  several  hours  without  catastrophic  arcing  in 
xenon  environments.  The  results  of  these  experiments  showed  that  xenon  adsorption 
does  not  change  the  cathode  work  function.  Performance  degradation  occurred  in  xenon 
environments  as  a  result  of  the  microtip  structure  being  sputtered  by  ion  bombardment. 
Experimental  results  showed  that  when  the  cathodes  were  operated  at  low  enough 
voltages,  ~50  V,  ion  energies  were  not  high  enough  to  cause  significant  sputter  damage  in 
one  hour.  Silicon  cathodes  emitting  100  mA  at  50  V  should  satisfy  current  and  lifetime 

limitation  requirements  necessary  for  the  successful  integration  of  Hall  thruster  and  FE 
cathode  technologies.  A  current  of  2.5  pA  has  been  measured  from  a  single  Si  microtip 

structure  with  the  gate  voltage  at  25  V.68  Therefore,  an  array  of  less  than  100,000  tips 
could  provide  100  mA  of  current  to  a  Hall  thruster  if  2.5  p  A/tip  can  be  extracted  from  the 

array  at  25  V. 


52 


P  (Torr) 


3.2.4  Mo  FEA  Cathode  Experimental  Methods  and  Results 

The  first  experiments,  conducted  on  the  Mo  cathode  from  SRI  International  with 
50,000  tips,  investigated  the  effect  of  xenon  adsorption  on  the  work  function  of  the 
cathode.  The  cathode  was  exposed  to  2x10“*  Torr  of  xenon  for  one  hour  to  provide  a  dose 
of  xenon  to  the  surface  of  the  cathode  greater  than  it  experienced  during  the  exposure 
tests  at  2xl0"5  Torr  with  the  cathode  operating.  I-V  traces  were  taken  at  UHV 
immediately  before  and  after  the  exposure;  these  data  are  shown  in  Figure  3.17.  The  I-V 
trace  taken  after  the  exposure  showed  slightly  better  performance  than  before  the 
exposure.  The  results  of  this  experiment  suggest  that  the  work  function  did  not  increase 
due  to  the  adsorption  of  xenon  or  any  contaminants  in  the  xenon  supply.  Any 
performance  degradation  observed  during  xenon  exposure  can,  therefore,  be  attributed  to 
changes  in  the  geometrical  parameters  of  the  cathode. 


Figure  3.17  I-V  data  taken  before  and  after  the  Mo  cathode  was  exposed  to  xenon. 

The  next  experiments  focused  on  studying  microtip  sputtering  by  ion 
bombardment  in  xenon  environments.  A  xenon  exposure  test,  with  operating  conditions 
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similar  to  the  operating  conditions  during  experiments  conducted  at  SRI  International 
with  a  gate  electrode  voltage  at  65  V  and  an  anode  voltage  of  80V  with  ~18  |iA  for  one 

hour  in  2x1  O'6  Torr.  Figure  3.18  shows  the  current  response  to  the  increase  in  xenon 
pressure.  The  current  dropped  to  one-third  of  its  original  value  in  one  hour.  No  current 
recovery  was  observed  in  13  hours  at  UHV  that  followed  the  xenon  exposure.  The 
cathode  was  permanently  damaged  from  ion  bombardment. 
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Figure  3.18  Mo  cathode  current  response  to  2xl0'4  Torr  of  xenon  with  Vg  at  65.6  V 
and  Va  at  80  V. 


An  experiment  was  conducted  to  determine  the  gate  voltage  that  corresponded  to 
the  onset  of  cathode  sputter  damage.  The  anode  voltage  in  this  experiment  was  100  V, 
and  the  xenon  pressure  was  2xl0‘5  Torr.  The  gate  voltage  was  increased  in  1  V 
increments  every  10  minutes,  starting  at  50  V.  Performance  decay  caused  by  tip 
sputtering  was  observed  within  10  minutes  when  the  gate  voltage  reached  58  V,  as  shown 
in  Figure  3.19.  With  further  increases  in  gate  electrode  voltage,  the  current  initially 
increases  and  then  decays  before  stabilizing.  Preliminary  experiments  showed  that  a  few 
minutes  are  required  for  stabilization  after  the  changes  in  gate  electrode  voltage.  Figure 
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3.20  shows  the  current  response  to  increases  in  gate  electrode  voltage  in  both  UHV  and 
10'5  Torr  xenon  environments.  This  cathode  had  already  been  exposed  to  xenon 
environments  for  several  hours  and  had  suffered  significant  sputter  damage  before  this 
experiment.  It  is  expected  that  other  cathodes  will  respond  somewhat  differently.  A 
cathode  with  sharper  tips  and  the  same  number  of  emitting  tips,  should  show  performance 
decay  within  ten  minutes  at  even  lower  gate  voltages,  because  the  energy  of  the  ions 
bombarding  the  emitting  area  of  the  tips  will  be  slightly  higher  for  sharper  tips.  This 
effect  is  a  result  of  the  higher  electric  fields  at  the  tips.  The  maximum  energy  of  the  ions 
bombarding  the  emitting  area  of  the  tips  slightly  increases  with  decreasing  tip  radius,  as 
will  be  shown  in  the  following  section. 
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Figure  3.19  The  Mo  cathode  current  response  to  increasing  gate  electrode  voltage 
with  Va  at  100  V,  2xl0'5  Torr  of  Xe. 
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Figure  3.20  A  comparison  of  the  Mo  cathode  current  response  to  increasing  gate 
electrode  voltage  in  2xl0'9  Torr  and  2xl0'5  Torr  of  Xe. 


Several  experiments  were  conducted  with  the  gate  voltage  at  50  V.  The 
maximum  ion  energy  in  this  configuration  is  46  eV,  as  predicted  by  the  models  discussed 
in  the  next  chapter  In  an  experiment  with  the  pressure  at  2xl0'5  Torr  that  lasted  for  5 
hours  with  the  gate  voltage  at  50  V  and  the  anode  voltage  at  60  V,  the  current  initially 
increased  and  then  decayed  monotonically.  The  results  of  this  experiment  are  shown  in 
Figure  3.21.  These  results  imply  that  the  E^,  for  Mo  sputtering  by  Xe  ions  is  below  46 
eV.  Another  experiment  was  conducted  at  the  same  pressure  with  Vg  at  50  V  to  study  the 
effect  of  increasing  anode  voltage.  It  was  anticipated  that  at  a  high  enough  anode  voltage 
the  current  would  begin  to  decay.  The  current  dropped  from  1 10  nA  to  85  nA  in  response 
to  the  initial  increase  in  Xe  pressure  up  to  2xl0’5  Torr,  and  then  increased  with  increasing 
anode  voltage  from  60  V  up  to  130  V  in  10  V  increments  every  10  min.  The  results  from 
this  experiment  are  shown  in  Figure  3.22.  They  show  that  increasing  the  anode  voltage 
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from  V  =V  with  V  ~E,h  can  actually  improve  the  performance  of  the  cathode.  In  an 

a  g  g  in 

experiment  at  the  same  pressure  with  Vg  at  50  V  and  Va  at  100  V,  the  current  only 
increased  during  a  two  hour  xenon  exposure  test  with  constant  voltages.  These  results 
are  shown  in  Figure  3.23. 
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Figure  3.21  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2xl0'5 
Torr  of  xenon  for  5  hours  with  Vg  at  50  V  and  Va  at  60  V. 
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Figure  3.22  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2xl0'5 
Torr  with  varying  Va  and  Vg  at  50  V. 
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Figure  3.23  Mo  cathode  current  response  to  a  Xe  pressure  increase  up  to  2x10  s 
Torr  for  2  hours  with  Vg  at  50  V  and  Va  at  100  V. 


The  performance  improvements  observed  at  low  gate  voltages  and  increasing 
anode  voltage  could  be  attributed  to  increases  in  the  electric  field  at  the  cathode  tips  from 
changes  in  the  tip  microstructure.  An  increase  in  anode  voltage  to  increase  the  electric 
field  at  the  tip  would  result  in  an  initial  increase  before  the  current  equilibrates; 
equilibration  typically  required  less  ~5  min.  During  exposures  with  -100  V  applied  to 
the  anode  and  50  V  applied  to  the  gate  electrode,  the  current  continually  increased.  If  the 
increase  in  anode  voltage  increased  the  energy  of  the  ions  bombarding  the  emitting  area 
of  the  tips,  a  current  decay  would  be  observed  because  the  r,  would  increase  from 
sputtering.  However,  not  only  is  the  emitting  area  being  sputtered,  the  cone  sidewalls  are 
also  being  sputtered  by  ion  bombardment.  While  the  emitting  area  is  being  bombarded 
with  ions  primarily  at  a  normal  incidence  to  the  surface,  the  cone  sidewalls  are  also  being 
bombarded  by  ions  at  smaller  angles.  Sputter  yields  increase  with  decreasing  incident  ion 
angle  from  90°  to  60°.  With  a  uniform  flux  of  ions  to  the  cone  surface,  the  sidewalls  will 
sputter  faster  to  decrease  the  half-cone  angle  and  increase  the  current  at  a  constant 
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voltage.  As  the  anode  voltage  increases,  ion  acceleration  into  the  cathode  tip  is  less 
focused  on  the  cathode  emitting  area  and  more  ions  hit  the  cathode  sidewalls  where  the 
sputter  yield  is  higher.  Cathode  sputtering  by  ion  bombardment  is  discussed  in  more 
detail  in  the  following  chapter. 

The  SRI  cathode  survived  for  21  hours  in  2xl0^-2xl0'5  Torr  of  xenon  without 
catastrophically  arcing.  The  cathode  performance  degraded  during  the  xenon  exposures 
because  the  cathodes  were  sputtered  by  ion  bombardment.  The  results  of  the  xenon 
exposure  experiments  showed  that  xenon  adsorption  does  not  change  the  cathode  work 
function.  They  also  showed  that  it  is  possible  to  avoid  sputter  damage  while  operating  in 
xenon  environments  if  operating  voltages  are  low  enough.  The  experimental  results 
showed  that  the  gate  voltage  should  be  below  50  V  to  avoid  sputter  damaging  the 
emitting  areas  of  the  tips.  It  was  learned  that  anode  voltage  does  affect  cathode 
performance  in  xenon  environments.  Increasing  the  anode  voltage  up  to  100  V  with  the 
gate  electrode  ~50  V  improved  the  cathode  performance,  possibly  by  reducing  cone  half¬ 
angles  or  sharpening  the  emitting  cones.  However,  operating  with  the  anode  voltage  at 
100  V  is  not  favored  because  the  material  sputtered  off  of  the  cone  sidewalls  will  then 
redeposit  on  the  insulator  walls,  eventually  shorting  the  tips  and  gate  electrodes. 

3.3  Carbon  Film  FE  Cathode  Experiments 

A  carbon  film  cathode  from  FEPET  was  used  to  determine  the  effect  of  xenon  ion 
bombardment  on  field  emission  cathodes  with  thin  films  instead  of  microtip  structures. 
These  cathodes  have  low  electron  affinity  and  require  lower  macroscopic  field  strengths 
than  FEA  cathodes.  Electron  emission  from  the  carbon  films  is  not  as  sensitive  to  ion 
bombardment  as  microtip  structures.  The  current  emitted  by  the  microtip  structures  is 
exponentially  sensitive  to  changes  in  tip  radius  from  sputtering  by  ion  bombardment. 
Limited  sputtering  of  carbon  film  may  not  affect  the  emission  current;  these  cathodes  can 
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withstand  more  sputter  damage  without  performance  degradation.  Carbon  film  cathodes 
are  more  sensitive  to  chemical  changes  in  the  surface  of  the  film;  inert  gas  environments 
should  not  affect  the  surface  chemistry. 

Carbon  film  cathodes  can  tolerate  much  higher  pressure  environments  than  micro 
tip  cathodes  when  operating,  and  can  operate  more  safely  at  much  higher  voltages.  The 
disadvantages  to  using  currently  available  carbon  film  cathodes  are  the  high  gate 
electrode  voltages  required  and  low  operating  efficiency.  The  FEPET  cathode  is  not  a 
microintegrated  structure  like  the  microtip  cathodes.  Gate  voltages  of  several  hundred 
volts  are  required  for  microamperes  of  current.  Low  cathode  efficiency  is  also  a  concern; 
the  gate  electrode  collects  up  to  80%  of  the  current.  Large  improvements  in  the  cathode 
efficiency  are  expected  when  a  microintegrated  film  and  gate  structure  is  fabricated. 

In  this  investigation,  the  focus  of  the  experiments  with  a  FEPET  cathode  was  on 
the  response  of  the  current  to  increases  in  xenon  pressure  up  to  2x1  O'6  Torr.  The  pressure 
was  only  increased  to  2xl0'6  Torr  of  xenon  because  of  the  close  proximity  of  the 
ionization  gauge  and  anode  and  the  high  anode  voltage  at  900  V .  Charged  particles  were 
created  in  this  configuration  with  high  pressures  and  were  collected  by  the  electrodes  at  a 
much  higher  rate  than  the  electrons  emitted  by  the  cathode,  making  it  impossible  to 
measure  the  emitted  electron  current.  With  high  operating  voltages  required  by  this 
cathode,  the  rate  of  xenon  ion  bombardment  will  be  much  higher  than  the  cathode  will 
experience  in  the  thruster  environment  where  operating  voltages  will  be  much  lower. 
Therefore,  any  positive  results  in  these  experiments  are  extremely  promising. 

3.3.1  Carbon  Film  FE  Cathode 

The  carbon  film  cathode  used  in  this  investigation  was  fabricated  at  Field 
Emission  Picture  Element  Technology  (FEPET).  The  FEPET  graphitic  and  amorphous 
carbon  films  are  grown  on  a  ceramic  substrate  by  hot  filament  chemical  vapor  deposition 
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(HFCVD).  The  ceramic  substrate  is  specially  treated  to  improve  the  emission  properties 
of  the  film  with  a  proprietary  substance.  Deposition  of  this  substance  can  be  controlled 
when  a  special  emission  pattern  or  emitting  area  is  required.  A  hydrogen  and  methane 
gaseous  mixture  is  introduced  in  a  reactor  chamber  which  is  activated  by  hot  filaments. 
Activated  species  a^  then  deposited  on  the  substrate  surface.  The  filament  and  substrate 
temperatures,  gas  flow  rate,  and  gas  consumption  are  carefully  controlled  to  obtain  films 
with  good  emission  properties.  Post-growth  processing  is  not  required. 

The  cathode  used  in  these  experiments  consisted  of  a  thin  carbon  film  deposited 
on  a  ceramic  substrate  with  a  stainless  steel  extraction  grid  spaced  100  jim  from  the 

cathode.  The  size  of  the  emitting  area  of  the  cathode  is  0.25  cm2.  The  cathode  is 
mounted  on  a  glass  stem  as  shown  in  Figure  3.24. 


Figure  3.24  A  FEPET  carbon  film  cathode  on  a  glass  stem. 

3.3.2  Experimental  Apparatus 

The  experimental  apparatus  is  described  in  section  3.2.2.  The  electrical 
configuration  is  shown  in  Figure  3.12.  Note:  This  figure  shows  a  cathode  with  a  tip 
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structure,  however,  the  FEPET  cathode  does  not  use  an  emitting  cone  structure;  the 
emitting  surface  is  a  thin  film. 

3.3.3  Experimental  Methods  and  Results 

The  efficiency  of  this  film  cathode much  lower  than  the  Spindt-type  cathodes. 
The  gate  electrode  intercepts  the  majority  of  the  current.  About  20  %  of  the  current 
emitted  from  the  carbon  film  was  typically  collected  at  the  anode  during  these 
experiments;  however,  the  efficiency  was  often  much  higher  at  currents  below  200  nA. 
The  performance  of  the  FEPET  cathode  used  in  these  experiments  in  UHV  is  shown  in 
Figure  3.25. 
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Figure  3.25  I-V  traces  showing  the  performance  of  a  FEPET  carbon  film  cathode  in 
a  UHV  environment. 
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The  cathode  demonstrated  impressively  stable  performance  in  the  xenon 
environment  as  shown  in  Figure  3.26.  This  figure  also  shows  that  fluctuations  in  the 
current  can  be  attributed  to  gate  voltage  fluctuations.  Operation  of  the  ionization  gauge 
caused  the  gate  electrode  voltage  to  fluctuate  between  835  V  and  850  V.  This  caused  2 


62 


|xA  fluctuations  in  the  emitted  current.  In  the  high  pressure  environments,  the  ionization 

gauge  affected  the  current  measurements  on  all  of  the  electrodes.  The  ionization  gauge 
was  responsible  for  -250  \iA  to  the  anode  and  15  \iA  to  the  gate  electrode.  With  a  1 

M£2  resistor  on  the  gate  electrode,  this  increase  in  current  changed  the  gate  voltage  by  15 
V.  The  ionization  gauge  was  turned  off  and  on  several  times  during  the  high  pressure 
exposure  experiments  to  record  the  pressure.  The  large  jumps  in  the  gate  voltage,  and 
responding  current  can  be  attributed  to  the  ionization  gauge  current.  The  results  of  these 
experiments  show  that,  even  at  such  high  electrode  voltages,  the  cathode  performance  is 
not  affected  by  ion  bombardment  and  sputtering  of  the  carbon  film  during  -10  hours. 
Two  of  these  exposure  test  were  conducted  at  the  same  operating  conditions.  The  first 
exposure  test  lasted  for  4  hours,  and  the  second  exposure  lasted  for  3  hours.  No 
performance  degradation  was  observed  during  either  of  the  experiments.  I-V  data  taken 
before  and  after  4  hours  of  xenon  exposure  were  identical,  as  shown  in  Figure  3.27. 
These  data  also  show  that  the  cathode  was  not  damaged  during  the  exposures. 
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Figure  3.26  The  cathode  current  response  to  increasing  the  Xe  pressure  up  to 
2x10^  Torr  for  4  hours.  The  anode  potential  was  900  V  above  the  potential  of  the 
carbon  film. 


Figure  3.27  I-V  curves  taken  in  UHV  before  and  after  the  carbon  film  cathode  was 
exposed  to  xenon. 
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Figure  3.28  shows  I-V  data  taken  in  UHV  and  xenon  environments.  The  xenon 
pressure  was  increased  to  2xl0"6  Torr  and  stabilized  while  the  cathode  was  off.  The 
cathode  was  then  slowly  turned  on  in  the  xenon  environment  for  -8  min.  and  then  slowly 
reduced  to  zero  current  to  obtain  two  of  the  data  sets  shown  in  Figure  3.28.  The  cathode 
start-up  in  2x1  O'6  Torr  was  no  different  from  the  start-up  in  UHV  environment.  This 
experiment  was  repeated  with  the  same  results.  This  characteristic  of  the  cathode  is 
extremely  valuable  because  the  base  pressure  of  the  facilities  used  in  Hall  thruster  testing 
is  ~10'7-10'6  Torr;  therefore  cathodes  must  comfortably  start-up  in  this  environment. 


Vg(V) 

Figure  3.28  Carbon  film  cathode  turn-on  and  turn-off  I-V  curves  with  an  anode 
voltage  of  900  V  in  2x10"*  Torr  of  xenon  compared  to  turn-on  I-V  data  in  UHV. 


The  FEPET  FE  cathode  was  operated  for  a  total  of  9  hours  in  2x1  O'6  Torr  xenon 
environment  with  850  V  on  the  gate  electrode  and  900  V  on  the  anode  to  get  ~28  \iA  of 

emission  current.  During  this  time,  the  cathode  showed  no  performance  degradation 
from  the  intense  ion  bombardment.  The  operating  voltages  and  efficiency  of  this  cathode 
must  be  significantly  improved  before  this  cathode  can  be  used  with  ion  and  Hall 
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thrusters;  however  it  is  believed  that  these  objectives  can  be  achieved  with  a 
microfabricated  integrated  gate  and  carbon  film  structure.  At  the  current  state  of  FEPET 
carbon  film  technology,  this  cathode  presents  itself  as  a  strong  candidate  for  ion  beam 
neutralization  with  FEEP  systems,  which  consume  approximately  1/4  W  for  1  mA  at  10 
kV  ion  accelerating  potentials. 
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CHAPTER  4 


MICROTIP  CATHODE  PERFORMANCE  MODELING 


4.1  FEA  Cathode  Performance  Degradation  Model 

To  expedite  FEA  cathode  development  for  EP  systems  and  reduce  development 
costs,  a  FEA  cathode  performance  decay  model  was  used  in  conjunction  with 
experimentation.  A  different  model  must  be  developed  to  predict  the  effect  of  a  plasma 
thruster  environment  on  the  performance  of  thin  film  cathodes.  The  model  discussed  in 
this  chapter  incorporates  3  previously  existing  models  with  slight  modifications  which  . 
predict  cathode  emission  currents,  flux  of  ions  to  the  surface,  and  sputter  yields  of 
cathode  materials  by  xenon  ions,  and  an  additional  model  to  predict  changes  in  tip  radius 
from  ion  bombardment.  With  these  models  and  experimental  results,  the  lifetime  of  Si 
and  Mo  FEA  cathodes  can  be  predicted  in  xenon  environments  of  electric  propulsion 
systems. 

A  flow  diagram  for  this  performance  degradation  model  is  shown  in  Figure  1.7. 

The  inputs  to  this  model  are  I-V  data  obtained  in  UHV,  cathode  work  function,  half-cone 
angle  of  the  microtip  structures,  |3C,  the  energy  threshold  for  sputtering  cathode  material 

by  ion  bombardment,  E^,  gate  aperture  diameters,  rg,  and  number  of  tips,  Ntips.  The  I-V 
data  are  used  in  a  model  developed  by  Jensen  to  estimate  the  effective  tip  radius,  r„  and 
distribution  of  radii  in  the  array.69  The  Jensen  model  is  also  used  to  determine  the 
emission  current  from  an  array.  To  determine  current  changes  in  time  caused  by  changes 
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in  tip  radius,  the  change  in  tip  radius  in  time  must  be  determined.  The  rate  at  which  the 


tip  radius  increases  in  time  is  approximated  by 


Eqn.  4.1 


dri^d 

dt  t. 


where  the  atomic  spacing  in  the  material  is  d  (A),  and  a  monolayer  of  material  is  removed 
from  the  tip  emitting  area  with  an  effective  tip  radius,  r„  in  t,.  (s).  A  model  developed  by 
Brodie  was  used  to  determine  te.  Material  sputter  yields  were  determined  using  the  model 
developed  by  Yamamura  et  al.67  The  tip  radius  increases  in  time  are  approximated  by 
assuming  that  each  layer  removed  from  the  emitting  area  of  the  tip  increases  the  tip 
radius  by  the  atomic  spacing  in  the  tip  material.  Tip  radii  are  represented  statistically  by 
an  effective  tip  radius  and  spread  in  tip  radius  across  an  array.  As.  The  radius  of  the 

sharpest  tips  in  the  array  will  change  at  the  fastest  rate  so  that  the  uniformity  in  the  array 
improves  in  time.  The  changes  in  the  spread  in  tip  radii,  in  time,  shown  in  Eqn.  4.2,  is  a 
function  of  the  change  in  effective  tip  radius  in  time  and  the  maximum  tip  radius  in  the 
array.  I(r(t),As(t»  can  then  be  determined  using  Jensen’s  emission  model  for  FEA 

cathodes. 

Eqn.  4.2  =  1- 

r,(t) 


The  following  three  sections  of  this  chapter  present  the  details  of  the  Jensen 
model  to  predict  cathode  parameters  from  I-V  data  and  cathode  emission  current,  the 
Brodie  model  to  predict  the  removal  rate  of  material  from  the  tip  emission  area,  and  the 
model  developed  by  Yamamura  et  al.67  which  predicts  the  sputter  yields  of  xenon  ions 
bombarding  Si  and  Mo  targets. 


4.1.1  Field  Electron  Emission  Model 

The  analytical  and  statistical  model  used  to  predict  the  emission  current  from  a 
field  emission  tip  and  gate  electrode  structure  as  shown  in  Figure  1.5,  was  developed  by 
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Jensen69.  The  model  is  described  in  much  more  detail  in  another  article.69  The  emission 
model  takes  a  statistical  approach  to  estimating  the  current  emitted  by  field  emission 
cathodes  because  the  emitting  area  of  a  tip,  current  distribution  in  that  area,  number  of 
tips  contributing  to  the  measured  current,  and  current  per  tip  are  impossible  to  determine 
exactly  for  any  FEA  cathode.  Work  function,  tip  radius,  emission  cone  half-angle,  and 
cone  structure  half-angle  assume  a  range  of  values  in  a  field  emission  array  cathode. 
Because  each  tip  is  processed  in  parallel  and  exposed  to  the  same  environments,  the  cone 
structure  half-angles  and  work  functions  for  the  all  of  the  tips  in  the  array  should  be  very 
similar.  The  experimental  results  discussed  in  the  previous  chapter  showed  that  the 


emitting  area  of  the  tip  depends  on  the  cleanliness  of  the  tip  and  its  radius  of  curvature. 
Other  studies  have  shown  that  field  emission  from  a  cathode  can  originate  from 
protrusions  of  atomic  dimensions,  even  on  a  cone  structure  with  a  radius  of  curvature  of 
several  nanometers/6  Emission  from  a  cathode  microtip  is  not  uniform  within  the 
emitting  area  as  shown  in  the  previous  chapter.  The  emitting  area  of  tips  in  an  array  is 
then  best  represented  statistically  by  an  effective  tip  radius  and  distribution  of  tip  radii  in 
the  array  of  emitters.  These  parameters  can  be  determined  from  experimentally  acquired 
I-V  data  and  the  Jensen  FE  current  model.  It  has  been  successfully  used  to  model  the 
performance  of  Si  and  Mo  FE  cathodes.71 

Jensen’s  model  represents  the  tip  with  a  hyperbolic  shape;  the  field  (V/A)  at  the 

apex  of  the  tip  can  then  be  approximated  by 


Eqn.  4.3 


Ftip  ~  PgVg 


where  Vg  is  the  voltage  on  the  gate  electrode  with  respect  to  the  emitting  cone.  The 
electric  field  at  the  tip  surface  increases  with  decreasing  tip  radius,  rt,  and  half-cone 
angle,  |5C,  by  the  field  enhancement  factor,  Pg(l/A), 


Eqn.  4.4 


Pi 


— ,  and 

rt 
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Eqn.  4.5 


Jfc  =  — 
54 


r„ 

86+— 

r. 


cot  (Pc). 


tj 


With  these  relationships,  Fdp  (eV/A)  can  be  approximated  with  an  analytical  expression. 

The  current  density  is  approximated  by  the  Fowler-Nordheim  (F-N)  equation  for 

emission  current  density  Jpufe/fs-A2), 20,21 

Eqn.  4.6  JFN(F)  =  afilF2exp(-bfn/ F), 

where  the  Fowler-Nordheim  coefficients,  b&(eV/A)  and  ajd  eV2fs),  are  given  in  atomic 
units,  by 

Eqn.  4.7  =  O.64203/2and 

Eqn.  4.8  =  (1160)'1  exp(l4.3994fo/h  /  <t>2). 

When  the  current  density  is  integrated  over  the  surface  of  the  emitter,  the  tip  current  in 

microamperes  can  then  be  determined  by 

Eqn.  4.9  = 


The  area  factor,  b^fA2),  is  the  ratio  between  the  current  from  the  tip  and  the  current 


density  on-axis,  and  depends  upon  the  magnitude  of  the  gate  voltage  as 

/"  T-.  A 

Eqn.  4.10  barta  =  2  nr2  cos2  (ft ) 


F» 


ft,+sin2(ft)ft, 


tip  J 


The  current  emitted  from  an  array  of  tips  can  be  defined  by 


Eqn.  4.11 


where  2  is  the  distribution  factor  associated  with  a  spread  in  tip  radii,  As.  The  tip  radii 
are  distributed  according  to  r(s)=rt(l+s),  where  0<s<As.  For  a  single  tip  emitting,  As=0. 
However,  this  value  governs  the  distribution  of  tip  radii  with  N,jpS,  and  therefore  accounts 
for  the  non-uniformity  of  an  array.  The  field  for  a  tip  of  radius  r(s)=rt(l+s)  is  then 
F(s)=(l+cos)F(0),  where  F(0)=Ftip,  and 
Eqn.  4.12  c0  =  -^-(tan2(ft)  +  Pgr,f-1. 
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Experimental  data  are  used  to  estimate  rt  and  As.  If  I-V  data  are  presented  on  a 

Fowler-Nordheim  plot  (ln(I/V2)  vs.  1/V),  A^A/V2)  is  the  y-intercept  and  B^CV)  is  the 
slope  of  the  curve.  From  these  parameters  and  an  estimated  work  function,  effective  tip 
radius,  r„  and  distribution  parameter.  As,  can  be  determined.  Once  the  effective  tip 

radius  and  spread  in  tip  radius  throughout  the  array  of  emitters  are  known,  the  current 
from  the  array  can  be  calculated.  The  radius  distribution  parameter,  As,  represents  the 

spread  in  tip  radii.  In  an  array,  As=rma</r,  -1,  where  r^^.fl+As)  is  the  maximum  radius 
and  rt  is  the  minimum  radius.  During  exposure  to  xenon  environments,  the  performance 
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of  the  sharpest  tips  will  degrade  the  fastest  as  dr/dt  will  be  highest  for  those  tips. 

Sputtering  of  the  cathode  tips  will  decrease  the  spread  in  tip  radii  to  improve  the 
uniformity  of  the  array.  It  is  assumed  that  w  is  constant  during  the  exposure  and  that  As 

changes  in  time  with  Xe  exposure  as  rt  changes,  as  Eqn.  4.2  shows.  The  change  in 
effective  tip  radius  in  time  depends  on  *he  removal  rate  of  material  from  the  tip  emission 
area.  The  time  required  to  remove  a  inonolayer  of  material  from  the  tip  emission  area  is 
determined  using  the  tip  sputtering  model  discussed  in  the  next  section. 


4.1.2  Tip  Sputtering  Model 


Brodie  developed  the  field  emission  cathode  tip  sputtering  model  described  in  this 
section.72  It  is  described  in  more  detail  in  another  article.72  This  model  is  capable  of 
estimating  sputter  rates  of  microtips  by  ions  generated  by  the  emitted  electrons.  It  is  used 
here  to  estimate  the  change  in  tip  radius  in  time  when  operating  in  xenon  environments. 
Brodie  used  average  sputter  yield  and  ionization  cross-sections  in  his  calculations  for 
high  operating  voltages,  >100  V.  This  study  is  primarily  investigating  low  voltage 
operation  where  ion  and  electron  energies  are  low  and  ionization  cross-sections,  Q,  and 
sputter  yields,  Y,  are  exponentially  dependant  on  ion  and  electron  energies.  Accurate 
sputter  yield  and  ionization  cross  sections  are  very  important  to  these  calculations; 
therefore,  average  values  were  not  used  in  the  calculations  of  the  removal  rate  of  material 


from  the  emission  area  of  the  tips,  ns.  The  time,  te,  required  to  erode  a  single  layer  of 

atoms  from  the  emission  area  is  given  by 

2t&;2(1-cos<t) 

Eqn.  4.21  te  = - , 


where  d=0.416  nm  for  the  atomic  spacing  of  evaporated  molybdenum,  d=0.313  nm  for 

silicon.  The  removal  rate  of  material  from  the  tip  emitting  area  is  ns(atoms/s).  The 
emission  half-cone  angle  is  cr;  it  is  estimated  as  37°  in  this  model.  The  half-cone 

emission  angle  varies  significantly  for  tips  in  an  array.  Experimental  results  discussed  in 
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the  previous  chapter  showed  that  this  value  ranged  from  8°  to  44°  for  eight  carbide  tips 
with  an  average  value  of  28°.  Gomer  claims  that  the  current  density  drops  to 
approximately  10%  of  the  value  on  axis  by  45°  off  axis.19  A  half-cone  emission  angle  of 
37°  is  used  because  it  significantly  simplifies  the  mathematics  of  the  sputter  rate 
calculations  and  is  a  reasonable  estimate  to  represent  tip  half-cone  emission  angle  in  an 
array. 

The  removal  rate  of  material  from  the  tip  emitting  area  depends  on  the  flux  of 
ions,  energy  of  the  ions,  and  the  dependence  of  material  sputter  yields  on  ion  energy. 

The  number  of  ions  striking  the  emitting  area  per  second,  n,  is  given  by 

r  / 

Eqn.4.22  n  =  \ -^-NQiV^dr , 

J  e 

n 

and  the  number  of  atoms  sputtered  from  the  emitting  area  per  second,  ns,  is  given  by  . 

rm  j 

Eqn.  4.23  ns  =  f  -s-NQi{Vr)Y{Vr)dr , 

J  e 


where  the  field  emitted  electron  current  is  Ie  and  the  number  of  molecules  per  cm3  in  the 
chamber  is  N=  3.55xl016P  molecules/cm3  (P  is  the  pressure  in  Torr  at  room  temperature). 
The  cross-section  for  ionization  of  gas  molecules  by  electrons  rf  energy  eVr  is  Q(Vr). 
The  sputter  yield  in  atoms  per  ion  of  energy  eVris  Y(Vr). 

In  this  cathode  sputtering  model,  the  cathode  configuration  is  approximated  as  an 
isolated  sphere  of  radius  rt  at  ground  potential,  as  shown  in  Figure  4.1,  and  the  gate 
electrode  as  a  concentric  sphere  with  radius  rg.  The  radial  potential  distribution  on  the 
axis  of  the  tip  is  Vr,  and  can  be  approximated  by 


Eqn.  4.24 


£-f  i-i 


— -il 

\r.  j 


This  simple  concentric  spheres  model  gives  remarkably  good  results  compared  to  exact 
predictions  obtained  by  computer  solutions  of  Laplace’s  equation  and  geometric 
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boundary  conditions  because  close  to  the  emitting  area  of  the  tip,  the  field  lines  are  radial 
and  fall  off  inversely  as  the  square  of  the  distance.73  The  Jensen  model  more  accurately 
predicts  the  electric  field  at  the  tip  surface,  which  is  more  important  when  calculating  the 
field  emission  current  from  a  tip.  However,  the  Brodie  model  is  sufficient  for  predictions 
of  the  potential  distribution  near  the  tip  to  determine  ion  production  near  the  cathode.  For 
this  analysis,  the  electric  field  and  current  over  the  emitting  area  of  the  tip  are  assumed  to 
be  uniform. 


rt 


actual  model 


Figure  4.1  The  actual  cathode  configuration,  and  the  configuration  used  in  the 
microtip  sputtering  model  to  predict  the  potential  field  near  the  microtip  and 
electric  field  at  the  microtip. 


The  ion  flux  to  the  tip  depends  on  the  local  gas  temperature  and  the  positions  at 
which  ionization  occurs.  If  an  ion  is  formed  with  initial  cross  radial  velocity,  v„ 


corresponding  to  a  temperature  of  Vt(eV),  then  only  the  ions  formed  within  a  radial 
distance  rm,  will  strike  the  emitting  area  of  the  tip,  where 


Eqn.  4.25 


rm  =  rt 


f  v  V 

s 

UJ ' 


Usually  the  ion  temperature  is  estimated  to  be  0.03  eY,  which  assumes  that  the  ions  are  in 
equilibrium  with  their  surrounding  vacuum  chamber  walls  with  the  same  temperature  as 
the  neutral  particles.  If  a  higher  temperature  of  ions  is  assumed,  the  flux  of  ions  to  the 
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emitting  area  of  the  tips  will  decrease.  Therefore  the  assumption  that  Vt=0.03  eV  will 
provide  an  upper  limit  on  dr/dt. 

The  significant  radial  positions  used  in  this  model  are  shown  in  Figure  4.2.  The 
radial  position  at  which  the  energy  of  the  emitted  electrons  reaches  the  threshold  of 
ionization  (Vr(r,)  =  Vi  =  ionization  threshold)  is  r*.  The  radial  position  from  which  an  ion, 
if  formed,  gathers  sufficient  energy  to  reach  the  threshold  for  sputtering  when  it  strikes 
the  tip  surface  (Vr(rs)  =  Vs  =  sputtering  threshold)  is  rs  The  radial  position  beyond  which 
the  average  ion  formed  will  not  strike  the  emitting  area  of  the  tip  is  rm. 


Figure  4.2  Significant  radial  positions  used  in  the  model  of  microtip  sputtering. 

The  calculations  made  by  Brodie  employed  average  values  of  Vr  and  then  used 
corresponding  Q  and  Y  values.  In  this  study,  data  curves  and  interpolation  were  used  to 
estimate  these  parameters  because  these  parameters  are  exponentially  sensitive  to  ion  and 
electron  energies  in  the  energy  range  relevant  in  this  investigation.  The  integral  for  ns 
was  then  evaluated  numerically  for  hundreds  of  points.  Ionization  cross-sections  from 
Brown74  used  in  the  calculations  are  shown  in  Figure  4.3.  Interpolation  was  used  to 
determine  values  between  data  points.  Sputter  rates  for  xenon  ions  on  molybdenum  and 
silicon  targets  were  used  in  the  model  because  the  thruster  environment  is  primarily 
xenon.  The  relationship  between  sputter  yield  and  ion  energy  is  discussed  in  the  next 
section. 
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Engery  (eV) 

Figure  4.3  Xe  ionization  cross-sections,  Q„  from  Brown.74 


4.1.3  Sputter  Yield  Model 

Sputter  yield  and  energy  threshold  values  used  in  this  study  for  xenon  ions 
bombarding  molybdenum  and  silicon  targets  with  a  zero  degree  angle  of  incidence  were 
calculated  using  a  model  developed  by  Yamamura  et  al.67  using  fit  parameters, Q,  from 
Matsunami  at  al.75  The  electric  field  within  rm  of  the  tip  is  primarily  radial.  Therefore  the 
ion  incidence  angle  can  be  approximated  by  zero.  Because  experimental  data  exist  only 
above  100  eV  ion  energies,  the  accuracy  of  the  model  is  unknown  at  energies  below  this 
level.  The  energy  of  ions  bombarding  the  emitting  area  of  the  cathodes  in  the 
configurations  of  interest  will  primarily  be  below  100  eV. 

In  the  model  developed  by  Yamamura  et  al.67,  the  sputter  yield  is  defined  as 

if- 8 
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and  the  energy  threshold  for  sputtering  is  defined  by, 


Eqn.  4.26 
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Eqn.  4.27 
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where  the  subscripts  I  and  II  refer  to  the  incident  and  target  particles,  respectively.  Q  is 
an  empirically  derived  parameter  for  a  specific  ion  target  combination,  Us(eV)  is  the 
sublimation  energy  of  the  target  materia1,  and  E(eV)  is  the  ion  energy.  Q  values  are 
presented  for  different  ion  target  combinations  by  Matsunami  at  al.75  Other  parameters 
used  in  the  model  are 
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The  parameters  used  in  this  investigation  for  molybdenum  and  silicon  targets  bombarded 
by  xenon  ions  are  shown  in  Table  4.1. 
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Table  4.1  Ion  and  target  material  parameters  used  in  the  sputter  yield  model. 


Xe 

Mo/MoXe 

Si/SiXe 

M 

131.3 

95.94 

28.1 

Z 

54 

42 

14 

Us(eV) 

6.82 

4.63 

Q 

0.84 

0.78 

Eh.  (eV) 

49 

91 

Table  4.1  shows  E^  values  predicted  by  the  model  for  Mo  and  Si  targets  being 
bombarded  by  xenon  ions.  Other  models  predict  even  higher  values  for  Ej®  while  some 
experimental  measurements  extrapolated  to  low  energies  yield  E*  estimates  for  a  Mo 
target  and  Xe  ion  combination  that  is  27  eV.77  Experimental  results  presented  in  the 
previous  chapter  show  that  cathodes  were  damaged  by  ion  bombardment  when  maximum 
ion  energies  were  lower  than  the  E*  values  shown  in  Table  4.1.  Also,  these  values  for  E^, 
led  to  great  underestimates  of  dl/dt  using  the  FEA  cathode  performance  and  sputtering 
models.  To  obtain  better  correlation  between  measured  and  theoretical  dl/dt,  lower  E^  . 
values  were  used  in  the  models.  These  Em  values  were  more  consistent  with  experimental 
results  presented  in  the  previous  chapter  and  generated  sputter  yield  curves  that  provided 
a  better  fit  through  measured  sputter  yield  values  at  ion  energies  less  than  300  eV.  This 
section  discusses  the  modeling  results  and  comparisons  to  experimental  results. 

The  model  developed  by  Yamamura  et  al.  was  used  to  predict  that  the  E^  values 
for  xenon  ions  sputtering  silicon  and  molybdenum  targets  are  91  eV  and  49  eV, 
respectively.  This  model  was  also  used  to  determine  the  relationship  between  ion  energy 
and  sputter  yield;  these  relationships  are  shown  in  Figure  4.4  and  Figure  4.5  for  Mo  and 
Si,  respectively.  The  cathode  sputtering  model  was  used  to  estimate  the  maximum 
energy  of  the  ions  bombarding  the  emitting  area  of  the  tips,  Eimat(rm),  at  any  operating 
condition.  Results  from  experiments  during  which  cathodes  were  exposed  to  xenon 
while  operating  showed  that  MCNC  Si  cathodes  were  damaged  when  gate  electrode 
voltages  were  at  85, 70,  and  60  V  (Eimax  =  -55  eV).  The  cathode  was  not  damaged  during 
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1  hr.  in  2x1  O'5  Torr  of  xenon  with  50  V  applied  to  the  gate  electrode  (E^  =  ~46  eV). 
These  results  are  preliminary  because  the  experiments  lasted  only  a  few  hours  and  some 
of  them  have  not  yet  been  repeated;  however,  they  suggest  that  the  E*  for  sputtering  of  Si 
by  xenon  ions  is  below  55  eV.  When  the  experimental  results  are  combined  with  results 
of  the  performance  model  discussed  in  this  chapter,  a  better  estimate  of  the  E^,  can  be 
obtained. 


Figure  4.4  Sputter  yields  for  a  Mo  target  being  bombarded  by  Xe  ions.  Shown  on 
the  graph  are  theoretical  curves  for  two  values  of  Ett,  and  experimental  values 
measured  by  Rosenberg  and  Wehner.7® 
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Figure  4.5  Sputter  yields  for  a  Si  target  being  bombarded  by  Xe  ions.  Shown  on  the 
graph  are  theoretical  curves  for  three  values  of  E^,  and  experimental  values 
measured  by  Rosenberg  and  Wehner.78 

The  results  of  molybdenum  cathode  testing  are  reported  with  more  confidence 
because  experiments  were  conducted  at  low  voltages  lasted  for  several  hours  and  were 
repeated.  The  model  developed  by  Yamamura  et  al.  predicted  that  £^=49  eV.  Cathode 
experiments  in  2xl0'5  Torr  of  xenon  showed  that  performance  changes  sometimes 
occurred  with  Vg  at  50  V  where  E,m„  is  46  eV,  therefore  the  must  be  below  49  eV . 
Experimental  results  combined  with  results  of  the  performance  model  discussed  in  this 
chapter,  is  required  for  a  better  estimate  of  the  E^. 

4.2  Modeling  Results 

The  performance  degradation  model  discussed  in  this  chapter  requires  input 
including  experimental  I-V  data  from  a  FEA  cathode,  cathode  geometry  parameters,  and 
sputter  yield  parameters  before  FEA  cathode  performance  degradation  can  be  predicted. 

It  incorporates  the  physics  of  field  electron  emission,  ionization  of  xenon,  and  sputtering 
of  the  tip  emitting  area  by  ion  bombardment  and  are  used  to  estimate  the  change  in  radius 
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of  curvature  of  cathode  microtips  and  corresponding  current  expected  during  electron 
emission  into  xenon  environments.  This  model  can  only  predict  changes  in  the  cone 
structure  by  increases  in  the  effective  tip  radius,  rt.  If  tip  sharpening  occurs  or  a  reduction 
in  cone  half-angle,  cathode  performance  will  improve.  It  is  possible  that  this  process  did 
occur  during  Xc  exposure  experiments  with  Va>80  V  and  is  responsible  for  the 
performance  improvements  discussed  in  the  previous  chapter.  Therefore,  the 
performance  degradation  model  is  expected  to  predict  cathode  performance  degradation 
in  Xe  environments  best  when  Va  is  low  and  damage  from  ion  bombardment  primarily 
occurs  to  the  emitting  area  of  the  tips.  This  situation  is  of  primary  interest  for  the  electric 
propulsion  application.  In  the  xenon  environment  of  a  Hall  thruster,  the  virtual  anode 
potential  will  only  be  ~20  V  above  the  gate  electrode  potential.  Therefore,  operating  with 
Va  at  20V  +Vg  best  represents  the  cathode  and  thruster  operating  configuration. 

4.2.1  Mo  SRI  International  Cathode 

The  performance  degradation  of  a  Mo  cathode  operating  in  a  Xe  environment 
with  a  low  Va  (80  V)  and  Vg  at  65  V  was  modeled  with  remarkable  results.  The 
measured  response  of  the  cathode  current  to  increases  in  xenon  pressure  is  shown  in 

Figure  4.6  as  I^.  I-V  data  taken  before  the  exposures  are  shown  in  Figure  4.7  as  1^. 
The  Jensen  performance  model  was  used  with  these  data  to  estimate  As  and  rt  shown  in 

Table  4.2  as  data  set  I.  Figure  4.7  also  shows  the  I-V  relationship  calculated  with  the 
Jensen  electron  emission  model  using  the  experimentally  determined  As  and  r„  as  IcaJc. 

The  cathode  performance  decay  caused  by  exposure  to  the  xenon  environment  was 
predicted  using  the  model.  If  an  E*  of  49  eV  was  assumed,  as  the  model  developed  by 
Yamamura  et  al.  predicted,  the  performance  decay  rate  was  greatly  underestimated. 

Using  Et,,  at  39  eV,  the  cathode  performance  degradation  model  predicted  the 
performance  decay  shown  in  Figure  4.6  as  1,^.  The  correlation  between  experiment  and 
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theory  is  remarkable.  The  sputter  yield  curve  generated  by  the  model  model  developed 
by  Yamamura  et  al.67  using  E*  at  39  eV  provides  a  better  fit  to  values  measured  by 
Rosenberg  and  Wehner78  at  ion  energies  below  300  eV,  as  shown  in  Figure  4.4.  Cathode 
parameters  obtained  from  I-V  traces  taken  before  and  after  the  Xe  exposure  are  shown  in 
Table  4.2.  Changes  in  cathode  parameters  predicted  by  the  performance  degradation 

model  are  shown  as  data  set  m  in  Table  4.2.  During  exposure  to  xenon,  rt  increased  and 
As  decreased;  the  sharpest  tips  were  damaged  at  the  fastest  rate  increasing  uniformity  in 

tip  radius  across  the  array,  and  decreasing  As. 


eV. 


82 


40 


Figure  4.7  Mo  cathode  I-V  data  measured  experimentally  before  the  Xe  exposure 
and  predicted  by  Jensen’s  performance  model. 


Table  4.2  Cathode  parameters  before  (I)  and  after  (II)  Xe  exposure  obtained  from 
I-V  data  and  predicted  (III)  by  the  performance  degradation  model. 


4>w(eV) 

Pc  (rad.) 

rg(A) 

Ntips 

Aro(A/V2) 

BpnCV) 

r,(A) 

As 

I 

4.35 

0.26 

0.0013 

819 

43.6 

50.0 

4.35 

0.26 

0.0024 

952 

53.5 

42.5 

4.35 

53.6 

40.8 

The  performance  degradation  model  was  also  used  to  study  the  effect  of  work 
function  on  E*  estimates.  The  work  function  typically  used  for  Mo  is  4.3-4.4  eV.79 
Performance  degradation  rates  were  predicted  assuming  three  different  work  function 

values  and  the  same  E*  value  of  39  eV.  The  results  showed  that  performance  degradation 
rates  predicted  using  each  <J)W  value  were  consistent  with  experimental  results  when  an  E* 

value  of  39  eV  was  used.  Therefore,  the  E*  estimates  do  not  seem  to  be  significantly 
affected  by  the  assumed  work  function  in  the  range  of  possible  values. 

The  performance  degradation  model  was  applied  to  predict  the  gate  electrode 
voltage  threshold  for  significant  performance  degradation  for  comparison  to  experimental 


83 


observations.  Experimental  results  presented  in  the  previous  chapter  in  Figure  3.19  show 
that  at  58  V  the  SRI  cathode  current  dropped  by  6%  during  10  min.  with  Va  at  100  V. 

The  cathode  parameters  before  the  exposure  are  shown  in  Table  4.3  as  data  set  I.  Figure 
4.8  shows  that  the  model  predicted  the  dl/dt  measured  with  Vg  at  64  V  for  this  cathode. 

A  comparison  of  cathode  parameters  before  and  after  the  exposure  to  xenon  in  Table  4.3 
shows  that  As  increased  and  r,  barely  changed.  Ion  bombardment  of  the  microtips  could 

have  changed  half-cone  angle  or  nanoprotrusion  emission  sites  while  sputtering  the  tip 
emitting  area.  The  same  cathode  was  used  in  an  experiment  which  showed  that  the 
cathode  performance  is  affected  by  the  voltage  applied  to  the  anode  in  xenon 
environments.  With  Vg  at  50  V,  the  cathode  performance  increased  as  Va  was  increased 
from  60  V  up  to  120  V.  With  increasing  anode  voltage,  more  ions  could  be  bombarding 
the  microtips  off-axis  increasing  the  sputtering  of  the  cone  sidewalls,  which  sharpens  the 
tips  and  reduces  the  cone  half-angle.  With  Va  at  100  V  there  is  a  combination  of  changes 
occurring  to  the  cathode  to  improve  and  degrade  cathode  performance.  At  lower  Va, 
modeling  results  are  more  consistent  with  experimental  results.  Longer  duration 
experiments  are  required  at  lower  anode  voltages  to  better  understand  discrepancies 
between  theoretical  and  experimental  results. 

Initial  cathode  performance  requirements  are  less  than  10%  current  reduction 
during  100  hours  in  2xl0'5  Torr  of  xenon.  The  performance  degradation  model  was  used 
to  investigate  the  effect  of  xenon  exposure  on  cathode  performance  during  100  hours. 

The  cathode  parameters  shown  in  Table  4.3,  data  set  I,  were  used  in  the  model.  The 
performance  degradation  model  predicted  that  if  this  cathode  is  operated  in  2xl0'5  Torr,  it 
will  demonstrate  a  current  reduction  of  5%  when  Vg  is  55  V  after  10  hours  of  operation. 
Performance  degradation  predicted  by  the  model  for  the  same  cathode  in  the  same 
environment  at  various  gate  voltages  during  100  hours  is  shown  in  Figure  4.9.  During 
100  hours  of  operation,  the  current  will  decrease  by  34  %  of  its  original  value  with  Vg  at 
55  V.  It  was  determined  that  the  gate  voltage  for  this  cathode  cannot  be  higher  than  52  V 
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to  limit  the  current  decay  to  less  than  10%  during  100  hours.  The  current  for  this  cathode 
at  52  V  is  only  10.4  nA.  At  this  voltage,  107  of  these  cathodes  would  be  required,  or 
5x10"  tips  to  emit  100  mA.  The  size  of  this  cathode  would  be  8  m2.  The  performance 
degradation  of  the  cathode  only  considered  the  ions  created  locally.  If  additional  ion 
populations  are  present,  operating  voltages  will  be  limited  to  even  lower  values  and 
cathode  dimensions  will  be  significantly  larger. 


Table  4.3  Cathode  parameters  before  (I)  and  after  (II)  xenon  exposure  which  are 
used  in  the  performance  decay  model  for  Va  at  100  V  in  2xl0'5  Torr  of  xenon. 


<MeV) 

Pc  (rad.) 

rg(A) 

Apn(A/V2) 

Bfn(V) 

r.(A) 

As 

I  4.35 

0.26 

0.0043 

1092 

64.9 

18.2 

H  4.35 

0.26 

0.00246 

1104 

66 

32 

Figure  4.8  Performance  decay  rate  predicted  by  the  performance  degradation 
model  for  a  Mo  cathode  operating  in  a  Xe  environment  with  a  pressure  of  2x10  s 
Torr  for  a  cathode  with  parameters  (I)  shown  in  Table  4.3. 
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Figure  4.9  Field  emission  array  cathode  performance  decay  predicted  by  the 
performance  degradation  model  for  a  Mo  cathode  operating  in  2x10  s  Torr  of  Xe  for 
100  hours  at  various  gate  electrode  voltages  with  initial  cathode  parameters  (I) 
shown  in  Table  4.3. 


In  the  cathode  environment  of  a  Hall  thruster  with  a  virtual  anode  potential  only 
20  V  greater  than  the  gate  electrode  potential,  ions  bombarding  the  cathode  that  originate 
near  the  tips  will  primarily  affect  changes  in  rt.  In  the  Hall  thruster  environment,  the 
CEX  ion  population  will  contribute  an  additional  ion  current  density  of  0.02  mA/cm2  to 
the  cathode  which  is  distributed  uniformly  to  the  surface  of  the  cathode  cones,  as 
discussed  in  Chapter  2.  The  significant  effect  of  the  CEX  population  on  performance 
degradation  of  the  Mo  cathode  with  parameters  (I)  shown  in  Table  4.2,  is  also  shown  in 
Figure  4.6.  The  cathode  performance  degradation  attributable  to  local  ions  is  negligible 
in  comparison.  Not  only  will  this  ion  population  degrade  cathode  performance  by 
sputtering  the  microtip  radius  of  curvature  to  degrade  performance  as  shown  in  Figure 
4.6,  it  will  also  sputter  material  off  of  the  cone  sidewalls  which  will  redeposit  on  the 
insulator  walls.  With  ion  energies  at  eVg+20eV=85eV,  and  xenon  pressure  at  2xl0'5 
Torr,  the  gate  could  be  shorted  to  the  tips  within  several  minutes.  For  thruster  and  Mo 
FEA  cathode  technologies  to  be  compatible,  either  the  potential  between  the  plasma  and 
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microtips  must  be  below  the  voltage  threshold  for  sputtering,  or  the  cathode  must  be 
shielded  from  this  population  of  ions.  To  shield  the  microtips  from  CEX  ions,  CLAIR 
was  developed.  CLAIR  is  a  cathode  lens  and  ion  repeller  (CLAIR)  configuration 
proposed  to  prevent  the  CEX  ions  from  bombarding  the  cathode  microtips.  CLAIR  has  a 
series  of  electrodes  similar  to  the  gate  electrode  that  electrostatically  retards  CEX  ions 
and  focuses  the  electrons  through  the  highest  voltage  electrode.  The  CLAIR 
configuration  is  shown  and  discussed  in  Appendix  B.  If  the  cathode  cannot  be  shielded 
from  these  ions,  then  molybdenum  cathodes  must  be  limited  to  operating  voltages  below 
-19  V  to  achieve  lifetimes  greater  than  100  hours. 

4.2.2  Si  MCNC  FE  Cathodes 

Silicon  FEA  cathode  performance  degradation  was  also  modeled,  and  the  results 
were  compared  to  experimental  results.  Results  of  modeling  and  experimental 
measurements  of  the  performance  degradation  of  Si  cathodes  in  Xe  environments  showed 
that  the  E*  is  also  much  lower  than  originally  predicted  by  the  model  developed  by 
Yamamura  et  al.  Figure  4.10  shows  the  measured  cathode  response  to  an  increase  in 
xenon  pressure  with  Vg  at  70  V  and  Va  at  70  V.  Cathode  performance  was  poor  becaus- 
this  cathode  had  been  damaged  in  several  other  experiments.  The  cathode  current  before 
the  exposure  with  Vg  at  70  V  was  3.2  nA  with  good  emission  stability.  The  cathode 
current  dropped  to  1.6  nA  during  1  hour  at  2xl0'5  Torn  Correlation  between  experiment 
and  theory  is  good  if  an  E^  value  of  48  eV  is  used  in  the  sputter  yield  model.  Sputter 
yield  curves  with  different  E^  values  are  shown  in  Figure  4.5.  With  E*  at  48  eV  a  better 
fit  to  experimental  measurements  of  sputter  yield  data  obtained  by  Rosenberg  and 
Wehner  is  obtained  than  with  E*  at  91  eV.  Performance  degradation  curves  were 
determined  using  a  range  of  work  function  values,  as  shown  in  Figure  4.10.  Results 
showed  that  an  E*  value  at  48  eV  produced  theoretical  results  consistent  with 
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experimental  results  for  a  range  of  possible  silicon  work  functions.  The  work  function  of 
silicon  depends  on  how  it  has  been  doped.  A  work  function  of  4.05  eV  is  assumed  for  the 
MCNC  Si  cathodes  tested.80  To  avoid  sputter  damaging  the  MCNC  Si  cathode  by  ions 
created  locally,  operating  voltages  must  be  limited  to  values  below  53  V.  At  this  voltage 
is  47  eV,  below  the  E*  for  Xe  ions  to  sputter  Si  targets. 
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Figure  4.10  Measured  and  predicted  Si  cathode  current  response  to  an  increase  in 
Xe  pressure  up  to  2xl0'5  Torr  with  Vg  and  Va  at  70  V. 


Modeling  results  were  compared  to  another  set  of  experimental  data  acquired 
during  Si  cathode  operation  in  a  xenon  environment.  This  experiment  was  conducted 
with  a  much  higher  voltage  applied  to  the  anode,  400  V,  and  higher  emission  current.  In 
experiments  conducted  with  Va>80  V,  it  is  believed  that  the  performance  degradation 
model  is  not  sufficient  to  account  for  all  of  the  changes  in  cathode  structure  caused  by 
xenon  ion  bombardment.  Figure  4.1 1  shows  the  response  of  a  Si  cathode  from  MCNC  to 
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an  increase  in  Xe  pressure.  The  cathode  parameters  before  the  exposure  are  shown  in 
Table  4.4.  Figure  4. 1 1  shows  the  performance  decay  predicted  by  the  model  using  the 
initial  cathode  parameters.  The  model  slightly  over  predicts  the  performance  degradation 
rate.  It  is  possible  that  this  result  occurs  because  of  the  high  voltage  applied  to  the  anode, 
400  V.  At  thi*:  voltage,  there  could  be  a  combination  of  tip  blunting  to  degrade 
performance  and  tip  sharpening  and  half  cone  angle  reduction  to  improve  cathode 
performance  as  shown  in  Figure  3.22.  Reductions  in  work  function  could  also  be 
occurring  as  the  xenon  ions  sputter  clean  the  cathode  surface.  The  increase  in  current  as 
the  xenon  pressure  is  increased  initially  could  be  attributed  to  decreases  in  cathode  work 
function.  Silicon  emitters  naturally  have  a  Si02  film  on  the  tips  that  increases  the  work 
function.  If  this  film  is  sputtered  off,  work  function  decreases  would  be  expected.  Any 
change  in  cathode  configuration  to  improve  performance  will  result  in  a  slower 
performance  degradation  rate  and  discrepancies  between  experimental  and  theoretical 
results. 


Figure  4.11  Measured  and  calculated  Si  cathode  performance  degradation  as  the 
Xe  pressure  increases  up  to  7X10"6  Torr  with  Vg  at  86  V  and  Va  at  400  V. 
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Table  4.4  MCNC  Si  cathode  parameters  before  the  exposure  to  xenon  during 
operation. 


4>w(eV) 

Pc  (rad.) 

rg(A)  Ntips  rt(A)  As 

4.0 

0.26 

10000  16000  45.2  185 

Experimental  and  performance  decay  modeling  results  were  used  together  to  gain 
a  better  understanding  of  how  Mo  and  Si  FEA  cathodes  will  respond  to  the  xenon 
environment  of  a  Hall  thruster.  Experimental  results  showed  that  the  E^  values  for  xenon 
ions  sputtering  Mo  and  Si  are  much  lower  than  the  model  developed  by  Yamamura  et 
al.67  predicted.  The  modeling  and  experimental  results  showed  that  the  gate  voltage 
should  be  below  43  V  and  53  V  to  avoid  sputter  damaging  the  emitting  area  of  the  Mo 
and  Si  tips,  respectively,  by  ions  created  locally.  It  was  shown  that  the  population  of 
CEX  ions  created  near  the  thruster  will  be  even  more  damaging  to  the  cathode  than  the 
ions  created  locally.  To  prevent  local  and  CEX  ions  from  sputter  damaging  the  Mo  and 
Si  tips  and  CEX  ions  from  shorting  the  gate  to  the  tips,  the  gate  voltage  should  be  lower 
than  19  V  and  28  V,  respectively. 

There  are  a  lot  of  assumptions  made  to  predict  electron  emission  and  tip 
sputtering  by  ion  bombardment  that  greatly  simplify  modeling  the  physics  of  these 
processes.  In  spite  of  such  simplifications,  the  models  are  able  to  capture  the  majority  of 
the  physics  so  that  modeling  results  show  impressive  correlation  with  experimental 
results  in  the  cathode  and  test  configurations  studied.  Experimental  and  theoretical 
results  also  lead  to  a  better  understanding  of  how  and  when  to  use  the  model.  The  results 
of  this  investigation  can  now  be  applied  in  parametric  studies  of  cathode  design  and 
performance  in  Xe  environments. 
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4.3  Application  of  Modeling  and  Experimental  Results 


The  cathodes  used  in  this  investigation  demonstrated  less  than  ideal  FE  cathode 
performance  because  they  did  not  represent  state-of-the-art  FE  cathode  technology  and 
the  anode  and  cathode  could  not  be  cleaned  in  the  vacuum  facility  prior  to  testing.  Any 
FEA  cathodes  are  very  difficult  to  acquire;  state-of-the-art  FE  A  cathodes  were  impossible 
to  acquire.  Cathodes  with  better  uniformity  and  smaller  dimensions  have  demonstrated 
better  performance  in  other  UHV  test  facilities.  Even  the  Mo  cathode  used  in  this 
investigation  demonstrated  much  better  performance  at  SRI  International  before  it  was 
tested  at  JPL.  However,  the  experiments  produced  the  results  required  in  this  study.  This 
model  can  now  be  applied  to  predict  the  performance  degradation  of  any  cathode  once 

the  cathode  geometry  is  known  and  an  I-V  trace  has  been  obtained  in  UHV  to  estimate  rt 
and  As.  The  initial  cathode  performance  objective  is  emission  of  100  mA  from  ~1  cm2  in 

2x1  O'5  Torr  of  xenon  for  100  hours  with  less  than  10  %  current  reduction  from  local  and 
CEX  ion  bombardment.  The  performance  of  the  cathodes  evaluated  at  JPL  indicated  that 
the  cathodes  used  in  this  study  could  not  meet  these  requirements,  however,  cathodes 
tested  elsewhere  may  satisfy  the  requirements.  This  section  presents  the  results  of  an 
investigation  into  the  theoretical  performance  of  more  ideal  c^’iode  configurations  in 
Hall  thruster  environments  to  determine  the  cathode  configurations  required  to  be 
compatible  with  EP  systems. 

The  performance  of  a  Mo  cathode  configuration  from  SRI  International  was 

modeled  in  a  Hall  thruster  environment.  This  cathode,  #52M,  was  tested  at  SRI 
International  where  an  I-V  curve  was  obtained  to  determine  As  and  rt.  The  actual  cathode 

parameters  are  shown  in  Table  4.5  in  the  first  row  and  first  5  columns  in  bold  type.  The 
predicted  performance  of  cathodes  with  similar  configurations  is  also  presented  in  that 
table  to  show  how  varying  some  parameters  would  affect  the  cathode  performance. 

Upper  limits  on  operating  voltages  were  determined  for  operating  conditions  with  and 
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without  a  CEX  ion  population  so  that  the  current  did  not  drop  by  more  than  10%  during 
100  hours.  The  number  of  tips  in  an  array  was  scaled  to  attain  approximately  100  mA. 
Table  4.5  shows  the  number  of  tips  required  to  emit  -100  mA  from  different  cathode 
configurations  with  lifetimes  greater  than  100  hours.  Higher  gate  voltage  operation,  -43 
V,  is  tolerable  when  the  CEX  ion  population  is  not  contributing  to  cathode  microtip 
sputtering.  With  the  original  cathode  configuration  (case  I)  and  a  gate  electrode  at  43  V, 
1.3xl08  tips  are  required.  The  packing  density  of  this  cathode,  #52M,  is  16,000,000 
tips/cm2,  therefore,  this  cathode  would  be  8  cm2  in  size.  Much  smaller  cathodes  are 
required  for  the  compatibility  of  Hall  thrusters  and  FEA  cathodes.  The  size  of  this 
cathode  would  be  acceptable  if  the  packing  density  is  increased  by  a  factor  of  10.  The 
number  of  tips  required  from  similar  cathode  configurations  is  also  shown  in  Table  4.5. 
Reducing  the  gate  aperture  radius  from  2000 A  to  1000 A  reduces  the  number  of  tips 
required  in  the  cathode  by  an  order  of  magnitude.  The  cathode  described  by  case  2 
should  be  compatible  with  a  Hall  thruster  without  CEX  ions  bombarding  the  microtips. 
Decreasing  the  tip  radius  and  increasing  the  uniformity  of  tip  radii  in  an  array  can  also 
significantly  improve  the  performance  of  the  cathode  so  that  the  cathode  size  is  less  than 
1  cm2.  Increasing  the  packing  density  in  an  array  by  an  order  of  magnitude  will  also 
reduce  the  size  of  the  cathode  required  to  <1  cm2. 

Presence  of  a  CEX  ion  population  significantly  reduces  the  upper  limit  on  the 
operating  voltage.  With  the  CEX  ions,  voltage  of  the  gate  electrode  must  not  be  higher 
than  20  V  for  a  Mo  cathode.  At  this  low  voltage,  a  prohibitively  large  number  of 
microtips,  4.5xl0ls-3.5xl0n,  are  required  to  attain  the  desired  current,  100  mA. 
Improvements  to  the  cathode  uniformity  of  microtip  radius  could  significantly  reduce  the 
number  of  microtips  required  for  the  desired  current.  This  cathode  would  still  be 
prohibitively  large  because  it  would  be  considerably  larger  than  the  propulsion  system. 
While  this  cathode  demonstrated  excellent  performance  in  UHV  environments,  it  could 
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only  be  used  with  a  Hall  thruster  if  a  filter  like  CLAIR  is  used  to  prevent  CEX  ions  from 
bombarding  the  cathode  tips  so  that  operation  ~42  V  is  an  option. 


Table  4.5  SRI  International  #52M  Mo  cathode  parameters  (in  bold)  and  predicted 
performance  using  different  parameters,  with  and  without  a  CEX  ion  bombardment 
during  100  hours. 


case 

<1>w 

(eV) 

Pc 

(rad.) 

(l) 

r, 

(A) 

As 

vg 

(V) 

•^tips 

Io 

(mA) 

If 

(mA) 

1-Mo 

4.41 

0.1974 

2000 

40.4 

1.98 

43 

1.3xl08 

108 

107 

2-Mo 

4.41 

0.1974 

1000 

40.4 

1.98 

42 

1.4xl07 

105 

101 

3-Mo 

4.41 

0.1974 

1000 

30 

1.98 

42 

1.2xl06 

108 

107 

4-Mo 

4.41 

0.1974 

2000 

40.4 

1 

43 

6.5xl07 

107 

106 

5-Mo 

4.41 

0.1974 

2000 

30 

1.98 

43 

7xl06 

100 

95 

6-Mocex 

4.41 

0.1974 

2000 

40.4 

1.98 

20 

4.5xl015 

108 

101 

7-Mocex 

4.41 

0.1974 

1000 

40.4 

1.98 

20 

4xl013 

104 

98 

8-MQcex 

4.41 

0.1974 

1000 

30 

1.98 

20 

3x10" 

107 

101 

Molybdenum  cathodes  fabricated  at  MIT/LL  have  demonstrated  record  low 
operating  voltages  because  of  incredibly  small  cathode  dimensions.  FEA  cathodes  at 
MIT/LL  have  been  fabricated  with  rg  only  800  A  and  packing  densities  of  109/cm2.  A 
current  of  1  uA  has  been  demonstrated  at  25  V  with  a  molybdenum  cathode  that  has  900 
tips  (2.8  nA/tip)  and  0.32  pm  tip-to-tip  spacing.  A  cathode  with  the  same  configuration 

and  quality  could  supply  100  mA  at  25  V  from  a  3x3  mm  cathode  with  9xl07  tips  and 
possibly  100  mA  from  less  than  1  cm2  at  20  V.  Scaling  up  the  cathode  in  number  of  tips 
with  the  same  performance  and  uniformity  is  the  challenge  yet  to  be  met  so  that  Mo 
cathodes  are  compatible  with  EP  systems.  Mo  cathodes  may  only  be  compatible  with  EP 
systems  if  CLAIR  is  used  to  prevent  CEX  ions  from  bombarding  the  microtips. 

FEA  cathodes  employing  other  materials,  which  are  currently  being  developed, 
could  demonstrate  performance  superior  to  the  Mo  cathodes.  A  Si  cathode  with 
dimensions  much  smaller  than  the  cathodes  used  in  the  experiments  discussed  in  Chapter 
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3  could  run  at  higher  gate  voltages  and  require  fewer  microtips  for  the  lifetime  and 
current  requirements,  as  shown  in  Table  4.6.  The  performance  of  Si  cathodes  could  be 
better  than  Mo  for  this  application  because  the  work  function  is  lower  and  the  for 
sputtering  is  higher.  Even  considering  the  CEX  ion  population,  silicon  cathodes  with  the 

dimensions  phuwn  in  cases  12  and  13  should  be  able  to  operate  at  100  mA  for  100  hours 
if  tips  are  fabricated  on  0.5  pm  centers  or  smaller  so  that  the  cathode  dimensions  will  be 

less  than  1  cm2.  A  current  of  2.5  nA  has  been  measured  from  a  single  silicon  microtip 
structure  at  25  V.68  If  silicon  FEA  cathodes  could  be  fabricated  with  great  uniformity  in 
tip  radius  across  an  array,  they  could  provide  100  mA  of  current  from  an  array  with  less 
than  106  microtips  with  an  emitting  area  <1  cm2.  Fabricating  FEA  cathode  with  the 
uniformity  required  is  currently  the  greatest  challenge  in  FEA  cathode  fabrication.  Other 
cathode  materials  can  further  improve  cathode  performance.  HfC  and  ZrC  have 
demonstrated  work  function  values  between  3.3  and  3.5  eV.81,82  Table  4.6  shows  the 
reduction  in  number  of  tips  required  if  the  work  function  is  reduced  to  3.3  eV  for  a 
particular  cathode  configuration.  The  number  of  tips  required  was  reduced  by  almost  two 
orders  of  magnitude  by  reducing  the  work  function  from  4.0  to  3.3  eV.  The  performance 
degradation  v-as  determined  using  the  sputter  yield  of  a  Si  target  by  xenon  ions,  however 
it  is  believed  that  the  sputter  yield  of  HfC  and  ZrC  will  be  much  lower  because  these 
materials  have  demonstrated  incredible  robustness  during  ion  bombardment.  HfC  and 
ZrC  cathodes  could  operate  at  much  higher  voltages  than  Si  and  Mo  cathodes  without 
being  damaged  by  ion  bombardment  which  will  further  result  in  fewer  tips  required  to 
meet  lifetime  and  current  requirements.  Table  4.6  shows  that  HfC  cathodes  will  be  able 
to  operate  in  EP  system  environments  with  or  without  CEX  ions  and  with  or  without 
CLAIR. 

Table  4.6  Performance  limitations  of  Si  and  HfC  cathodes  with  and  without  a 
charge-exchange  (CEX)  ion  population. 
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Case 

4>w(eV) 

pc  (rad.) 

re(A) 

r,(A) 

As 

VE(V) 

^tips 

Io(mA) 

If  (mA) 

9-Si 

4.0 

0.26 

2000 

40.4 

i 

52 

6x10s 

105 

104 

10-Si 

4.0 

0.26 

2000 

40.4 

2 

52 

1.2xl06 

106 

106 

11-Sicex 

4.0 

0.26 

2000 

40.4 

2 

30 

8.3xl09 

109 

101 

12-Sicex 

4.0 

0.26 

2000 

30 

2 

30 

2.6x1 08 

100 

92 

13-Si^ 

4.0 

0.26 

1000 

40.4 

2 

30 

2.9x1 08 

109 

102 

14-HfQex 

3.3 

0.26 

1000 

40.4 

2 

30 

3xl06 

106 

100 

15-HfC„x 

3.3 

0.26 

2000 

40.4 

2 

30 

4.6x1 07 

105 

100 

HfC  and  ZrC  cathodes  have  also  demonstrated  greater  stability  in  contaminating 
environments  and  at  higher  currents  than  Mo  and  Si  cathodes.  Chemical  inertness  will  be 
very  important  during  operation  in  Hall  thruster  test  facilities  where  base  pressures  are 
typically  greater  than  10'7  Torr.  With  cathodes  that  are  capable  of  >100  mA/cm2,  the 
emission  current  will  be  limited  by  space-charge  effects.  The  cathode  must  be  designed 
with  the  proper  tip-to-tip  spacing  to  achieve  the  desired  emission  current,  —100  mA. 
Space-charge  limited  emission  is  discussed  in  the  following  chapter. 
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CHAPTER  5 


SPACE-CHARGE  LIMITED  ELECTRON  EMISSION  FROM  FIELD  EMISSION 

CATHODES  INTO  A  PLASMA 

In  this  chapter,  the  development  of  a  sheath  model  is  discussed  with  results  used 
to  predict  sheath  structure  and  emission  current  limitations  for  different  cathode  operating 
configurations.  The  sheath  was  modeled  in  one  dimension  only,  therefore  the  effect  of 
the  potential  of  surrounding  electrodes  on  sheath  structure  was  not  considered.  The 
model  is  based  on  a  planar  sheath  model  that  was  originally  developed  by  Crawford  and 
Cannara54  and  Prewett  and  Allen55  with  a  sheath  boundary  criterion  derived  by  Andrews 
and  Allen56  This  model  and  criterion  were  used  and  augmented  by  Goodfellow  to 
consider  the  initial  electron  thermal  energies  at  the  surface  of  a  hot  cathode.57  In  the 
derivation  of  the  sheath  model  used  in  this  investigation,  initial  electron  energies  were 
used  which  corresponded  to  the  accelerating  potential  between  the  electron  emitting 
surface  and  gate  electrode.  This  sheath  model  was  then  developed  for  cylindrical  and 
spherical  sheath  configurations  also.  One  of  the  objectives  of  this  study  was  to  determine 
the  potential  of  the  gate  electrode,  with  respect  to  the  plasma  potential,  required  to 
transmit  the  desired  electron  current  through  the  sheath  with  various  electron  energies. 
The  details  of  the  pre-sheath  were  not  explored;  however,  it  was  assumed  that  the  ions 
were  accelerated  through  the  pre-sheath  up  to  high  enough  velocities  to  ensure  monotonic 
potential  profiles  in  the  cathode  sheath.  The  Andrews  and  Allen  criterion  was  used  to 
determine  the  ion  energy  at  the  sheath  boundary.  It  is  possible  that  a  virtual  cathode 
forms  in  the  sheath,  where  a  potential  minimum  in  the  sheath  is  lower  than  the  cathode 
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potential.  Virtual  cathode  formation  plasma  instabilities  could  increase  the  upper  limits 
on  the  space-charge  limited  emission  current  and  change  the  structure  of  the  sheath; 
however  these  scenarios  cannot  be  considered  by  the  sheath  model.  This  1-D  sheath 
model  was  used  to  determine  upper  and  lower  limits  on  current  emitted  by  a  planar  FE 
cathode  due  to  space-charge  effects  and  limiting  sheath  geometries  and  is  limited  to 
stable  monotonic  potential  profiles  in  the  sheath.  The  results  of  this  model  are  compared, 
to  the  results  of  a  Particle-in-Cell  (PIC)  model  capable  of  considering  virtual  cathode 
formation. 

This  sheath  model  was  developed  for  planar,  cylindrical,  and  spherical  cathode 
and  sheath  geometries.  A  cathode  and  sheath  potential  diagram  is  shown  in  Figure  1.6 
and  a  possible  potential  profile  is  shown  in  Figure  5.1.  FE  cathodes  are  planar,  however 
the  structure  of  the  sheath  depends  on  the  size  of  the  cathode  and  on  the  plasma 
environment.  With  cathode  dimensions  greater  than  several  electronic  Debye  lengths,  the 
sheath  structure  is  primarily  planar.  If  cathode  dimensions  are  much  smaller  than  the 
sheath  thickness,  the  sheath  is  better  represented  by  a  spherical  sheath  model.  Two 
different  cathode  and  sheath  scenarios  are  illustrated  in  Figure  5.2.  The  cathode 
dimensions  and  environment  significantly  affect  the  emission  current  limitations,  as 
shown  in  the  following  sections  of  this  chapter. 


cathode 


Figure  5.1  Cathode  sheath  configuration  which  was  modeled.  The  boundary 
between  the  sheath  and  pre-sheath  is  shown  at  rsh. 
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rsh 


Figure  5.2  A  planar  cathode  sheath  scenario  (left)  with  r^r^  and  a  spherical 
cathode  sheath  scenario  (right)  with  rc<rsh. 

The  three  environments  considered  in  this  study  for  space-charge  limited  emission 
from  a  FE  source  of  electrons  into  a  plasma  are  those  generated  by  a  Hall  thruster,  a 
mesoscale  ion  thruster,  and  a  tethered  spacecraft  in  low  Earth  orbit.  Since  a  large  Hall 
thruster  could  use  a  FE  cathode  internal  or  external  to  the  thruster  as  shown  in  Figure  1.3, 
both  conditions  were  considered  in  the  analysis.  A  small  or  mesoscale  Hall  thruster  will 
require  an  external  cathode.  Electron  temperatures  in  the  environment  of  a  Hall  thruster 
were  obtained  by  Domonkos  et  al.83  in  the  plasma  plume  of  a  1.4  kW  Hall  thruster. 
Plasma  number  densities  were  approximated  from  Faraday  probe  measurements 
discussed  in  Chapter  2.  These  parameters  were  used  in  the  model  because  initial  testing 
of  these  cathodes  with  Hall  thrusters  will  most  likely  occur  with  such  medium  scale 
systems.  An  ion  engine  requires  two  cathodes;  one  cathode  in  the  discharge  chamber  for 
propellant  ionization,  and  another  cathode  external  to  the  thruster  for  ion  beam 
neutralization.  Plasma  parameters  used  to  represent  the  environment  of  an  ion  engine 
discharge  chamber  are  estimates  only  for  a  mesoscale  ion  engine.  Charged  particle 
number  densities  could  be  much  higher  in  mesoscale  ion  engines  to  get  the  required 
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performance.  There  is  also  interest  in  using  FE  cathodes  on  tethers  in  low  Earth  orbit  to 
emit  electrons  -250-1000  km  above  the  earth.  The  parameters  used  to  represent  each  of 
these  environments  are  presented  in  Table  5.1.  Ambient  electron  number  densities  are 
given  at  an  altitude  of  300  km;  at  higher  altitudes,  electron  number  densities  will  be  even 
lower. 


Table  5.1  Characteristics  of  the  thruster  and  tether  plasma  environments 
investigated  using  the  cathode  sheath  model. 


neo(/cm3) 

Te(eV) 

XD(mm) 

Hall/Ion  thruster  side 

00 

o 

i-H 

X 

00 

5 

0.60 

Hall  thruster  center3 

8xl010 

1 

0.03 

Ion  engine  discharge  chamber 

3x10" 

2-3 

0.02 

Tether  at  250  kmb 

5x10s 

0.1 

3.3 

a  Domonkos,  M.  T.  Marrese,  C.  M.,  Haas,  J.  M.,  Gallimore,  A.  D.,  “Very  Near-Field  Plume  Investigation  of  the  D-55 
Anode  Layer  Thruster”  AIAA-97-3062  Joint  Propulsion  Conference,  Seattle,  WA,  July  1997. 

b  F.  F.  Chen,  Plasma  Physics  and  Controlled  Fusion,  Plenum  Press,  New  York,  1984. 


A  sheath  model  described  in  this  chapter  was  used  to  estimate  space-charge  current  limits 
that  can  be  expected  for  FE  cathodes  in  Hall,  ion  thruster,  and  electrodynamic  tether 
environments  described  in  Table  5.1.  The  effect  of  cathode  operating  voltages,  sheath 
geometry,  and  cathode  dimensions  on  the  current  limit  was  investigated.  Modeling 
results  are  presented  at  the  end  of  this  chapter. 

5.1  Sheath  Model 


5.1.1  Planar  Sheath  Geometry 

A  sheath  model  with  a  planar  cathode  and  sheath  geometry  predicts  conservative  limits 

on  cathode  emission  current  due  to  space-charge  effects.  All  potentials  represented  with 
<b  are  defined  with  respect  to  the  potential  at  the  sheath  boundary  where  quasineutrality  is 
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assumed.  The  pre-sheath  accelerates  the  ions  through  <j)0  to  the  sheath  boundary.  The 
potential  difference  between  the  plasma  and  cathode  gate  electrode  is  <|)c.  Virtual  cathode 

formation  can  not  be  predicted  by  this  model;  it  can  only  predict  monotonic  potential 
profiles.  The  maximum  current  is  associated  with  a  zero  electric  field  at  the  cathode 
surface.  Full  recombination  is  assumed  for  the  ions  at  the  cathode  surface.  Potential 
magnitudes  are  used  in  the  input  aud  output  of  this  model  for  simplicity,  however,  their 
true  signs  are  inherent  to  the  following  equations.  In  this  case,  the  1-D  form  of  Poisson’s 


equation  is 

Eqn.  5.1 


_  P  _  e  t 
dx2  e0  ec  K 


nt  -n 


ee 


The  continuity  and  energy  equations  for  the  ions  and  electrons  in  the  sheath  are  used  to 
determine  particle  number  densities.  The  ion  number  density  at  the  sheath  boundary  is  ni0. 
The  ion  number  density,  ni;  is 


Eqn.  5.2 

ii=w=vvA), 

Eqn.  5.3 

=e(([>  +  <t>0). 

Eqn.  5.4 

(  0  V5 

ni  =  nio  l  +  -r  • 

V  0J 

Plasma  electron  densities  are  defined  by  the  Boltzmann  relationship  for  a  Maxwellian 
distribution  of  electrons 

Eqn.  5.5  ne  =  nt0  ex 


It  is  assumed  that  the  electrons  are  Maxwellian  with  temperature,  Te,  and  the  number 
density  at  the  sheath  boundary  is  neo.  The  beam  electron  number  densities  are  defined  as 
Eqn.  5.6  7„  = 

Eqn.  5.7  i mev2ee  =  e(Vg  +  0c-0), 
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Eqn.  5.8 


2e 


The  condition  of  quasineutrality  at  the  sheath  boundary  is  used  to  define  the  ion  number 
density  at  the  sheath  boundary  in  terms  of  electron  number  densities.  At  <t>=0 


Eqn.  5.9 


nio=neo+nee  ,snd 


Eqn.  5.10 


n.  , 
-!2-  =  l  +  . 


Poisson’s  equation  can  then  be  written  as 


Eqn.  5.11 


d2{-<j>)  _  p  _ g  n  r1+± 


dx  £„  £0 


I7t(^  +<Pc~  ^)"12  “  nto  exPl  ~j^T 


The  parameters  used  to  normalize  Eqn.  5. 1 1  are 


Eqn.  5.12 


ed> 

71  ~W' 


Eqn.  5.13 


£  =  — 
^  2  9 


Eqn.  5.14 


x  Ji&\2 

Uy  ’ 


Eqn.  5.15 


P--2- 

en. 


Eqn.  5.16 


J  2JkT4 1  . 

]c  =  -  .  and 
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Eqn.  5.17 


The  electric  field  at  the  sheath  boundary  is  approximated  at  zero.  This  approximation  is 
common  because  the  electric  field  is  very  small  in  this  region  and  has  little  effect  o ;  the 
structure  of  the  sheath.57,84  The  ion  energy  at  the  sheath  boundary  can  be  approximated 
using  the  criterion  derived  by  Andrews  and  Allen,56 

Eqn.  5.22  -7-  =  0 . 

dn  „=0 

The  ion  energy  at  the  sheath  boundary,  r|0,  is  determined  using  Eqn.  5.18  and  Eqn.  5.22. 
Eqn.  5.23  J70  =  3 

2 +VcP 
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Notice  that  in  the  limit  of  zero  cathode  emission  r|0  is  1/2.  This  result  is  consistent  with 

the  model  developed  by  Bohm.53  The  upper  limit  on  the  emission  current  is  defined  as 
the  current  emitted  when  the  electric  field  at  the  cathode  surface  becomes  zero, 


f 


Eqn.  5.24  0  =  4 r\L 


( 


1  +  ^ 


-1 


+*4  J  (%¥  -  K  +  *»c)5  +  2(exp(-77c)  - 1) . 


With  this  assumption,  the  maximum  emission  current  possible  is 


f 


Vc 


Eqn.  5.25 
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1+a.f 

Vo 


(vg+Vcf-(vg) 
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Eqn.  5.23  and  Eqn.  5.25  can  then  be  solved  simultaneously  to  estimate  J^mar  and  r\0  for 
any  combination  of  t\c  and  qg.  Jeemax  is  the  emission  current  density  that  corresponds  to  a 

zero  electric  field  at  the  cathode  surface.  This  value  represents  the  lower  limit  on  the 
emission  current  density  capability  of  FE  cathodes  emitting  electrons  into  a  plasma.  If 
electron  beam  expansion  occurs  in  the  sheath,  virtual  cathodes  form  in  the  sheath,  or 
plasma  instabilities  develop,  larger  emission  currents  can  be  possible.  Emission  current 
limitations  predicted  by  this  model  agree  with  results  of  the  Goodfellow  sheath  model, 
which  provided  the  basis  for  this  model.  In  the  case  considered,  Tjc  was  10,  T|g  was  0.3, 

and  both  models  predicted  that  JMmax  was  0.7.  Goodfellow  used  a  similar  sheath  model 
combined  with  thermal  models  to  predict  cathode  temperatures  in  a 
magnetoplasmadynamic  thruster  system,  and  achieved  good  correlation  between 
experiment  and  theory.  The  results  of  the  1-D  sheath  model  derived  in  this  thesis  work 
were  also  compared  with  the  Prewett  and  Allen  model  results.  A  comparison  of  the 
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results  from  both  models  shown  in  Figure  5.3.  The  Prewett  and  Allen  model  assumes 
negligible  initial  velocities,  and  an  rig  value  of  10‘8  is  used  in  the  model  developed  in  this 

thesis  work  to  obtain  the  results  shown  in  Figure  5.3.  The  results  of  this  model  are 
consistent  with  results  from  the  Andrews  and  Allen  model  and  the  Goodfellow  model. 


Figure  5.3  Results  of  the  cathode  sheath  model  compared  to  the  Prewett  and  Allen 
sheath  model  results. 

5.1.2  Cylindrical  Geometry 

The  same  sheath  model  was  also  developed  for  a  cylindrical  cathode  and  sheath 
geometry.  This  model  assumes  a  1/r  expansion  of  the  emitted  electron  current  in  the 
cathode  sheath  and  consequent  decrease  in  ion  number  density.  As  the  sheath  dimensions 
become  comparable  to  the  cathode  dimensions  some  radial  beam  expansion  occurs  and 
the  cylindrical  sheath  model  becomes  more  appropriate  than  the  planar  model.  In  the 
cylindrical  sheath  model,  the  cathode  is  emitting  electrons  radially  outward  into  a  plasma 
where  the  electron  number  density  is  ne0.  Plasma  electron  number  density  is  given  by  the 
Boltzmann  equation,  assuming  a  Maxwellian  distribution  of  energies, 
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r  ±\ 


Eqn.  5.26 


n<=n<0Qm 


\u.j 


The  number  density  of  electrons  emitted  from  the  cathode  is  derived  using  the  continuity, 
Eqn.  5.27  jjpre)  =  naevte{nr) , 


and  energy  equations, 
Eqn.  5.28 


v«  = 


+ &  -  <p) 


i 

■M 


m. 


The  electron  number  density  for  the  beam  electrons  is  then 


Eqn.  5.29 


2e 


The  continuity  and  energy  equations  are  used  to  determine  the  ion  number  density, 


Eqn.  5.30 
Eqn.  5.31 

Eqn.  5.32 


h  =  nfivfar)  =  nioevio(7trsh). 


v.  - 


2  g(0o  +  0) 


V  mi  J 
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f  ^  Ur  ^ 
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1+-I  I*, 

V  t)  Vr; 


The  initial  ion  energy  at  the  sheath  boundary  is  e<J>0,  and  the  ion  number  density  at  the 
sheath  boundary  is 


Eqn.  5.33 


nio=neo+nee=neo 


1+ — — + 


'2l 


(  2e3V2 


where  quasineutrality  is  assumed.  The  ion  number  density  in  the  sheath  is 
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Eqn.  5.34 


»,=«„ 1+  /" 
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Poisson’s  equation  can  then  be  written  for  this  cathode  configuration  and  environment  as 

Eqn.  5.35  _v  (  (t>)~~\^  dj.i  +  r  dr  J  £g 


“f  41  +  £  • 


Eqn.  5.12-Eqn.  5.17  are  used  to  get  the  normalized  form  of  Poisson’s  equation 

describing  the  potential  profile  in  the  cylindrical  sheath, 

d2r]  1  dr]  - 

Eqn.  5.36  d^+^d^~P 


=  i + Jm(v,  +  ffcpf  fj  [j )  -  +7?c-  n)\ j  j  ~  e*p(-v)- 


The  Andrews  and  Allen  criterion56  is  evaluated  at  tho  sheath  boundary  to  derive 
the  relationship  for  the  normalized  potential  drop  across  the  pre-sheath, 

i +JM(nt  +  ncp{j- 

_  „  \SshJ 


Eqn.  5.37 


2-4(^  +  nc)l^ 


Eqn.  5.36  is  evaluated  numerically  using  a  Runge-Kutta  fourth  order  method.  The  initial 

conditions  at  the  sheath  boundary  are  given  by 
Eqn.  5.38  *?(£  =  L)  =  0  and 


Eqn.  5.39 


|f« =&.>=»• 
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The  cylindrical  cathode  sheath  problem  is  not  very  well  defined,  since  only  the  potential 
at  the  cathode  and  sheath  boundaries,  and  the  electric  field  at  the  sheath  boundary  are 
known.  The  position  of  the  sheath  boundary  is  an  input  to  the  model,  however,  it  has  to 
be  estimated  initially.  The  much  simpler  nature  of  the  equations  for  charged  particle 
densities  in  the  planar  sheath  geometry  does  not  require  this  input  parameter.  The  sheath 
thickness  is  a  result  of  the  planar  sheath  model;  therefore,  the  structure  of  the  sheath  can 
be  determined  for  any  J^,  r\c,  and  rig  combinations.  In  cylindrical  coordinates  the  Poisson 

equation  is  more  complicated,  requiring  more  input  parameters.  Initial  estimates  are 
made  for  the  thickness  of  the  sheath,  ^h.  It  is  assumed  that  the  maximum  current  is 

emitted  when  the  electric  field  at  the  cathode  surface  is  zero  and  the  potential  at  the 
cathode  surface  is  equivalent  to  the  cathode  potential.  The  sheath  thickness  is  adjusted 
until  these  conditions  are  met  at  the  cathode  surface  to  estimate  JKmax.  A  potential  profile 
for  the  sheath  can  only  be  obtained  for  limiting  cases.  Potential  profiles  and  Jeemax  are 
obtained  in  this  model  for  the  limiting  case  of  a  zero  electric  field  at  the  cathode  surface. 


5.1.3  Spherical  Geometry 


The  sheath  model  with  a  spherical  geometry  was  derived  like  the  cylindrical 
model  except  that  the  ion  and  electron  number  densities  change  with  a  1/r2  relationship  in 
the  sheath.  Poisson’s  equation  in  spherical  coordinates  for  this  cathode,  sheath 


configuration,  and  environment  is 


Eqn.  5.40 
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The  pre-sheath  accelerates  the  ions  through  a  normalized  potential  r)0 
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Eqn.  5.41 
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as  derived  using  the  Andrews  and  Allen  criterion.  Eqn.  5.40  is  evaluated  numerically 
using  a  Runge-Kutta  fourth  order  method  for  a  second  order  differential  equation.  The 
initial  conditions  at  the  sheath  boundary  are  Eqn.  5.38  and  Eqn.  5.39.  An  upper  limit  on 
jge  is  determined  using  the  spherical  sheath  model.  Initial  guesses  are  made  for  the 
thickness  of  the  sheath,  r^.  It  is  assumed  that  the  maximum  current  is  emitted  when  the 
electric  field  at  the  cathode  surface  is  zero  and  the  potential  at  the  cathode  surface  is 
equivalent  to  the  cathode  potential.  The  sheath  thickness  is  adjusted  until  these 
conditions  are  met  at  the  cathode  surface. 


5.2  Modeling  Results 


Space-charge  limited  FE  was  investigated  for  three  geometrical  sheath 
configurations  in  one  dimension.  The  planar  geometry  was  the  simplest  geometry 
considered.  The  results  of  this  model  are  most  accurate  when  the  sheath  thickness  is 
much  smaller  than  the  cathode  dimensions.  As  the  cathode  dimensions  become  smaller 
or  the  ambient  ion  and  electron  number  densities  decrease,  the  sheath  becomes  less  planar 
and  some  radial  beam  expansion  begins  to  occur  in  the  cathode  sheath.  When  the  sheath 
thickness  becomes  several  times  larger  than  the  cathode  radius,  the  sheath  approaches  a 
spherical  geometry.  Lower  energy  emitted  electrons  will  support  a  more  spherical  sheath 
configuration  because  the  electron  beam  will  be  less  collimated  in  the  sheath.  Local 
radial  electric  fields  generated  by  high  negative  charge  densities  in  the  sheath  will 
generate  more  expansion  for  a  lower  energy  beam.  With  cylindrical  and  spherical  sheath 
geometries,  ion  number  densities  increase  in  the  sheath  as  1/r  and  1/r2,  respectively, 
providing  much  better  negative  space-charge  neutralization  at  the  cathode  surface  by  the 
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plasma  ions  than  the  planar  geometry  to  permit  higher  emitted  electron  currents.  The 
cylindrical  and  spherical  sheath  models  were  used  to  study  the  effect  of  cathode 
dimensions  on  FE  current  limitations,  an  investigation  not  possible  with  the  planar  model. 

The  planar  model  provided  conservative  estimates  of  upper  cathode  current  limits 
set  by  space  charge  effects  for  a  monotonic  potential  distribution.  Virtual  cathode 
formation  and  some  plasma  instabilities  will  increase  the  upper  limits  on  the  emission 
current.  The  results  of  this  model  are  shown  in  Table  5.2.  Current  density  objectives  for 
tether  and  thruster  environments  are  -100  mA/cm2  from  a  1  cm2  cathode.  The  results  of 
the  planar  model  show  that  the  conservative  estimates  of  lower  limits  on  the  emission 
current  greatly  exceed  the  requirements  when  the  cathodes  are  used  in  the  discharge 
chamber  of  an  ion  thruster  or  at  the  internal  position  of  a  Hall  thruster.  The  sheath  model 
results  show  that  FE  cathodes  emitting  into  a  tether  environment  are  prohibitively  limited 
by  space-charge  effects  if  a  planar  sheath  geometry  is  assumed. 


Table  5.2  Space-charge  limited  currents  estimated  by  the  planar  sheath  model  for 
thruster  and  tether  environments. 


neo  (/cm3) 

Te  (eV) 

v,(V) 

4>c(V) 

jeemax  (mA/CHl2) 

Hall  jiruster-extemal 

8x10s 

5 

8x10s 

5 

68 

8xl08 

5 

34 

8x10s 

1 

20 

17 

Hall  thruster-internal 

8xl010 

1 

8xl010 

1 

mM 

2,884 

Ion  engine  discharge 
chamber 

3xlOn 

2 

30 

20 

8,800 

Tether 

5x10s 

0.1 

30 

20 

0.003 

5x10s 

0.1 

100 

20 

0.005 

5x10s 

0.1 

1000 

20 

0.016 

5x10s 

0.1 

30 

100 

0.002 

The  effects  of  Vg  and  <J)C  on  the  current  limitations  were  also  considered.  The 
results  of  this  investigation  are  also  shown  in  Figure  5.4  and  Figure  5.5.  According  to 
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these  model  results,  cathodes  will  require  operating  voltages,  -200  V,  which  greatly 
exceed  the  sputter  threshold  of  the  Si  or  Mo  cathodes  if  used  external  to  the  Hall  or  ion 
thruster  where  n^  is  8x10s  /cm3  and  Te  is  5  eV.  The  effect  of  increasing  the  electron 
energy  with  the  gate  voltage,  Vg,  is  shown  in  Figure  5.4.  Increasing  the  electron  energies 

at  the  gate  electrode  does  increase  the  emission  current  from  the  cathode.  Increasing  the 
potential  between  the  cathode  and  plasma,  <j>c,  either  increases  or  decreases  the  emission 

current  limit,  depending  on  the  initial  electron  energies,  as  represented  by  the  gate 

electrode  voltage, Vg.  Increasing  this  potential  only  decreased  the  cathode  current  limit 
for  T|g  >2,  as  shown  in  Figure  5.5.  Increased  potential  drops  in  the  cathode  sheath  do  not 

increase  the  electric  field  proportionally  because  increasing  <)>c  also  increases  the  size  of 
the  sheath.  Increasing  <|>c  with  a  planar  sheath  geometry  also  decreases  the  ion  number 

density  at  the  gate  electrode  because  the  ions  are  accelerated  through  a  larger  potential 
between  the  sheath  boundary  and  gate  electrode.  This  effect  reduces  negative  space- 
charge  neutralization  by  the  ions.  FE  cathode  gate  electrode  potentials  should  not  be 
much  more  than  20  V  below  the  plasma  potential  to  maximize  the  current  emitted  from 
the  cathode  and  minimize  the  thermalized  electron  current  collected  by  the  gate  electrode. 
A  gate  electrode  which  is  no  more  than  20  V  below  plasma  potential  is  also  required  to 
avoid  excessive  cathode  sputtering  by  the  local  ions,  as  shown  in  the  results  of  the 
cathode  sputtering  model  discussed  in  the  previous  chapter. 
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Figure  5.4  The  effect  of  Vg  on  for  the  planar  sheath  geometry.  The  numeric 
markers  represent  the  value  of  t]0  for  each  solution. 


Figure  5.5  The  effect  of  <t>c  on  for  the  planar  sheath  geometry.  The  numeric 
markers  represent  the  value  of  Tl0  for  each  solution. 


Ill 


Potential,  electric  field,  and  charge  density  profiles,  as  predicted  by  the  planar 
cathode  sheath  model,  are  shown  in  the  following  figures  for  the  Hall  thruster 
application.  Figure  5.6  shows  potential  and  electric  field  profiles  in  the  cathode  sheath 
for  some  operating  configurations.  Charge  density  profiles  are  shown  in  Figure  5.7.  The 
effect  of  Jge  on  the  sheath  structure  is  shown.  The  charge  density  profiles  show  the 
double  sheath  structure  for  both  Jee=1.5,  and  2.0,  where  the  sheath  consists  of  beta 
positive  and  negative  space-charge  regions. 


Figure  5.6  a,  b  Potential  (a)  and  electric  field  (b)  profiles  for  different  emission 
currents  in  the  cathode  sheath. 


112 


Figure  5.7  Charge  density  profiles  in  a  planar  cathode  sheath  with  T]c=4  and  t]g=6. 

In  the  Hall  thruster  configuration  where  the  cathode  is  external  to  the  thruster,  the 
sheath  dimensions  are  comparable  to  the  cathode  dimensions,  which  are  assumed  to  be 
lxl  cm.  Therefore,  the  sheath  would  be  better  represented  by  a  cylindrical  geometry  that 
considers  some  radial  electron  beam  expansion.  In  the  tether  environment,  plasma 
densities  are  so  low  that  the  cathode  sheath  dimensions  can  be  much  larger  than  the 
cathode  dimensions,  necessitating  the  use  of  a  sheath  model  that  is  more  spherical  in 
geometry  with  a  1/r2  change  in  number  densities  across  the  sheaths.  A  comparison  of  the 
results  of  the  models  of  three  different  sheath  geometries  is  presented  in  Table  5.3.  The 
planar  model  can  be  used  to  predict  lower  limits  on  the  cathode  current.  The  spherical 
model  should  predict  the  upper  limits  on  the  current  that  can  be  emitted  from  the  cathode 
because  this  sheath  configuration  will  assume  the  highest  negative  space-charge 
neutralization  at  the  cathode  surface.  The  ion  collection  area  of  the  sheath  boundary  will 
be  much  larger  with  a  spherical  geometry  than  with  a  planar  or  cylindrical  geometry.  It  is 
uncertain,  however  improbable,  that  virtual  cathode  formation  can  not  supply  the  same 
advantage  to  electron  beam  emission  as  a  spherical  sheath  geometry. 
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Table  5.3  Space-charge  limited  currents  predicted  by  the  planar,  cylindrical,  and 
spherical  sheath  models  in  normalized  units. 


geometry 

Sc 

Ssh 

Tic 

Tlo 

Jeemax 

planar 

20.0 

4 

6 

cylindrical 

40 

54.0 

4 

6 

0.87 

2.9 

cylindrical 

20 

32.0 

4 

6 

0.87 

3.5 

cylindrical 

10 

19.4 

4 

6 

0.85 

4.0 

cylindrical 

4 

13.0 

4 

6 

0.85 

6.6 

cylindrical 

1 

8.8 

4 

6 

0.80 

15.3 

cylindrical 

1 

18.4 

20 

6 

0.55 

9.8 

spherical 

40 

48.6 

4 

20 

1.34 

31.0 

spherical 

40 

47.0 

4 

6 

1.34 

6.5 

spherical 

20 

4 

6 

1.22 

7.4 

spherical 

10 

4 

6 

1.10 

9.9 

spherical 

4 

4 

6 

0.95 

17.6 

spherical 

1 

6.8 

4 

6 

0.76 

70.0 

spherical 

0.4 

5.4 

4 

6 

0.65 

168.0 

The  environment  investigated  for  the  results  shown  in  Table  5.3  was  that  of  a  Hall 
thruster  with  Te=5  eV  so  that  Vc=20  V  and  Vg=30  V.  Vc  was  assumed  to  be  20  V  because 
it  is  common  for  hollow  cathodes  to  float  20  V  below  the  local  plasma  potential  to  retard 
the  collection  of  plasma  electrons.  The  sheath  and  cathode  geometry  strongly  influences 
the  resuits  of  the  models,  as  data  in  Table  5.3  show.  The  planar  model  predicts  lower 

jL 

current  density  limits.  The  cylindrical  and  spherical  geometry  sheath  models  predict 
current  limits  that  approach  the  planar  model  results  for  large  cathode  configurations. 

Results  from  the  sheath  models  for  cathodes  operating  in  the  external  cathode 
position  of  a  Hall  thruster  are  presented  in  Table  5.4  for  different  sheath  geometries. 
While  cathodes  can  emit  much  higher  current  densities  into  spherical  sheaths,  spherical 
sheath  formation  requires  that  rc«rsh.  With  such  small  emitting  cathode  areas,  total 
cathode  current  is  sacrificed  for  greater  current  density  as  shown  in  Table  5.4.  The  area 
of  the  cathode  is  determined  by  assuming  a  rectangular  (cylindrical  sheath)  or  square 
(spherical  sheath)  and  that  one-half  of  the  length  of  one  side  is  rc.  The  most  favorable 
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cathode  configuration  for  this  application  is  a  ~1.2  cm  cathode  with  jK  at  82  mA/cm. 
Therefore,  with  a  1.2  cm  cathode  length,  Iw  will  be  98  mA.  Many  of  these  cathodes 
could  be  used  to  provide  a  few  Amperes  of  current  required  by  a  1.35  kW  system  or  1-2 
could  be  used  for  a  50  W  thruster. 


Table  5.4  Field  emission  current  limitations  predicted  by  the  sheath  model  in 
dimensional  parameters  for  planar,  cylindrical,  and  spherical  geometries  for  the 
Hall  thruster  environment  with  an  external  cathode.  Vg  is  at  30  V  and  Vc  is  at  20  V. 


geometry 

rc  (mm) 

jee  (A/Cm2) 

jM  (mA/cm) 

Iee(mA) 

planar 

34 

82 

cylindrical 

10 

6 

68 

4 

2.4 

112 

54 

1 

0.6 

260 

31 

spherical 

4 

2.4 

299 

69 

1 

0.6 

1190 

17 

Space-charge  current  limitations  are  also  predicted  using  the  sheath  model  to 
study  a  cathode  emitting  from  an  electrodynamic  tether  system  into  a  Low  Earth  Orbit 
(LEO)  environment.  Results  with  non-dimensional  parameters  for  a  spherical  sheath  are 

shown  in  Table  5.5.  This  table  shoyvs  that  the  sheath  dimensions  are  so  large  that  the 
sheath  geometry  could  be  spherical  even  when  ^  is  10.  In  the  tether  environment  Te  is 

0.1  eV,  therefore  r|g  is  300  and  r|c  is  200  when  Vg  is  30  V  and  Vc  is  20  V.  The  current 

density  emitted  by  the  cathode  greatly  increases  with  decreasing  cathode  dimensions. 
Table  5.6  shows  how  the  total  current  emitted  by  the  cathode  is  affected  by  the  size  of  the 
cathode.  Smaller  cathodes  are  capable  of  larger  current  densities,  however,  the  emitting 
area  is  also  considerably  smaller.  Therefore,  there  is  no  inherent  advantage  in  reducing 
cathode  size  in  terms  of  emission  characteristics.  Also  shown  in  Table  5.5  and  Table  5.6 
is  the  current  gain  attributable  to  increases  in  the  energy  of  the  emitted  electrons  by 
increasing  the  potential  of  the  gate  electrode.  Increasing  the  gate  electrode  from  30  V  up 
to  100  V  only  increases  the  current  by  a  factor  of  -2.  To  satisfy  the  -2  A  current 
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requirements  of  the  tether  application,  several  hundred  cathodes  will  be  required 
according  to  these  model  results. 


Table  5.5  Field  emission  current  limitations  predicted  by  the  sheath  model  in 
normalized  parameters  for  a  spherical  sheath  in  a  tether  environment  where  Te  is 
0.1  eV. 


4c 

4* 

Tic 

Tie 

Jeemax 

10 

64 

100 

100 

0.53 

40 

10 

100 

200 

200 

0.52 

91 

9 

103 

200 

300 

0.52 

145 

5 

100 

200 

300 

0.52 

430 

1 

97 

200 

300 

0.52 

8650 

1 

99 

200 

400 

0.52 

10400 

1 

101 

200 

600 

0.52 

13400 

1 

104 

200 

1000 

0.52 

18400 

1 

109 

200 

10000 

0.53 

63550 

1 

60.5 

100 

300 

0.53 

4925 

1 

37 

50 

300 

0.55 

2750 

1 

38 

50 

1000 

0.56 

5350 

1 

38.5 

50 

10000 

0.56 

17400 

Table  5.6  Field  emission  current  limitations  predicted  by  the  spherical  sheath 
model  for  the  tether  environment. 


4c 

rc(mm) 

VC(V) 

V„(V) 

I«  (mA) 

5 

16.5 

20 

30 

/.  0 

1 

3.3 

20 

30 

5.6 

1 

3.3 

20 

40 

6.8 

1 

3.3 

20 

60 

8.8 

1 

3.3 

20 

100 

12.0 

1 

3.3 

20 

1000 

41.5 

1 

3.3 

10 

30 

3.2 

1 

3.3 

5 

30 

1.8 

1 

3.3 

5 

100 

3.5 

1 

3.3 

5 

1000 

11.3 

Unfortunately,  these  modeling  results  have  yet  not  been  substantiated  by  experimental 
observations.  The  trends  observed  are  consistent  with  modeling  results  obtained  in  other 
sheath  studies  that  increasing  sheath  potential  increases  the  size  of  the  sheath,  and  that 
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increasing  the  electron  current  emitted  from  the  cathode  increases  the  size  of  the  sheath.57 
The  enormous  size  of  the  cathode  sheath  in  the  tether  environments  was  especially 
surprising;  however,  the  CHARGE  2  sounding  rocket  experiment  coupled  with 
NASCAP/LEO  computer  simulations85  did  predict  sheath  dimensions  similar  to  the 
dimensions  predicted  by  the  sheath  model  discussed  in  this  chapter.  '  ", 

The  sheath  structure  around  a  rocket  payload  charged  up  to  460  V  negative  to  the 
ambient  ionospheric  plasma  at  245  km  altitude  was  studied  experimentally  and  by 
computer  simulations  of  the  CHARGE  2  sounding  rocket  experiment.  The  payload  was 
split  into  two  sections  connected  by  a  conducting  tether.  The  mother  section  was  biased 
to  potentials  as  much  as  450  V  below  the  daughter  section,  which  was  close  to  ambient 
plasma  potential.  The  mother  section  carried  four  floating  probes  on  a  boom,  which 
extended  out  from  the  spacecraft  with  probes  at  25,  50,  75,  and  100  cm.  In  this 
investigation,  electron  number  densities  were  assumed  to  be  4xl04  cm'3,  electron 
temperatures  0.086  eV,  and  ion  temperatures  0. 1  eV.  The  sheath  boundary  in  that  study 
was  defined  at  -0.4  V.  Theoretical  and  experimental  results  showed  that  the  cylindrical 
sheath  thickness  was  60  cm  when  the  spacecraft  potential  was  -30  V.  The  cylindrical 
sheath  model  discussed  in  this  chapter  cannot  be  used  to  estimate  sheath  thickness 
without  the  electric  field  at  the  cathode  surface  being  known  a  priori ;  however,  the  planar 
sheath  model  could  be  used,  assuming  zero  electron  emission,  and  the  thickness  of  the 
sheath  could  be  adjusted  for  a  cylindrical  sheath.  The  planar  sheath  model  was  used  to 
predict  a  sheath  thickness  of  79  cm  if  the  sheath  boundary  was  defined  at  -0.4  V  and  the 
sheath  voltage  was  -30  V  with  Te  at  0.086  eV  and  neo  at  4xl04cm3.  This  correlation  was 
fairly  good  considering  that  a  planar  sheath  model  will  over-predict  the  thickness  of  a 
cylindrical  sheath.  A  cylindrical  sheath  around  a  negatively  charged  electrode  requires  a 
smaller  thickness  than  a  planar  sheath  because  of  the  amount  of  positive  charge  contained 
within  the  sheath  volume.  A  cylindrical  sheath  with  a  thickness  of  68  cm  will  contain  the 
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same  volume  of  charge  as  a  planar  sheath  with  a  thickness  of  79  cm.  The  sheath  model 
effectively  predicted  the  sheath  dimensions  in  this  case  with  an  error  of  13%. 

The  modeling  results  discussed  in  this  chapter  were  also  compared  to  the  results 
of  a  three-dimensional  particle-in-cell  simulation  model  developed  by  Wang86.  The  one¬ 
dimensional  sheath  model  can  be  used  to  obtain  almost  immediate  results,  while  the 
three-dimensional  PIC  model  is  computationally  expensive.  The  Wang  simulation  model 
considers  three-dimensional  beam  expansion,  determines  the  true  structure  of  the  sheath, 
and  can  evaluate  the  effect  of  virtual  cathode  formation  on  space-charge  limited  electron 
emission  into  a  plasma.  His  model  was  originally  developed  to  study  virtual  anode 
formation  in  ion  beam  emission  in  space87  and  was  later  modified  to  investigate  electron 
beam  emission  into  a  plasma.86  It  is  capable  of  resolving  interactions  between  beam 
electrons  and  ambient  plasma  with  a  full  particle,  three-dimensional  electrostatic  PIC 
code,  which  follows  both  the  beam  electrons  and  the  plasma  ions  and  electrons. 

Results  of  the  three-dimensional  PIC  simulation  model  and  the  one-dimensional 
sheath  model  were  compared  for  scenarios  with  cathodes  operating  in  thruster  and  tether 
environments.  The  results  from  the  Wang  model  for  a  FE  cathode  operating  in  a  thruster 
environment  are  shown  in  Table  5.7.  The  current  emitted  from  the  cathode  is  1^,  while 
the  current  transmitted  through  the  sheath  is  1^.  The  minimum  potential  in  the  cathode 

sheath  with  respect  to  the  potential  at  the  sheath  boundary,  the  virtual  cathode  potential, 
is  Note  how  a  virtual  cathode  potential  at  -50  V  in  the  sheath  does  not  prevent  the 

total  emitted  current  from  escaping  through  the  sheath.  However,  with  Ire  at  52  mA/cm2 
in  this  case,  a  virtual  cathode  forms  with  a  potential  at  -75  V  which  limits  the  transmitted 
current  to  34  mA.  In  this  scenario,  18  mA  would  be  collected  by  the  cathode  gate 
electrode,  which  is  intolerable  with  FEA  cathodes.  Some  of  this  current  will  be  absorbed 

by  the  flux  of  ions  to  the  gate  electrode;  however,  the  current  to  this  thin  film  gate 
electrode  is  limited  to  only  -10-20  pA.  The  electron  emission  current  density  upper  limit 

predicted  by  the  Wang  model  in  this  cathode  and  plasma  configuration  is  between  125 
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and  236  mA/cm2.  The  current  density  limitations  predicted  by  the  1-D  sheath  model  of 
planar,  cylindrical,  and  spherical  sheath  geometries  are  34, 112  and  299  mA/cm  , 
respectively.  The  Wang  model  showed  that  virtual  cathode  formation  does  not  provide  a 
current  advantage  greater  than  that  provided  by  the  ideal  spherical  sheath  geometry ,  for  a 
planar  cathode  configuration.  The  3-D  Wang  model  predicted  upper  limits  on  the 
emission  current  density  that  were  between  the  upper  and  lower  limits  predicted  by  the 
1-D  sheath  model,  closest  to  the  current  limit  predicted  using  the  cylindrical  sheath 
model. 


Table  5.7  Current  limitations  and  virtual  cathode  potentials  predicted  by  the  Wang 
simulation  model  for  a  field  emission  cathode  in  a  Hall  thruster  environment  with 
r  =  4  Xo,  Vg  at  30  V,  <(>c  at  -20  V,  Tt  at  5  eV,  and  nTO  at  8x107cm3. 


j«  (mA/cm2) 

Iee(mA) 

*™(V) 

I^fmA) 

^trans^ce 

34 

7.5 

7.5 

1 

84 

18.5 

-34 

18.5 

l 

125 

27.5 

27.5 

1 

236 

52 

-75 

34 

0.67 

Wang  also  used  his  3-D  PIC  model  to  study  electron  emission  from  FE  cathodes 
into  a  space  tether  environment.  The  cathode  and  spacecraft  potential  configuration 
affects  the  emission  current  limitations;  however  these  details  were  not  investigated  with 
this  model.  In  the  configuration  investigated,  the  spacecraft  was  at  the  cathode  potential. 
The  results  are  shown  in  Table  5.8.  In  this  case,  a  current  density  of  13  mA/cm2  is 
unacceptable  because  of  the  prohibitively  high  current  that  will  be  collected  by  the  gate 
electrode.  The  current  density  upper  limit,  as  predicted  by  that  model  is  between  6.5  and 
13  mA/cm2.  The  upper  limit  on  the  current  density  predicted  by  the  spherical  1-D  sheath 
model  is  13  mA/cm2.  Considering  virtual  cathode  formation  in  cathode  sheaths  with  a  3- 
D  PIC  model  leads  to  current  density  limit  estimates  below  those  made  by  the  1-D 
spherical  sheath  model.  According  to  the  3-D  PIC  model,  thel-D  spherical  and  planar 
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sheath  models  do  provide  upper  and  lower  limits  on  the  electron  current  density  that  can 
be  expected  from  planar  FE  cathodes  in  plasma  environments,  which  was  the  objective 
for  the  simple  model  development.  The  advantage  of  the  1-D  sheath  model  is  that  results 
can  be  obtained  almost  immediately.  The  3-D  PIC  model  is  computationally  expensive, 
requiring  approximately  two  hours  to  execute.  To  study  the  effect  of  the  potentials  of 
surrounding  electrodes,  the  interaction  of  multiple  electron  beams,  and  to  obtain  more 
accurate  predictions  of  the  emission  current  limitations  from  FE  cathodes  in  plasma 
environments,  the  3-D  PIC  simulation  model  is  required. 

Table  5.8  Current  limitations  and  virtual  cathode  potentials  predicted  by  the  Wang 
simulation  model  for  a  field  emission  cathode  operating  in  a  tether  environment 
with  rc=  Xo,  Vg  at  30  V,  at  -20  V,  Te  at  0.1  eV,  and  nro  at  5xlOs/cm3. 


j«  (mA/cm2) 

Iee(mA) 

(V) 

Ilians  (mA) 

Itrans^ee 

6.5 

2.86 

-42 

2.86 

1 

13 

5.7 

-49 

4.9 

26 

11.4 

-60 

5.8 

EH 

This  chapter  presented  a  sheath  model  which  is  capable  of  predicting  upper  and 
lower  limits  on  the  emission  currem  density  that  can  be  expected  from  a  FE  cathode  in 
plasma  environments.  It  can  be  used  to  estimate  the  maximum  and  minimum  number  of 
cathodes  required  to  supply  the  desired  current  in  a  plasma  environment.  This  model  was 

used  to  study  the  effects  of  initial  electron  energies,  eVg,  gate  electrode  potential  below 
the  potential  at  the  plasma  sheath  boundary,  <f>c,  and  rc  on  emission  current  limitations. 

Increasing  Vg  increases  the  jcc  limit.  While  decreasing  rc  increases  jce,  Icc  from  the  cathode 
decreases.  The  optimal  Vc  value  depends  on  the  cathode  operating  configuration  and 
environment;  could  increase  or  decrease  with  increasing  Vc. 

The  model  was  used  to  assess  the  compatibility  of  FE  cathode  and  electric 
propulsion  and  tether  systems.  Modeling  results  showed  that  the  FE  cathode  current 
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should  not  be  space-charge  limited  within  the  required  current  regime  if  operated  in  the 
discharge  chamber  of  an  ion  engine  or  internal  position  of  a  Hall  thruster.  In  the  external 
position  of  a  Hall  or  ion  thruster,  a  cathode  area  of  40  cm2  is  required  if  Vg  is  30  V  for  4.5 
A.  This  current  could  be  emitted  from  a  segmented  cathode  ring  around  the  thruster. 
Cathode  area  can  be  decreased  by  increasing  Vg.  Much  sma’ler  cathode  areas  are 
required  for  smaller  power  systems.  The  size  of  the  cathodes  required  is  not  prohibitively 
large.  The  results  of  the  sheath  model  for  some  possible  cases  show  that  the  tether 
application  will  require  a  prohibitively  large  number  of  cathodes,  or  cathode  area,  if 
electrons  are  emitted  at  30  eV.  A  reasonable  number  of  cathodes  is  required  if  Vg  is  1000 
V.  However,  it  is  not  possible  to  operate  FEA  cathodes  with  Vg  at  1000  V,  because  of  the 
breakdown  strength  of  the  insulator  material  between  the  gate  electrode  and  tip. 
Acceleration  of  the  electron  beam  from  30  eV  up  to  1000  eV  could  occur  downstream  of 
Vg  by  another  electrode.  The  CLAIR  structure  disclosed  in  Appendix  B  could  be  used  to 
decouple  the  electron  energy  from  the  gate  electrode  potential  and  accelerate  the  electron 
beam  by  an  electrode  downstream  of  the  gate  electrode. 


121 


CHAPTER  6 


CONCLUSIONS 

The  objective  of  this  thesis  was  to  gain  a  better  understanding  of  the  effects  of  EP 
system  environments  on  FE  cathode  performance,  address  compatibility  issues,  and 
develop  models  to  predict  performance  limitations  which  can  then  be  used  in  cathode 
design  and  accelerated  life  tests.  Sheath  models  developed  were  also  used  to  study 
space-charge  limitations  of  FE  cathodes  in  tether  environments.  The  results  of  this 
research,  possible  cathode  configurations  for  the  thruster  and  tether  applications  of 
interest,  and  recommendations  for  future  research  are  discussed  in  this  concluding 
chapter. 

6.1  Summary  of  Major  Results 

Initial  cathode  experiments  in  xenon  environments  were  based  on  experimental 
results  obtained  in  research  laboratories  at  SRI  International  and  MCNC  and  cathode 
environment  characterization  of  the  PEPL-70.  The  results  obtained  from  experiments 
with  single  HfC  and  ZrC  cathode  tips  in  ~10'6Torr  Ar  and  02  environments  showed  that 
these  materials  are  incredibly  robust.  The  cathodes  were  not  destroyed  in  hostile 
environments;  often  performance  improved  from  the  high-pressure  exposures  with 
operating  voltages  at  a  few  thousand  volts.  Experiments  with  Mo  and  Si  FEA  cathodes, 
and  carbon  film  FE  cathodes  showed  that  no  work  function  increases  resulted  from  xenon 
exposure.  The  efficiency  of  these  cathodes  was  not  affected  by  increases  in  pressure  up 
to  2x1  O'5  Torr  until  the  tips  were  damaged  by  ion  bombardment.  Si  and  Mo  FEA 
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cathodes  were  operated  in  xenon  environments  with  pressures  as  high  as  10 5  Tore  for 
several  hours  without  measurable  performance  degradation.  Experimental  results  showed 
that  FEA  cathode  performance  in  xenon  pressure  environments  similar  to  Hall  thruster 
environments  is  very  sensitive  to  operating  voltages.  Energy  threshold  values  for 
sputtering  Mo  and  Si  targets  by  xenon  ions  were  determined  to  be  much  lower  than  the 
model  developed  by  Yamamura  et  al.  predicted.  Experimental  results  also  showed  that 
the  SRI  International  cathode,  fabricated  on  a  high  resistivity  wafer,  was  more  resistant  to 
arcing  in  high-pressure  environments  than  other  cathode  configurations  tested.  The 
carbon  film  cathode  was  incredibly  robust,  showing  no  performance  degradation  during 
several  hours  in  2x1  O'6  Tore  of  xenon  with  the  gate  electrode  at  800  V.  This  cathode 
turned  on  no  differently  in  2xl0‘6  Tore  of  xenon  and  in  UHV  environments. 

Combining  the  performance  degradation  model  with  experimental  results  enabled 
much  more  accurate  predictions  of  E*  than  ever  before  available.  The  for  a  Mo  target 
being  sputtered  by  Xe  ions  is  now  estimated  at  39  eV  as  opposed  to  49  eV  predicted  by 
the  sputter  yield  model.  The  E*  for  a  Si  target  being  sputtered  by  Xe  ions  is  estimated  at 
48  eV  as  opposed  to  91  eV  as  predicted  by  the  sputter  yield  model.  These  results  are 
consistent  with  experimental  observations  of  gate  electrode  thresholds  for  sputtering  and 
measured  sputter  yield  values.  Predictions  of  performance  degradation  rates  in  xenon 
environments  correlate  well  with  experimental  results.  Using  the  improved  estimates  for 
Efc,  this  model  can  now  be  used  to  estimate  the  performance  degradation  rates  of  any  Mo 
or  Si  FEA  cathode  operating  in  neutral  xenon  environments  at  any  pressure.  Operating 
voltages  are  limited  to  values  near  E^;  the  operating  voltage  limits  depend  on  the  cathode 
configuration.  This  model  was  also  used  to  determine  the  effect  of  the  CEX  ion 
population  on  the  cathode  performance  degradation  rate.  It  was  shown  that  this  CEX  ion 
population  is  so  detrimental  to  the  cathode  performance  that  upper  limits  on  operating 
voltage  were  further  reduced  by  20  V  when  this  ion  population  is  bombarding  the  field 
emission  microtips.  In  this  situation,  the  population  of  ions  generated  locally  is 
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negligible;  CEX  ions  dictate  performance  degradation  rates.  To  meet  cathode  lifetime 
requirements  of  greater  than  100  hours  in  2x10 5  Torr  of  xenon,  Si  cathodes  will  be 
limited  to  operating  voltages  below  30  V  and  Mo  cathodes  will  be  limited  to  operating 
voltages  below  20  V.  These  values  will  increase  with  decreasing  pressure  from  2xl0'5 
Torr.  Even  it  these  voltages,  current  limiting  cathode  architectures  will  be  required  to 
prevent  catastrophic  arcing  between  the  tips  and  gate  electrodes. 

Modeling  of  space-charge  current  limits  in  the  applications  of  interest  showed  that 
this  issue  could  significantly  limit  the  performance  of  FE  cathodes  in  thruster  and  tether 
environments.  It  was  shown  that  cathode  dimensions  and  environment  play  a  major  role 
in  determining  their  performance  limitations  because  they  affect  the  configuration  of  the 
sheath.  Planar  sheath  modeling  results  showed  that  cathodes  emitting  electrons  into  the 
discharge  chamber  of  an  ion  thruster  will  have  a  lower  current  density  limit  which  is  >8 
A/cm2,  and  exceeds  the  current  density  requirements  in  this  application.  If  positioned  in 
the  center  of  the  PEPL-70  Hall  thruster,  a  field  emission  cathode  emitting  electrons  at  30 
eV  will  have  a  space-charge  current  density  limit  greater  than  ~1.7  A/cm2,  independent  of 
cathode  dimensions,  as  determined  by  the  planar  model.  A  1.2  x  1.2  cm  cathode  in  the 
external  portion  of  a  Hall  thruster  or  ion  engine  emitting  electrons  at  30  eV  should  be 
able  to  emit  near  -68  mA/cm2,  assuming  a  1/r  electron  beam  expansion  in  the  cathode 
sheath.  Either  a  few  cathode  arrays  could  be  used  in  the  center  of  the  thruster  to  provide 
up  to  4.5  A  or  a  segmented  cathode  ring  could  be  used  in  the  external  cathode  position  of 
a  Hall  thruster  with  a  much  larger  area.  The  external  cathode  position  will  be  required 
for  small  and  microscale  propulsion  systems.  Current  requirements  for  small  and 
microscale  propulsion  systems  will  be  less  than  a  few  hundred  milliamperes  so  that  a 
cathode  area  of  a  few  cm2  will  suffice.  If  electrons  are  emitted  at  100  eV,  the  lower  limit 
on  current  density  for  the  1.2  x  1.2  cm  cathode  with  Vg  at  100  V  and  <(>c  at  20  V  should  be 

-120  mA/cm2.  Of  course,  operating  voltages  limit  cathode  lifetime.  While  silicon 
cathodes  are  limited  to  operation  with  gate  electrode  voltages  no  higher  than  30  V, 
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molybdenum  cathode  gate  electrodes  are  limited  to  potentials  no  higher  than  19  V.  If 
electrons  must  be  emitted  at  higher  energies,  the  electron  energy  must  be  decoupled  from 
the  potential  of  the  gate  electrode.  CLAIR  would  provide  additional  electrodes  so  that 
electron  acceleration  up  to  100  eV  or  higher  would  occur  downstream  of  the  gate 
electrode. 

In  a  tether  environment,  much  higher  cathode  operating  voltages  will  be  tolerated 
because  of  low  pressures;  however,  there  are  fewer  ions  to  neutralize  the  negative  space- 
charge  of  the  electron  beam  so  that  the  space-charge  current  limits  will  be  significantly 
lower  than  the  limits  in  the  thruster  environment.  The  current  density  limits  in  this 
environment  are  significantly  affected  by  the  size  of  the  cathode.  Only  a  few  operating 

configurations  were  explored.  A  6.6  x  6.6  mm  cathode  should  not  be  able  to  emit  more 
than  5.6  mA  if  <|>c  is  20  V  and  Vg  is  30  V,  according  to  the  spherical  sheath  model.  If  a 

current  of  2  A  is  required,  then  several  hundred  cathodes  will  be  necessary  to  supply  that 
current.  This  configuration  of  cathodes  may  be  prohibitively  complicated.  With  the  gate 

electrode  at  100-1000  V,  it  is  possible  that  fewer  than  100  cathodes  will  be  required  to 
emit  2  A  in  some  operating  configurations.  It  was  also  shown  that  increasing  <>c  also 

increases  the  emission  current  from  i  tether  system. 

Cathode  dimensions  for  the  applications  discussed  must  consider  both  space- 
charge  limitations  and  cathode  performance  at  the  voltages  limited  by  lifetime 
requirements.  Even  though  current  densities  >1.7  A/cm2  are  possible  in  the  Hall  thruster 
internal  position,  that  performance  is  very  optimistic  for  a  cathode  that  must  survive  for 
>100  hours.  If  the  cathode  microtips  are  not  shielded  from  CEX  ions,  then  Si  and  Mo 
cathodes  will  be  limited  to  operating  voltages  ~30  and  ~19  V  respectively.  Such  low 
operating  voltages  severely  limit  the  current/tip  so  that  the  array  of  microtips  could  be 
prohibitively  large  to  provide  the  required  current.  The  size  of  the  cathode  and  number  of 

microtips  required  in  the  different  applications  depends  on  the  packing  density  of  the 
microtips,  <|)w,  rt,  rg,  and  As.  Results  presented  in  Chapter  4  show  how  the  cathode 
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parameters  affect  the  number  of  tips  required  for  the  desired  current  and  operating 
voltages  limited  by  cathode  lifetime.  Some  possible  cathode  configurations  are  suggested 
in  the  next  section  for  the  applications  of  interest. 

6.2  Recommended  Cathode  Designs 

The  number  of  microtips  or  thin  film  cathodes  necessary  to  emit  the  required 
current  depends  on  the  environment  and  cathode  configuration  and  materials.  The 
environment  places  lower  limits  on  operating  voltages,  as  discussed  in  Chapters  3  and  4, 
and  current  densities,  as  discussed  in  Chapter  5.  The  cathode  materials  and  dimensions 
also  affect  the  number  of  cathodes  required  and  sizes  of  the  arrays  as  discussed  in 
Chapter  4.  The  design  of  FE  cathodes  for  the  applications  discussed  requires  careful 
consideration  of  the  environment,  current  and  lifetime  requirements  for  the  applications, 
and  performance  of  available  cathodes.  The  modeling  tools  developed  for  this 
dissertation  can  be  used  to  optimize  the  cathode  configuration  and  electrode  potentials  for 
the  application. 

The  models  developed  were  used  to  suggest  cathode  configurations  that  would 
satisfy  the  requirements  of  the  applications  discussed.  The  ’Imits  on  current  density  were 
determined  from  the  sheath  model.  Then  data  presented  in  Table  4.5  and  Table  4.6  were 
used  to  determine  which  cathodes  could  supply  that  current  density  at  sufficiently  low 
voltages  to  ensure  cathode  lifetimes  greater  than  100  hours  without  CLAIR  to  retard  CEX 

ions.  One-half  of  the  side  length  of  a  square  cathode  is  given  as  rc.  The  cathode  packing 
density  was  assumed  to  be  108/cm2.  Other  cathode  parameters  were  assumed  to  be  As=2, 

rg=1000  A,  and  rt=40  A.  This  cathode  configuration  is  not  overly  optimistic.  Typically, 

cathode  parameters  are  not  this  good,  however,  higher  packing  densities  and  lower  r„  rg, 
and  As  have  been  demonstrated.  Cathode  configurations  suggested  with  these  assumed 

parameters  are  shown  in  Table  6.1.  Cathode  configurations  are  also  suggested  at  higher 
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operating  voltages  assuming  a  CLAIR  structure  is  used  to  improve  cathode  lifetime  at 
higher  voltages.  Cathodes  with  higher  packing  densities  and  smaller  As,  r„  and  rg  will 

require  smaller  cathode  arrays  or  fewer  arrays  to  supply  the  required  current  than 
configurations  suggested  in  Table  6.1  and  Table  6.2.  The  number  of  cathodes  can  be 
scaled  down  for  much  smaller  thrusters  to  meet  lower  current  requirements.  Modeling 
results  suggest  that  Mo  cathodes  could  not  be  used  without  a  CLAIR  structure,  with  the 
assumed  parameters,  because  they  could  only  supply  the  required  current  from  reasonable 
cathode  areas  with  a  gate  electrode  voltage  that  is  at  least  at  42  V.  In  the  cases  where  Vg 
is  30  V,  Si,  HfC,  or  ZrC  cathodes  could  provide  the  currents  required  from  the  maximum 
cathode  size  noted  in  Table  6.1. 


Table  6.1  Recommended  number  of  cathodes,  operating  voltage,  cathode  material, 
and  emitting  area  required  in  different  applications  if  packing  density  is  108/cm2, 
As=  2,  rg=1000,  and  rt=40. 


Application 

I 

(A) 

rc 

(cm) 

vg 

(V) 

I/cathode 

(mA) 

N 

cathodes 

Area 

(cm2) 

Material 

PEPL-70-ext. 

4.5 

0.6 

30 

>68 

<66 

<95 

Si, HfC, ZrC 

PEPL-70-int. 

4.5 

1.5 

30 

>4500 

1 

<9 

Si,HfC,ZrC 

Ion  engine-int. 

<6 

0.5 

30 

>6000 

1 

<1 

HfC,ZrC 

Tether 

2 

0.33 

30 

<4.3 

>460 

>460 

Mo, Si, HfC, ZrC 

Tether 

2 

0.33 

100 

<9.1 

>220 

>220 

Mo, Si, HfC, ZrC 

Table  6.2  Recommended  number  of  cathodes,  operating  voltage,  cathode  material, 
and  emitting  area  required  in  different  applications  if  packing  density  is  108/cm2, 
As=  2,  rg=1000,  rt=40,  and  CLAIR  is  used. 


Application 

I 

(A) 

rc 

(cm) 

vg 

(V) 

I/cathode 

(mA) 

N 

cathodes 

Area 

(cm2) 

Material 

PEPL-70-ext. 

4.5 

0.6 

100 

>136.0 

<33 

<47.5 

Mo,Si,HfC,ZrC 

PEPL-70-int. 

4.5 

1.5 

42 

>4500 

1 

<9 

Mo, Si, HfC, ZrC 

Tether 

2 

0.33 

100 

<9.1 

>220 

>220 

Mo,Si,HfC,ZrC 

Tether 

2 

0.33 

1000 

<41.0 

>49 

>21 

Mo,Si,HfC,ZrC 
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6.3  Recommendations  for  Future  Research 


This  research  topic  is  rich  and  has  only  begun  to  be  explored.  This  thesis  is  the 
first  in  a  series  of  studies  investigating  the  compatibility  of  FE  cathodes  and  EP  systems 
and  tethers.  Many  of  the  experiments  were  preliminary.  More  of  the  same  experiments 
at  low  operating  voltages  are  recommended  to  support  the  E*  estimates.  The  effect  of 
anode  potential  on  performance  degradation  rates  with  low  gate  electrode  potentials  in 
higher-pressure  environments  should  be  further  investigated.  HfC  and  ZrC  cathodes 
should  be  tested  in  xenon  environments  and  E*  should  be  determined  for  these  materials 
also. 

Cathode  configuration  improvements  have  been  prescribed  to  successfully  merge 
these  technologies.  Cathodes  should  have  smaller  dimensions  and  more  microtips  in  the 
cathode  arrays  with  better  uniformity  than  those  cathodes  tested  in  this  study.  The  next 
step  is  to  fabricate  the  recommended  Mo,  Si,  HfC,  and  ZrC  FEA  cathodes.  Cathodes 
with  VECTL  architectures  are  required  to  prevent  arcing  between  tips  and  gate  electrodes 
and  maximize  cathode  efficiency.  These  cathodes  should  be  tested  at  current  levels  ~1 
mA  in  xenon  environments,  and  life  tests  should  be  conducted  to  show  which  field 
emission  cathodes  can  operate  in  2xl0'5  T  for  -100  hours  and  validate  theoretical  results. 
I-V  traces  in  UHV  environment,  cathode  geometry,  and  work  function  are  now  sufficient 
to  estimate  cathode  lifetime  in  xenon  environments. 

Carbide  cathodes  are  incredibly  robust;  however,  they  operate  at  high  voltages 
and  low  efficiencies.  The  gate  electrode  and  carbon  film  should  be  microfabricated  as  an 
integrated  structure  to  significantly  reduce  operating  voltages  and  current  collected  by  the 
gate  electrode. 

CLAIR  should  be  developed  to  improve  the  compatibility  of  FEA  cathodes  and 
EP  systems.  The  multi-electrode  structure  must  be  microfabricated  and  tested  at  the 
necessary  operating  voltages.  Microtips  must  then  be  deposited  in  the  back  of  the 
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channels  in  CLAIR.  CLAIR  should  be  tested  in  UHV  environments,  xenon 
environments,  and  then  with  EP  systems. 

Field  emission  cathodes  must  be  tested  with  EP  systems.  Initially  they  will  share 
current  with  hollow  cathodes  to  evaluate  lifetime  and  space-charge  limitations  in  a 
thruster  env:  onment.  Eventually,  these  cathodes  could  provide  the  total  discharge 
current.  New  start-up  sequences  must  be  developed  for  FE  cathodes  and  Hall  and  ion 
thrusters.  In  the  typical  start-up  sequences,  the  propellant  breakdown  into  a  plasma  is 
initiated  at  the  cathode.  This  start-up  could  destroy  the  delicate  cathodes.  If  field 
emission  cathodes  are  used,  propellant  breakdown  should  be  initiated  at  the  anode. 

Another  issue  to  explore  is  the  deposition  of  vacuum  chamber  material  on  the 
cathodes.  High  energy  ion  beams  emitted  from  the  thruster  sputter  material  off  of  the 
vacuum  chamber  walls,  which  could  then  redeposit  on  the  cathode,  shorting  the  tips  to 
the  gate  electrodes.  Cathode  positioning  could  play  an  important  role  in  the  effect  of  this 
issue  on  cathode  lifetime. 

Field  emission  cathodes  must  also  be  tested  in  simulated  LEO  environments. 
Modeling  results  of  space-charge  current  limits  must  be  validated  experimentally.  The 
modeling  results  predict  very  large  cathode  sheaths  in  LEO  environments.  Plasma 
potential  measurements  of  space-charge  limited  currents  in  a  simulated  LEO  environment 
will  determine  the  accuracy  of  the  sheath  modeling  results. 

6.4  Other  Applications 

There  are  several  other  applications  of  the  results  of  this  research  including  those 
discussed  in  the  introduction  and  throughout  this  thesis.  These  cathodes  could  be  used 
for  spacecraft  discharging  to  control  the  potential  of  different  surfaces.88  FE  cathodes 
could  also  be  used  with  Field  Emission  Electric  Propulsion  (FEEP)  systems  to  neutralize 
cesium  ion  beams.  Because  cesium  ions  are  accelerated  though  several  keV,  more  robust 
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cathode  configurations  are  recommended,  like  the  carbon  film  cathodes.  It  has  been 
proposed  to  use  FE  cathodes  as  pressure  sensors.89  The  models  developed  in  this  thesis 
can  be  used  to  design  these  sensors  and  estimate  lifetime.  The  indirect  methods  used  to 
determine  E*  for  sputtering  of  Mo  and  Si  by  xenon  ions  could  be  applied  to 
experimentally  determine  these  vah«es  for  other  FEAC  materials. 
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APPENDIX  A 


NOMENCLATURE 


Fowler-Nordheim  term  (e(VA)2/fs) 

Afn  Fowler-Nordheim  I  term  (A/V2) 

a0  Bohr  radius  (A) 

b  term  in  radius  distribution  factor  for  cathode 

b,^  area  factor  (A2) 

bfc  Fowler-Nordheim  term  (eV/A) 

Bjn  Fowler-Nordheim  I  term  (V) 

c  speed  of  light  (A/fs) 

c0  term  in  radius  distribution  factor  for  cathode 
d  atomic  spacing  (nm) 

e  electron  charge  (C) 

E  ion  energy  (eV) 

Efo  energy  threshold  for  sputtering  (eV) 

Ftip  electric  field  on  axis  of  the  microtip  structure  (V/cm) 
h  Plank’s  constant  (eV-fs) 

Ia  current  collected  by  the  anode  (A) 

Id  Hall  thruster  discharge  current  (A) 

Ie  field  emission  current  (A) 

IM  cathode  current  (A) 

If  final  cathode  current  after  100  hour  xenon  exposure  (A) 

Ig  current  collected  by  the  gate  electrode  (A) 

I0  initial  cathode  current  (A) 
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Itip  current  emitted  from  a  single  microtip  structure  (A) 

Imns  current  transmitted  through  the  sheath  (A) 

J  ion  current  density  measured  by  the  Faraday  probe  (A/cm2) 
je  plasma  electron  current  density  (mA/cm2) 
jM  electron  beam  current  density  (mA/cm2) 

T  normalized  electron  beam  current  density 

ec 

Jeemax  normalized  upper  limit  on  electron  beam  current  density 
Jm  Fowler-Nordheim  field  emission  electron  current  density  (A/cm2) 
k  Boltzmann  constant  (J/K) 

k  electric  field  factor 

K  term  in  sputter  yield  model 

m  mass  flow  rate  of  propellant  (mg/s) 

M,  incident  ion  atomic  weight 

Mn  target  atom  atomic  weight 

Mj  initial  mass  of  a  rocket  with  fuel  (kg) 

Mf  final  mass  of  a  rocket  (kg) 

n  number  of  ions  striking  the  emitting  area  of  a  microtip  during  one  second 

ne  plasma  electron  number  density  (electrons/cm3) 

ntt  electron  beam  electron  number  density  (electrons/cm3) 

neo  electron  number  density  in  the  plasma  (electrons/cm3) 

rij  ion  number  density  (ions/cm3) 

ni0  ion  number  density  at  the  sheath  boundary  (ions/cm3) 

ns  number  of  atoms  sputtered  from  the  emitting  area  of  a  microtip  per  second 

N  neutral  xenon  number  density  (atoms/cm3) 

N, jps  number  of  tips  in  a  field  emission  array  cathode 

xenon  pressure  measured  by  the  NPF  probe  (Torr) 
xenon  pressure  measured  by  the  vacuum  gauge  (Torr) 


133 


Qi  cross-section  for  ionization  (7ta02) 

Q  term  in  sputter  yield  model 

r  radial  distance  from  microtip  structure  (cm) 
r.  radial  position  corresponding  to  ionization  potential  (cm) 
rc  cylindrical  or  spherical  cathode  radius  (cm) 
r,  radius  of  curvature  of  a  microtip  structure  (cm) 
rg  radius  of  aperture  in  gate  electrode  (cm) 

rm  maximum  distance  from  a  microtip  within  which  an  ion  formed  will  hit 
the  emitting  area  of  the  tip  (cm) 

rs  radial  position  corresponding  to  the  sputter  threshold  voltage  (cm) 
sn  term  used  in  sputter  yield  model 

Se  term  used  in  sputter  yield  model 

T  temperature  (eV) 

te  time  to  erode  a  single  layer  of  atoms  from  a  microtip  emitting  area  (s) 

Te  electron  temperature  (eV) 

Tn  thrust  (N) 

U.  sublimation  energy  (eV) 

Va  anode  voltage  with  respect  to  ground 

Vd  discharge  voltage  of  the  Hall  thruster  (V) 
vM  velocity  of  electrons  emitted  from  the  cathode  (m/s) 

Vg  potential  of  the  gate  electrode  with  respect  to  the  cathode  emitting 
surface  (V) 

Vj  ion  velocity  (m/s) 

Vj  ionization  potential  (eV) 

VM  modulation  voltage  (V) 

maximum  voltage  in  an  I-V  trace  (V) 
minimum  voltage  in  an  I-V  trace  (V) 
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Vpk  peak  voltage  (V) 

Vs  voltage  corresponding  to  energy  threshold  for  sputtering  (V) 

V,  ion  temperature  (eV) 

x  position  from  the  cathode  surface  (m) 
x„  term  in  Bw 

Y  sputter  yield  (atoms/ion) 

a  term  in  sputter  yield  model 

pc  microtip  cone  angle  (radians) 

Pg  field  enhancement  factor 

5  term  in  Bj*, 

As  tip  radius  spread  factor 

E  term  in  sputter  yield  model 

e0  permittivity  of  free  space  (F/m) 

T)  normalized  voltage 

T|c  normalized  potential  of  gate  electrode  below  the  potential  at  the  sheath 
boundary 

r|eff  thruster  efficiency 

rig  normalized  potential  between  the  gate  electrode  and  cathode  potential 
r|0  normalized  ion  accelerating  potential  of  the  pre-sheath 
lD  electronic  Debye  length  (cm) 

o  half  cone  emitting  angle  of  a  cathode  tip  (°) 

S  distribution  factor  associated  with  a  spread  in  tip  radii  in  a  field  emission 
cathode  array 

<]>  potential  with  respect  to  the  potential  at  the  sheath  boundary  (V) 

<j)0  accelerating  potential  of  the  sheath  (V) 
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potential  of  the  gate  electrode  with  respect  to  the  potential  at  the  sheath 
boundary  (V) 

electron  emission  half  cone  angle  on  a  microtip  (radians) 
work  function  (eV) 
normalized  position 

normalized  spherical  or  cylindrical  cathode  radius 
normalized  position  of  sheath  boundary 
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APPENDIX  B 


CLAIR 

A  Cathode  Lens  and  Ion  Repeller  (CLAIR)  has  been  designed  to  significantly 
improve  the  compatibility  of  Hall  thrusters  and  FEA  cathodes.  The  performance  of  this 
structure  as  a  lens  and  ion  filter  was  demonstrated  using  a  particle-in-cell  (PIC)  code. 
CLAIR  is  a  series  of  electrodes,  similar  to  a  gate  electrode  in  FEA  cathodes,  in  an  einzel 
lens  configuration  that  accelerates  or  decelerates  an  electron  beam  and  retards  charge- 
exchange  (CEX)  ion  flow.  The  CLAIR  configuration  is  shown  in  Figure  B.l.  Three 
electrodes  are  used  in  addition  to  the  gate  electrode.  Electrode  V2  is  at  a  high  potential, 
-100  V,  to  retard  the  flow  of  CEX  ions  to  the  microtip  structures.  V3  is  approximately  20 
V  below  plasma  potential.  The  other  electrodes  are  biased  with  respect  to  V3,  depending 
on  the  cathode  operating  voltage  and  the  accelerating  or  decelerating  function  of  CLAIR. 
The  thickness  of  the  electrodes  and  insulators  determines  the  lens  optics.  It  is  more 
difficult  to  focus  the  electron  beam  when  it  is  being  decelerated  by  CLAIR  than  when  it 
is  being  accelerated.  One  possible  cathode  configuration  is  represented  by  the  parameters 
shown  in  Table  B.l.  The  design  and  performance  of  CLAIR  was  evaluated  using 
MAGIC,  a  PIC  code,  that  determined  equipotential  configurations  and  electron 
trajectories.  MAGIC  was  developed  by  Mission  Research,  Inc.  MAGIC  output  is  shown 
Figures  B.2  and  B.3.  The  cathode  was  assumed  to  be  a  planar  surface  with  electrons 
emitted  uniformly  from  it  with  a  velocity  vector  directed  along  the  axis  of  CLAIR  only. 
This  cathode  configuration  should  be  capable  of  filtering  out  60  eV  ions  and  focusing  the 
electron  beam  through  the  electrodes,  as  verified  using  MAGIC. 


137 


~1.2  jim 


Figure  B.l  Cathode  Lens  and  Ion  Repeller  (CLAIR)  congifuration. 


Table  B.l  Electrode  potentials,  thicknesses,  and  spacing  for  the  CLAIR 
configuration  simulated  by  MAGIC. 


Vc 

Vg 

V, 

V2 

V3 
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*g 
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CLAIR  can  be  used  to  significantly  increase  the  space-charge  current  limit  on  the 
cathode  by  increasing  the  energy  of  the  electrons  emitted  from  the  cathode  without 
increasing  the  potential  of  the  gate  electrode.  CLAIR  essentially  decouples  the  electron 
energy  from  the  gate  electrode  potential.  It  can  significantly  increase  the  lifetime  of  the 
cathode  by  filtering  the  CEX  ions  from  the  thruster  out  of  the  flux  of  particles 
bombarding  the  microtips.  Without  CLAIR,  cathode  gate  electrode  operating  voltages 
are  limited  to  values  near  30  V.  Ions  generated  between  the  gate  electrodes  and  microtips 
and  CEX  ions,  which  sputter  the  microtip  structures,  are  responsible  for  limiting  the 
cathode  lifetime.  If  CLAIR  is  used,  operating  voltages  can  be  increased  by 
approximately  20  V  from  the  30  V  limit  for  Si  cathodes.  At  this  higher  operating  voltage, 
the  current  is  increased  by  orders  of  magnitude  from  the  same  array  of  tips.  With  higher 
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voltage  gate  electrodes,  space-charge  current  limits  are  increased  for  a  specific  current 
requirement,  as  shown  in  Figure  4.6. 


Figure  B.3  Electron  trajectories  predicted  by  MAGIC  for  the  CLAIR  configuration 
described  in  Figure  B.l  and  Table  B.l. 
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APPENDIX  C 


FEA  Cathode  Performance  Degradation  Model  Macro 

This  macro  is  used  in  IGOR  to  calculate  the  performance  degradation  rate  of  FEA 
cathodes  in  xenon  environments.  The  model  uses  pre-existing  arrays  for  radius, 
Qiintxenon,  and  xenoneVi.  The  radius  array  has  all  of  the  radial  positions  (cm)  used  in 
the  calculation  of  the  ion  flux  to  the  microtips.  Ionization  cross-section  data  obtained 
from  Brown  are  in  the  arrays  Qiintxenon  (cross-section)  and  xenoneVi  (energy).  The  tim 
array  is  used  to  calculate  the  change  in  radius  in  time  increments  tim  (hours).  Those  data 
arrays  are  given  after  the  macro. 


Macro  calclt  0 
Silent  1 

variable  ag,  phi,  be,  Ntips,  Vg,  a,  b,  n,  elect,  vs,  pressure,  vt,  d,  yeti,  rtmax 
ag= 10000  //(A) 

phi=3.5  //(eV) 

bc=0.26  //(radians) 

Ntips= 16000 

ds(0)=340 

Vg  =86  //(V) 

YEti=.017  //  sputter  yield  value  for  a  Si  target  and  xenon  CEX  ions 

rt(0)  =59.1  //initial  effective  tip  radius  (A) 

rtmax=  (ds(0)+l)*rt(0)  //initial  maximum  tip  radius 

elect  =  1 .6022E- 1 9  //(Coulombs) 

pressure  =  7E-6  //(Torr) 

vt=  .03  //(eV) 

d  =  3.13  //(A) 

variable  Us,  MI,  MO,  ZI,  ZQ,  Que,  eth,  alpha,  kay,  lkay,  i,  ethmf 

//duplicate/O  rt  Irtip,  Irarray,  Ftip,  JFN,  barea,  bb,  SIG,  kc,  bg,  co,  delr,  ns,  te 

duplicate/O  radius  voltage,  Quei,  Sxesi,  epsilon,  se,  sn 

duplicate/O  radius  dr,  nsi,  qydr 

differentiate  dr 

n  =  0 

ns(n)  =0 
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do 


ds(n)  =  rtmax/rt(n)-l 
kc(n)=l/54  *(86+ag/rt(n))*  l/tan(bc) 

Ftip(n)=((pi/(ln(kc*ag/rt(n))))-(tan(bc))A2)*Vg/rt(n) 

a=(.008621/phi)*exp(9.246/phiA.5)  //Fowler-Nordheim  a  coefficient 

b=0.6421 *phiAl  .5  //Fowler-Nordheim  b  coefficient 

JFN(n)=a*Ftip(n)A2*exp(-b/Ftip(n)) 
bg(n)=  ((pi/(ln(kc(n)*ag/rt(n))))-(tan(bc))A2)/rt(n) 

barea(n)=  2*pi*rt(n)A2*(cos(bc))A2*(bg(n)*Vg)/(b+(sin(bc))A2*bg(n)*Vg) 
co(n)=(  1  /(pi  *bg(n)*rt(n)))  *  ((tan(bc))A2+bg(n)*rt(n))A2- 1 
bb(n)=((co(n)*b+(3*co(n)+2)*bg(n)*Vg)/(bg(n)*Vg»- 
(co(n)*bg(n)*Vg/b)*(tan(bc))A2 

SIG(n)=(exp(ds(n)*bb(n))-l)/(ds(n)*bb(n))  //=2 

Irtip(n)  =  barea(n)*JFN(n)*  1 .6E-4  //current  from  a  single  tip  with  rt in  Amperes 

Irarray(n)  =  Ntips*SIG(n)*barea(n)*JFN(n)*  1 .6E-4  //current  in  Amperes 

IrarrayuA  =  Irarray*  le6  //current  in  microamperes 

//calculate  voltage  as  a  function  of  radial  distance  from  the  tips  (r  in  cm) 
rm(n)  =  rt(n)*lE-8*(vg/vt)A.3333 
vm(n)  =  vg*((rt(n)*  lE-8/rm)-  l)*((rt(n)/  ag)- 1  )A- 1 
voltage  =  vg*((rt(n)*  lE-8/radius)-  l)*((rt(n)/ ag)- 1  )A- 1 

//calculate  sputter  yield  for  the  voltage  corresponding  to  each  radial  position 

Us=  4.63 

MI=131.3 

MH=28.1 

Que=  0.78 

71=  54 

ZH=  14 

Ethmf  =  .53  //energy  threshold  multiplication  factor 
Eth=(1.9  +  3.8*(MII/MI)A-1  +  0.134*(MII/MI)A1.24)*Us*Ethmf 
alpha  =  0.08+0. 164*(MH/MI)A0.4+0.0 145*(MII/MI)A  1 .29 
Kay=  8.478*(ZI*ZII/((ZIA.66667  +  ZnA.66667)A0.5))*(MF(MI+Mn)) 
epsilon=(0.03255/(ZI*ZII*(ZIA.66667+ZnA.66667)A0.5))* 
(MH/(MI+MII))*voltage 

sn=  3.44*(epsilonA0.5)*ln(epsilon+2.718)/(l+6.355*(epsilonA0.5)+epsilon* 
(-1.708+6.882*(epsilon)A0.5)) 
lkay=  0.79*(((MI+Mn)A.66667)/(MIA.66667*MnA1.5))* 
(ZIA.66667*ZIIA.5)/(ZIA.66667+ZnA.66667)A.75 
se  =  lkay*epsilonAl/2 

Sxesi  =  0.42*  alpha*Que*Kay*sn*(Us*(l+0.35*Us*se))A-l* 

( 1  -(eth/voltage)  A0. 5)A2 . 8 

//calculate  te 

Interpolate/T=l/N=28/I=3/A=0/J=2/Y=Quei/X=voltage  Qiintxenon  /X=xenoneVi 
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ns(n)=0 

i=0 

do 

if  (voltage(i)  >  vg ) 

voltage(i)  =  vg 

endif 

i=i+l 

while  ( i  <  numpnts(nsi) ) 
i=0 

do 

if  (voltage(i)  >  eth ) 

if  (voltage(i)  <  vm(n)) 

qydr(i)  =dr(i)*0.5*(Sxesi(i)*Quei(i)+Sxesi(i+l)* 
Quei(i+1))*8.79E-17 

nsi(i)  =  (Irtip(n)/elect)*pinterp(n)*3.55E16*qydr(i) 

ns(n)=  ns(n)  +  nsi(i) 

endif 

endif 

i=i+l 

while  ( i  <  numpnts(nsi)  ) 

//  nsti(n)  =  (2e-5/1.6e-19)*2*pi*rt(n)A2*le-16*0.33*Yeti  //CEXions 
nsti(n)  =  0 

te(n)  =  pi/3600*2*cos(37*2*pi/360)*rt(n)A2/((d)A2*(ns(n)+nsti(n))) 
delr(n)  =  tim(n)*d/te(n)  //change  in  tip  radius  during  time  tim. 
rt(n+l)  =  rt(n)  +delr(n) 

n=n+l 

while  (n<numpnts(tim)) 
end  macro 


radius 

xenon  Vi 

Qiintxenon 

tim 

4.25e-07 

12 

0.118 

0.1 

4.5e-07 

13 

0.6 

0.1 

5e-07 

14 

1.12903 

0.1 

5.5e-07 

15 

1.45161 

0.1 

6e-07 

16 

1.77419 

0.1 

6.5e-07 

17 

2.05556 

0.1 

7e-07 

18 

2.24074 

0.1 

7.5e-07 

19 

2.42593 

0.1 

8e-07 

20 

2.64286 

0.1 

8.5e-07 

30 

4.45918 

0.1 

9e-07 

40 

5.10515 

9.5e-07 

50 

5.32886 

le-06 

60 

5.55257 

l.le-06 

70 

5.77629 
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1.2e-06 

80 

6 

1.3e-06 

90 

5.93333 

1.4e-06 

100 

5.86667 

1.5e-06 

200 

5.2 

1.6e-06 

300 

4.8125 

1.7e-06 

500 

4.0375 

1.8e-06 

700 

3.2625 

1.9e-06 

900 

2.4875 

2e-06 

1000 

2.1 

2.1e-C6 

2.2e-06 

2.3e-06 

2.4e-06 

2.5e-06 

2.6e-06 

2.7e-06 

2.8e-06 

2.9e-06 

3e-06 

3.1e-06 

3.2e-06 

3.3e-06 

3.4e-06 

3.5e-06 

3.6e-06 

3.7e-06 

3.8e-06 

3.9e-06 

4e-06 

4.1e-C6 

4.2e-06 

4.3e-06 

4.4e-06 

4.5e-06 

4.6e-06 

4.7e-06 

4.8e-06 

4.9e-06 

5e-06 

5.1e-06 

5.2e-06 

5.3e-06 

5.4e-06 

5.5e-06 

5.6e-06 
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5.7e-06 

5.8e-06 

5.9e-06 

6e-06 

6.1e-06 

6.2e-0 6 

6.3e-06 

6.4e-06 

6.5e-06 

6.6e-06 

6.7e-06 

6.8e-06 

6.9e-06 

7e-06 

7.1e-06 

7.2e-06 

7.3e-06 

7.4e-06 

7.5e-06 

7.6e-06 

7.7e-06 

7.8e-06 

7.9e-06 

8e-06 

8.1e-06 

8.2e-06 

8.3e-06 

8.4e-06 

8.5e-06 

8.6e-06 

8.7e-06 

8.8e-06 

8.9e-06 

9e-06 

9.1e-06 

9.2e-06 

9.3e-06 

9.4e-06 

9.5e-06 

9.6e-06 

9.7e-06 

9.8e-06 

9.9e-06 

le-05 

l.le-05 
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1.2e-05 

1.3e-05 

1.4e-05 

1.5e-05 

1.6e-05 

1.7e-05 

1.8e-05 

1.9e-05 

2e-05 

2.1e-05 

2.2e-05 

2.3e-05 

2.4e-05 

2.5e-05 

2.6e-05 

2.7e-05 

2.8e-05 

2.9e-05 

3e-05 

3.1e-05 

3.2e-05 

3.3e-05 

3.4e-05 

3.5e-05 

3.6e-05 

3.7e-05 

3.8e-05 

3.9e-05 

4e-05 

4.1e-05 

4.2e-05 

4.3e-05 

4.4e-05 

4.5e-05 

4.6e-05 

4.7e-05 

4.8c-05 

4.9e-05 

5e-05 

5.1e-05 

5.2e-05 

5.3e-05 

5.4e-05 

5.5e-05 

5.6e-05 
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5.7e-05 

5.8e-05 

5.9e-05 

6e-05 

6.1e-05 

6.2e-05 

6.3e-05 

6.4e-05 

6.5e-05 

6.6e-05 

6.7e-05 

6.8e-05 

6.9e-05 

7e-05 

7.2e-05 

7.3e-05 

7.4e-05 

7.5e-05 

7.6e-05 

7.7e-05 

7.8e-05 

7.9e-05 

8e-05 

8.1e-05 

8.2e-05 

8.3e-05 

8.4e-05 

8.5e-05 

8.6e-05 

8.7e-05 

8.8e-05 

8.9e-05 

9e-05 

9.1e-05 

9.2e-05 

9.3e-05 

9.4e-05 

9.5e-05 

9.6e-05 

9.7e-05 

9.8e-05 

9.9e-05 

0.0001 
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APPENDIX  D 


Cathode  Sheath  Model  Macro 

This  appendix  presents  IGOR  macros  used  to  calculate  potential  and  electric  field 
profiles  m  cathode  sheaths  assuming  planar,  cylindrical,  or  spherical  cathode  and  sheath 
geometries.  Poisson’s  equation  is  solved  using  a  fourth  order  Runge-Kutta  method. 


Macro  PlanarsheathO 
Silent  1 

variable  h,  n,  f,  etav,  Te,Vb,  Vm,  Vg,  bk,  mi,  me,  elec,  neo,  sqrtalf,  je,  ljip,  ljee,  Jip,  Jee 

variable  etam,  etag,  etap,  etab,  etao,  Vo,  coni,  con2,  con3,  con4,  Vc,  etac,ef,  etabo 

neo  =  4E4 

Te  =  0.086 

elec  =  1.6E-19 

bk  =  1.38E-23 

me  =  9.11E-31 

Vg  =  5 

ljee  =1 

Vc=30 

Uim=140 

h=0.2j  //  step  size  in  spatial  position, 

je=  neo*elec*100*(2*bk*Te*11600/me)A0.5 

etag  =  Vg/Te 
etac=  Vc/Te 
Jee  =  ljee/je 

etao  =  (l+Jee*(etag+etac)A-0.5)/(2-Jee*(etag+etac)A-1.5) 

exIII(O)  =  0  //  electric  field  at  sheath  boundary 

etam(0)=  etac  //potential  at  cathode  surface 


n=0 

do 


etav  =  etain(n) 
exIII(n)  =  exIII(O)  +h*n 


149 


ef  =  4*Jee*((etag+etac-etav)A0.5  -(etag+etac)A0.5)  + 

2*(exp(-etav)-l)+4*etao*(l+Jee*((etag+etac)A-0.5))* 

(((1  +etav/etao)  A0.5)- 1 ) 

f  =  -(4*(Jee*((etag+etac-etav)A0.5  -(etag+etac)A0.5)  + 

0.5*(exp(-etav)-l)+etao*(l+Jee*((etag+etac)A-0.5))* 
(((l+etav/etao)A0.5)-l)))A0.5 
kayini(n)=  h*f 
etav  =  etalLI(n)  +  kayini(n)/2 
f  =-  (4*(Jee*((etag+etac-etav)A0.5  -(etag+etac)A0.5)  + 

0.5*(exp(-etav)-l)+etab*(l+Jee*((etag+etac)A-0.5))* 
(((l+etav/etao)A0.5)-l)))A0.5 
kayIH2(n)=  h*f 
etav  =  etaIH(n)  +  kayIII2(n)/2 

f  -  -(4*(Jee*((etag+etac-etav)A0.5  -(etag+etac)A0.5)  + 

0.5*(exp(-etav)-l)+etao*(l+Jee*((etag+etac)A-0.5))* 
(((l+etav/etao)A0.5)-l)))A0.5 
kayIII3(n)  =  h*f 
etav  =  etalH(n)  +  kayIII3(n) 

f  =  -(4*(Jee*((etag+etac-etav)A0.5  -(etag+etac)A0.5)  + 

0.5  *  (exp(-etav)- 1  )+etao*  ( 1  +Jee*((etag+etac)A-0.5))  * 

((( 1 +etav/etao)  A0.5)- 1  )))A0.5 
kaym4(n)  =  h*f 

etam(n+l)  =  etalll(n)  +  l/6*(kayIHl(n)  +  2*kayIII2(n)  +  2*kayIII3(n)  + 
kayIEI4(n)) 


n=  n+1 
while  (n<llim) 
exIII(n)  =  ex]H(0)  +h*.i 

//calculate  charged  particle  number  densities 
ninom=(l+Jee*((etag+etac)A-0.5))*((l+etam/etao)A-0.5) 

ni=  (l+ljee*((Vg+Vc)A-.5)/(neo*elec*100*(2*elec/me)A.5))*(l+etain/etao)A-0.5 
nenoin=  exp(-etalll) 
ne=  exp(-etalll) 

neenoIH  =  Jee*(etag-etain+etac)A-0.5 

nee=  ljee/neo*((Vg+Vc-etain*Te)A-0.5)/(elec*  100*(2*elec/me)A0.5) 
netotin=  nenoIH+neenoin 
ntotm  =  ninoIH  -  netotm 
ntot  =  ni-ne-nee 

efieldm  =  (4*(Jee*((etag+etac-etaIII)A0.5  -(etag+etac)A0.5)  + 
0.5*(exp(-etaIII)-l)+etao*(l+Jee*((etag+etac)A-0.5))* 
(((l+etam/etab)A0.5)-l)))A0.5 

efieldi=  4*etao*(l+Jee*(etag+etac)A-0.5)*(((l+etaIII/etao)A0.5)-l) 
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efielde=  2*(exp(-etaIII)-l) 

efieldee=  4*Jee*((etag+etac-etam)A0.5  -(etag+etac)A0.5) 
end  macro 


Macro  CylindricalsheathO 
Silent  1 

variable  h,  n,  zrsh,  llim,  etac,  etao,  etag,  ef,  efl,  ef2,  Jee,  zrc,  etav,  efld,  rad 
h=  0.2  //  step  size  in  spatial  position,  £ 

zrsh=13  //normalized  sheath  radius  ,rsb/XD 

zrc=4  //normalized  cathode  radius,  r <Ad 

llim=45 

etac=4 

etag=6 

Jee=6.6 

Make/N=46/D/0  eta,  efield,  rr 
eta(0)  =  0 

efield(O)  =  0  //efield  at  sheath  boundary 

rr(0)=  zrsh  //position  of  sheath  boundary 


etao=(l+Jee*(etag+etac)A-0.5*(zrc/zrsh))/(2-Jee*(etag+etac)A-1.5*(zrc/zrsh)) 

print  "etao=",  etao 

n=0 

do 

etav=  eta(n) 
efld=efield(n) 
rad=  rr(n) 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh))*(l+etav/etao)A-0.5* 

(zrsh/rad)-Jee*(etag+etac-etav)A-0.5*(zrc/rad)-exp(-etav) 
efl=efld 
ef2=ef-efl/rad 
kayl(l)=  h*efl 
kay2(l)=  h*ef2 

etav=  eta(n)+0.5*kayl(l) 
efld=efield(n)+0 . 5  *kay  2(  1 ) 
rad=  rr(n)-h/2 

ef=  ( 1  +  Jee  *  (etag+etac )  A-0 . 5  *  (zrc/zrsh))*  ( 1  +etav/ etao)  A-0 . 5  * 

(zrsh/rad)-Jee*(etag+etac-etav)A-0.5*(zrc/rad)-exp(-etav) 
efl=efld 
ef2=ef-efl/rad 
kayl(2)=  h*efl 
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kay2(2)=  h*ef2 


etav=  eta(n)+0.5*kayl(2) 
efld=efield(n)+0.5*kay2(2) 
rad=  rr(n)-h/2 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh))*(l+etav/etao)A-0.5* 

(zrsh/rad)-Jee*(etag+etac-etav)A-0.5*(zrc/rad)-exp(-etav) 
efl=efld 
ef2=ef-efl/rad 
kayl(3)=  h*efl 
kay2(3)=  h*ef2 

etav=  eta(n)+0.5*kayl(3) 
efld=efield(n)+0.5*kay2(3) 
rad=  rr(n)-h 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh))*(l+etav/etao)A-0.5* 

(zrsh/rad)-Jee*(etag+etac-etav)A-0.5*(zrc/rad)-exp(-etav) 
efl=efld 
ef2=ef-efl/rad 
kayl(4)=  h*efl 
kay2(4)=  h*ef2 

n=n+l 

eta(n)  =  eta(n-l)  +( l/6)*(kay  1  (1  )+2*kay  1  (2)+2*kay  1  (3)+kay  1  (4)) 
efield(n)  =  efield(n-l)  +(l/6)*(kay2(l)+2*kay2(2)+2*kay2(3)+kay2(4)) 
rr(n)  =  rr(n-l)-h 

while  (ncllim) 

//calculate  charged  particle  number  densities  in  the  sheath 
Duplicate/O  efield  ntot,  nion,npelect,nbelect 
ntot  =  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh))*(l+eta/etao)A-0.5* 
(zrsh/rr)-Jee*(etag+etac-eta)A-0.5*(zrc/rr)-exp(-eta) 
nion=  (l+Jee*((etag+etac)A-0.5)*(zrc/zrsh))*(l+eta/etao)A-0.5*(zrsh/rr) 
npelect=Jee*(etag+etac-eta)A-0.5*(zrc/rr) 
nbelect=exp(-eta) 
end  macro 


Macro  SphericalsheathO 
Silent  1 

variable  h,  n,  zrsh,  Him,  etac,  etao,  etag,  ef,  efl,  ef2,  Jee,  zrc,  etav,  efld,  rad 
h=  .2 
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zrsh=10.4 

zrc=4 

llim=32 

etac=4 

etag=6 

Jee=17.6 

Make/N=33/D/0  eta,  efield,  rr 

eta(O)  =  0 
efield(O)  =  0 
rr(0)=  zrsh 

etao=(l+Jee*(etag+etac)A-0.5*(zrc/zrsh)A2)/(2-Jee*(etag+etac)A-1.5*(zrc/zrsh)A2) 


n=0 

do 

etav=  eta(n) 
efld=efield(n) 
rad=  rr(n) 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh)A2)*(l+etav/etao)A-0.5* 

(zrsh/rad)A2-Jee*(etag+etac-etav)A-0.5*(zrc/rad)A2-exp(-etav) 
efl=efld 
ef2=ef-2*efl/rad 
kayl(l)=  h*efl 
kay2(l)=  h*ef2 

etav=  eta(n)+0.5*kayl(l) 
efld=efield(n)+0.5*kay2(l) 
rad=  rr(n)-h/2 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh)A2)*(l+etav/etao)A-0.5* 

(zrsh/rad)A2-Jee*(etag+etac-etav)A-0.5*(zrc/rad)A2-exp(-etav) 
efl=efld 
ef2=ef-2*efl/rad 
kayl(2)=  h*efl 
kay2(2)=  h*ef2 

etav=  eta(n)+0.5*kayl(2) 
efld=efield(n)+0.5*kay2(2) 
rad=  rr(n)-h/2 

ef=  ( 1  +  Jee  *  (etag+etac )  A-0 . 5  *  (zrc/zrsh)  A2)  *  ( 1  +etav/etao)  A-0. 5  * 

(zrsh/rad)A2-Jee!|c(etag+etac-etav)A-0.5*(zrc/rad)A2-exp(-etav) 
efl=efld 
ef2=ef-2*efl/rad 
kayl(3)=  h*efl 
kay2(3)=  h*ef2 


153 


etav=  eta(n)+0.5*kayl(3) 
efld=efield(n)+0.5*kay2(3) 
rad=  rr(n)-h 

ef=  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh)A2)*(l+etav/etao)A-0.5* 

(zrsh/rad)A2-Jee*(etag+etac-etav)A-0.5*(zrc/radA2)-exp(-etav) 

efl=efld 

ef2=ef-2*efl/rad 
kayl(4)=  h*efl 
kay2(4)=  h*ef2 

n=n+l 

eta(n)  =  eta(n-l)  +(l/6)*(kayl(l)+2*kayl(2)+2*kayl(3)+kayl(4)) 
efield(n)  =  efield(n-l)  +(l/6)*(kay2(l)+2*kay2(2)+2*kay2(3)+kay2(4)) 
rr(n)  =  rr(n-l)-h 

while  (n<llim) 

//calculate  charged  particle  number  densities 

Duplicate/O  efield  ntot,  nion,npelect,nbelect 
ntot  =  (l+Jee*(etag+etac)A-0.5*(zrc/zrsh)A2)*(l+eta/etao)A-0.5* 
(zrsh/rr)A2-Jee*(etag+etac-eta)A-0.5*(zrc/rr)A2-exp(-eta) 
nion=  (l+Jee*((etag+etac)A-0.5)*(zrc/zrsh)A2)*(l+eta/etao)A-0.5*(zrsh/rr)A2 
npelect= Jee*  (etag+etac-eta)A-0 . 5  *  (zrc/rr)  A2 
nbelect=exp(-eta) 

end  macro 
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